Introduction a la geodynamique 
Notes pour le cours SCT4003 
version beta January 12, 2011 


Jean-Claude Mareschal 


exosup.com 


page facebook 


exosup.com 


page facebook 


Contents 


Introduction page 1 

References 3 

1 Historique 5 

1.0. 1 Revolution scientifique et paradigme. 5 

1.0. 2 Les grands debats en Geologie 6 

1.0. 3 La contribution de Wegener: la derive des 

continents 7 

1.0. 4 Reponse a l’hypothese de Wegener. 9 

1.0. 5 Exploration des fonds marins. Paleomagnetisme. 10 

1.1 Remarques 11 

2 Structure interne de la Terre 13 

2.1 Structure sismique de la Terre 13 

2.2 Isostasie 14 

2.3 Structure rlreologique 18 

2.3.1 Rheologie 19 

2.3.2 Deformation continue des roches 19 

2.3.3 Loi de Byerlee 20 

3 Expansion des fonds oceaniques 23 

3.1 Fonds oceaniques 23 

3.1.1 Morphologie des fonds oceaniques 23 

3.1.2 Structure sismique de la croute oceanique 24 

3.1.3 Gravite 25 

3.1.4 Flux de chaleur 25 

3.1.5 Geologie marine 26 

3.2 Les inversions du champ magnetique terrestre. 26 

3.2.1 Anomalie magnetique 26 


exosup.com 


page facebook 


IV 


Contents 


3.2.2 Proprietes magnetiques des roches et 

mineraux 27 

3.2.3 Magnetisme inverse dans les roches 28 

3.2.4 Chronometrie des inversions 28 

3.2.5 Mecanisme des inversions 30 

3.3 Anomalies magnetiques marines: Hypothese de 

Vine-Matthews-Morley. 31 

3.4 Applications 33 

3.5 Problemes en suspens 35 

3.6 Notes supplementaires 37 

3.6.1 1: La methode sismique refraction 37 

3.6.2 Notes supplementaires: Flux de chaleur 39 

3.6.3 Notes supplementaires: Anomalies du 

champ magnetique 40 

3.6.4 Champ d’un dipole. 40 

3.6.5 Proprietes magnetiques des roches et des 

mineraux. 41 

3.7 Notes supplementaires: Chronometrie K-Ar. 41 

4 Les plaques: Cinematique 44 

4.1 Failles transformantes 44 

4.2 Zones de subduction 46 

4.3 Tectonique des plaques 46 

4.4 Cinematique des plaques. 51 

4.4.1 Description du mouvement dans l’espace 

des vitesses 51 

4.4.2 Mouvement sur une sphere. Poles d’Euler 52 

4.4.3 Lois de conservation 53 

4.4.4 Stability des points triples 53 

4.5 Applications 54 

4.5.1 Tectonique de la cote Pacifique de l’Amerique 

du Nord 54 

5 Point Chauds-Cycle de Wilson 61 

5.1 Points Chauds 61 

5.1.1 Points chauds actifs 63 

5.1.2 Fixite des points chauds 63 

5.1.3 Application: Ouverture de l’Atlantique 63 

5.2 Cycle de Wilson 65 

6 Sismologie et tectonique des plaques 66 

6.1 Observations des tremblements de Terre 66 


exosup.com 


page facebook 


Contents v 

6.2 Sismicite et definition des frontieres de plaques 67 

6.3 Mecanisme au foyer des tremblements de terre 68 

6.3.1 Description des failles 68 

6.3.2 Determination des mecanismes au foyer 68 

6.3.3 Application 70 

6.3.4 Sismicite intraplaque 75 

6.4 Notes supplementaires: Tenseur des contraintes. 77 

6.4.1 Contraintes principales 78 

6.4.2 Effet de la friction 80 

6.4.3 Changement de systeme d’axes de coor- 

donnees 81 

7 Paleomagnetisme 82 

7.1 Introduction 82 

7.2 Champ magnetique terrestre 82 

7.3 Proprietes magnetiques des roches 85 

7.4 Determination des paleopoles 87 

7.5 Applications du paleomagnetisme 89 

7.6 Appendice: Determination du paleopole 93 

8 Applications de la geodesie spatiale a la tectonique 

des plaques 95 

8.1 Telemetrie satellite. SLR (Satellite Laser Ranging) 96 

8.2 GPS-Global Positioning System 97 

8.3 Interferometrie a large base: VLBI (Very Long 

Baseline Interferometry) 98 

8.4 Altimetrie satellite 100 

8.5 Interferometrie radar. (SAR-Synthetic Aperture 

Radar) 100 

8.6 Application en planetologie 102 

8.7 Applications environnementales 103 

8.8 Methodes 103 

8.8.1 Effet Doppler 103 

8.8.2 Interferometrie 103 

8.8.3 Mesure de la distance par GPS 104 

8.8.4 Satellite Laser Ranging 105 

8.8.5 La mission des satellites GRACE 106 

9 Evolution de la lithosphere oceanique. 

Mecanismes de la tectonique des plaques. 112 

9.1 Modeles devolution thermique de la lithosphere 

oceanique 112 


exosup.com 


page facebook 


VI 


Contents 


9.1.1 Variations du flux de chaleur avec l’age 112 

9.1.2 Variation de la bathymetrie des oceans 

avec l’age 113 

9.1.3 Effet des points chauds 115 

9.1.4 Applications 116 

9.2 Forces agissant sur les plaques 116 

9.2.1 Poussee de la dorsale (ridge push) 116 

9.2.2 Traction de la zone de subduction (slab pull) 116 

9.2.3 Forces sur les failles transformantes 118 

9.2.4 Friction a la base de la plaque (viscous drag) 118 

9.2.5 Contraintes intraplaques 119 

9.2.6 Stability des chaines de montagne 120 

9.3 Mecanismes de la Tectonique des plaques 121 

9.3.1 Organisation de la convection dans le 

manteau 122 

9.4 Flux de chaleur 126 

9.5 Modele de refroidissement de la lithosphere 127 

9.6 Subsidence thermique 128 

9.7 Convection thermique 128 

9.7.1 Viscosite d’un fluide 128 

9.7.2 Nombre de Rayleigh critique 129 

10 Rifts intracontinentaux. 

Marges continentales passives. Bassins sedimentaires. 136 

10.1 Rifts intracontinentaux 136 

10.1.1 Comparaison des rifts Est Africain et du 

Rio Grande 136 

10.1.2 Mecanisme de formation d’un graben 139 

10.1.3 Modeles devolution 140 

10.1.4 Rift oceanique 141 

10.2 Rifts avortes. Aulacogenes. 142 

10.2.1 Rifts avortes en Amerique du Norcl 142 

10.2.2 Aulacogenes 142 

10.3 Marges continentales passives. Formation des 

bassins sedimentaires. 142 

10.3.1 Subsidence des marges continentales 142 

10.3.2 Bassins sedimentaires intracontinentaux 143 

10.4 Chronologie des rifts de 1’Afrique de l’Est et du 

Rio-Grande 145 

10.4.1 Rio-Grande 145 


exosup.com 


page facebook 


Contents vii 

10.4.2 East Africa 145 

11 Zones de convergence 1: Subduction 146 

11.1 Subduction oceanique 146 

11.1.1 Morphologie 147 

11.1.2 Sismicite 147 

11.1.3 Gravite 147 

11.1.4 Structure sismique profonde 148 

11.1.5 Volcanisme 148 

11.1.6 Bassins marginaux 148 

11.1.7 Prismes d’accretion 148 

11.1.8 Metamorphisme des zones de convergence 149 

11.2 Marges continentales actives 149 

11.2.1 Cordillere des Andes 149 

11.2.2 Cordillere nord americaine 150 

11.2.3 Bassins d’avant pays (Foreland basins) 150 

11.3 Megathrust earthquakes in the Cascadia subduc- 
tion zone 151 

11.4 Notes supplementaires: Flexure de la lithosphere. 152 

12 Collision continentale 155 

12.1 Zones de convergence 2: Collision continentale 155 

12.2 Himalaya. Tibetan plateau 160 

13 Evolution tectonique des planetes internes 162 

13.1 Proprietes des planetes internes 162 

13.2 Atmosphere et Climat 162 

13.3 Age apparent de la surface 164 

13.4 Gravite et Topographie 165 

13.5 Activite et evolution tectonique 169 

13.6 Le systeme Terre-Lune. Origine de la Lune 169 

13.7 Conclusion 171 

13.8 Lois de conservation 172 

14 Remarques finales 173 

Appendix A Anomalies de gravite. 175 

Index 180 

References 181 


exosup.com 


page facebook 


exosup.com 


page facebook 


Introduction 


Objectifs du cours 

Ce cours resulte de la fusion de deux cours: Tectonique des plaques et 
Physique du Globe. Une introduction serieuse a ces deux sujets exigerait 
au minimum un cours complet pour chacun. Dans l’espace de ce cours, 
nous allons poursuivre trois objectifs principaux: 

• Examiner les methodes qui permettent de mettre en evidence et determiner 
le mouvement relatif des plaques tectoniques. 

• Etudier comment le mouvement des plaques est lie a la structure in- 
terne et a la dynamique du manteau terrestre. 

• Etudier comment la tectonique des plaques explique et integre l’ensemble 
des phenomenes geologiques dans le cadre du cycle de Wilson, cycle 
d’ouverture et fermeture de bassins oceaniques. 

Un objectif accessoire du cours est de montrer Fertilisation de modeles 
physiques simples. Le but n’est pas de developper des modeles mathematiques 
complexes, mais de montrer comment comprendre et resoudre des problemes 
geologiques a l’aide de simples raisonnemments physiques. Ce sont les 
mesures geophysiques qui ont permis de l’avenement de la tectonique 
des plaques, il y a maintenant 40 ans. Au cours des 20 dernieres annees, 
le developpement des methodes de tomographie sismique et de geodesie 
spatiale a permis d’explorer aux limites de la tectonique des plaques. 


References 

II existe de nombreux et excellcnts ouvrages de reference en anglais. 
Malheureusement, il n’existe aucun ouvrage en frangais qui puisse servir 
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de reference pour ce cours. Parmi les textes en anglais, les plus utiles 
seront Cox and Hart (1991), Kearey et al. (2009) et Davies (1999). 
Bott (1982) est superbement ecrit et agreable a lire, malheureusement 
pas assez quantitatif et maintenant un peu depasse. Condie (2004) est 
une nouvelle edition de Condie (1997) et reste une bonne introduction 
pour geologue, rnalgre le changement opportuniste du titre. Les livres 
de Fowler (2005), Schubert et al. (2001), Stacey and Davis (2008) sont 
de niveau plus avances. . 
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Liens WWW 


De nombreuses informations et des donnees sont directement accessibles 
par l’intermediaire du web. La liste des sites ci dessous n’est pas exhaus- 
tive. La plupart de ces sites offrent aussi un interface pour visualiser les 
cartes avec Google Earth. 

National Geophysical Data Center. NOAA. Boulder (CO), 
http: / /www. ngdc.noaa.gov 

Ce site offre differentes donnees globales (topographie, bathymetrie, gravite, 
climatologie, etc.) qui peuvent etre dechargees ou visualisees. 

USGS National Earthquake Information Center HOME PAGE 
http: / / earthquake.usgs.gov 

Ce site presente toutes les donnees d’activite sismique en temps reel, 
les mecanismes au foyer des tremblements de Terre les plus recents sont 
aussi affiches. 

USGS National Earthquake Information Center World Data Center 
A for Seismology 

http : / /earthquake . usgs . gov/regional / world / seismicity / 

USGS Volcano hazards program 
http: / /volcanoes. usgs. gov/ 

Presente l’etat de l’activite volcanique et des liens a d’autres sites. 

NASA National Aeronautics and Space Administration 
http://www.nasa.gov/ 

Le site presente des images de la Terre vue de l’espace ainsi que les 
programmes d’exploration planetaire de NASA. 

Modele tomographique global de vitesse des ondes S dans manteau 
superieur. http://ciei.colorado.edu/ nshapiro/MODEL/ 

Des liens avec certains de ces sites ainsi que de nombreux autres sites 
d’interet se trouvent sur le site WEB de la commission geologique du 
Canada. 

http: / / ess.nrcan.gc.ca/index_e.php 

GPlates is desktop software for the interactive visualisation of plate- 
tectonics . http: / /www .gplates.org/ 
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Historique 


La theorie de la tectonique des plaques a permis de placer tous les 
phenomenes geologiques dans un cadre explicatif tres simple qui peut 
se resumer en quelques phrases. 

L’enveloppe externe de la Terre est formee de plaques rigides. Le mou- 
vement relatif de ces plaques produit des deformations concentrees a 
leurs frontieres. La separation de deux plaques conduit a la formation 
du plancher oceanique qui est detruit dans les zones de convergences. 
Tous les phenomenes geologiques et revolution des continents sont ainsi 
lies a des cycles d’ouverture et fermeture de bassins oceaniques (cycles 
de Wilson). Le mouvement des plaques est entretenu par et participe au 
refroidissement du manteau terrestre par convection. 


1.0.1 Revolution scientifique et paradigme. 

De nombreuses etudes ont tente de formaliser la demarche de la recherche 
scientifique. Au vingtieme siecle, il faut surtout citer les travaux de Pop- 
per, Lakatos, Kuhn, et Feyerabend. 

La tectonique des plaques a ete presentee comme un exemple de 
revolution scientifique selon Kuhn (Kuhn, 1962). D’apres lui, la sci- 
ence progresse surtout au cours de grandes crises, les revolutions sci- 
entifiques. Entre ces revolutions, la science normale precede a l’interieur 
de paradigmes. Le paradigme est Lensemble des regies acceptees par la 
communaute. En temps normal, ces regies sont celles que l’on trouve 
dans les livres de classe. La recherche scientifique fonctionne dans le 
cadre de ces regies et ne les remet pas en question. De temps en temps, 
des anomalies sont observees. La premiere reponse consiste a amenager 
l’ancienne theorie. Si les anomalies se multiplient et que la theorie sem- 
ble irreparable, il se produit une crise qui ne peut se resoudre que par 
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un changement de paradigme. La nouvelle theorie doit d’une part con- 
tend l’ancienne comme cas particular, et d’autre part expliquer les 
phenomenes qui etaient inexplicables par l’ancienne. 

L’explication de Kuhn ne fait pas l’unanimite parmi les historiens et 
philosophes de la science. Sir Karl Popper avait essaye de montrer que 
la recherche scientifique obeit a une veritable logique (Popper, 1934). II 
a montre qu’une theorie ne peut jamais etre prouvee rnais simplcment 
etre rendue plus probable. Pour la rendre plus probable, il faut la tester 
experimentalement. La demarche scientifique consiste done a faire des 
experiences critiques qui permettent de tester la theorie et de la rendre 
plus probable. 

Pour Imre Lakatos, la demarche scientifique est definie par la mise en 
competition de programmes scientifiques. Un programme de recherche 
consiste en un noyau dur (les paradigmes) et une ceinture protectrice 
d’hypotheses supplementaires (Lakatos, 1978). La demarche consiste a 
explorer ces hypotheses peripheriques. Le choix de la communaute sci- 
entifique de poursuivre un programme est un choix rationnel qui est 
d’abord base sur sa rentabilite. Un programme deperit et est abandonne 
lorsqu’il cesse de produire des resultats nouveaux. 

Finalement, Feyerabend a une vision beaucoup moins rationaliste. Ainsi, 
il afhrme que la science est une entreprise essentiellement anarchiste qui 
ne precede pas d’une demarche aussi organisee et rigoureuse que celle 
proposee par Popper. Il laisse beaucoup plus de place au hasard dans la 
decouverte scientifique ( serendipity ) (Feyerabend, 1975). 


1.0.2 Les grands debats en Geologie 

La tectonique des plaques n’etait pas la premiere tentative d’une theorie 
generale de la geologie et la controverse sur la derive des continents 
n’etait pas le premier debat a agiter les sciences de la Terre. A la fin du 
XVIII siecle, un premier grand debat a oppose le neptunisme de Werner 
au plutonisme de Hutton. Sans entrer dans le detail, les neptunistes 
proposaient que toutes les roches avaient ete formees dans les oceans 
et que l’histoire de la Terre se resumait a l’emergence progressive des 
continents et au refroidissement des oceans. Ceci explique revolution de 
la nature des roches depuis les roches formees dans des oceans chauds 
(granite, basalte) aux roches formees dans des oceans refroidis (calcaire, 
argile) . Les roches les plus anciennes se trouvent ainsi au coeur des con- 
tinents. Hutton a reconnu la difference entre roches ignees (primaires) 
et les roches sedimentaires (secondaires) qui en sont derivees. Ses obser- 
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vations d’intrusions et d’inconformites lui ont fait suggerer que l’histoire 
de la Terre est une serie de cycles qui se repetent. No vestige of a begin- 
ning, No prospect for an end. 

Les differences entre plutonisme et neptunisme sont a la fois methodologiques 
et philosophiques. Le neptunisme etait base sur une analyse globale a 
l’echelle du continent europeen. La theorie est basee sur les tendances 
que l’on peut interpreter sur la carte geologique de l’Europe. Le pluton- 
isme, au contraire est basee sur des observations locales. La deuxieme 
difference concerne la nature du temps. Pour les neptunistes, le temps 
a une direction definie avec une evolution irreversible. Pour les pluton- 
istes, au contraire, les phenomenes geologiques changent pas au cours 
du temps et les cycles se repetent indefiniment. 

Au cours du XIX siecle un second debat a oppose le catastrophisme 
des paleontologues frangais a V uniformitarisme de Lyell. Pour les catas- 
trophistes, les extinctions ne peuvent s’expliquer que par des catastro- 
phes majeures qui expliquent toute l’histoire de la Terre. Pour les uni- 
formitaristes, la Terre a evolue suivant les mernes processus qui agissent 
et que Ton peut observer aujourd’hui. 

Les reponses aux questions soulevees par tous ces debats sont moins 
tranchees qu’il peut paraitre. La theorie de Wegener a fait renaitre 
l’opposition entre modeles cycliques et unidirectionnels. Recement, le 
catastrophisme est reapparu lorsqu’on a reconnu l’importance des im- 
pacts meteoritiques tant pour les extinctions (Fin du Cretacee) que pour 
la formation de certaines structures geologiques (Sudbury). Aujourd’hui, 
de nombreux debats sont encore alimentes par les contradiction entre 
modeles bases sur des donnees globales et ceux bases sur des donnees 
locales. 

1.0.3 La contribution de Wegener: la derive des 

continents 

Bien que la similarity des cotes de part et d’autre de l’Atlantique Sud ait 
provoque de nombreuses speculations auparavant, ce n’est qu’au debut 
du XX e siecle que les consciences geologiques d’une derive des con- 
tinents ont ete formulees de fagon precise. En 1910, Taylor a suggere 
que les chaines de montagnes sont causees par le deplacement des con- 
tinents. D’apres son schema, le deplacement de l’Eurasie vers le Sud et 
le deplacement des Ameriques vers l’ouest causent les chaines alpines et 
himalayennes et des cordilleres americaines. 

En 1915, Wegener a publie la premiere edition de son livre L’origine 
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des continents et des oceans (Wegener, 1929). Cette contribution est im- 
portante parce que Wegener a propose de faire une serie d’observations 
qui seraient expliquees par sa theorie. Prises individuellement, ces ob- 
servations ne suffisent pas a demontrer la validity? de l’hypothese de We- 
gener. Mais c’est l’ensemble de ces observations qui a donne du poids 
a la theorie de Wegener. Les observations les plus importantes sont les 
suivantes: 

• La similarity? des cotes. 

• Les arguments paleontologiques: l’apparition et la disparition d’especes 
simultanement de part et d’autre de l’Atlantique. De nombreux paleontologues 
suggeraient que des ” ponts” continentaux auraient existe et se seraient 
effondres pour donner naissance a l’ocean atlantique. 

• Wegener a fait grand cas d’arguments bases sur l’isostasie (voir sec- 
tion 2.2) pour demontrer qu’il n’etait pas possible que les oceans se fer- 
ment par effondrement des continents. Cet argument est faible parce 
qu’il est possible de transformer des continents en oceans par d’autres 
mecanismes qui preseveraient l’isostasie (changements de phase). De 
plus, ces arguments montrent que les continents ne peuvent se trans- 
former en oceans mais ne demontrent pas le deplacement des conti- 
nents. 

• Les arguments paleoclimatiques sont bases sur les anciennes glacia- 
tions en Inde, les anciens climats desertiques en Angleterre. 

• Les arguments stratigraphiques sont bases sur la similarity? de sequence 
de sediments de part et d’autre de l’Atlantique sud. 

• La similarity? des provinces geologiques de part et d’autre de l’Atlantique 
nord et sud a egalement servi d’argument a Wegener. 

• Wegener a aussi propose que de precises mesures astronomiques per- 
mettraient de demontrer le mouvement relatif des continents. Mal- 
heureusement, Wegener avait surestime le taux d’ouverture des oceans 
(20m/a) . II a ainsi publie des mesures qui auraient montre un deplacement 
de plusieurs centaines de m par an du Greenland par rapport a 1’Europe. 

L’approche de Wegener, qui favorise les observations globales et revolution 
a sens unique de la Terre depuis 1’eclatement de Pangee, a ete critiquee 
comme la renaissance de l’approche neptuniste. La plupart des propo- 
sitions de Wegener sont aujourd’hui verifiees. Malgre son exceptionnelle 
intuition, Wegener avait laisse de grossieres erreurs dans ses travaux. 

L’une de ces erreurs provient du manque de chronometrie absolue et done 
de la surestimation des taux de deplacement. Une autre erreur tient au 
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mecanisme propose par Wegener. Wegener proposait que la croute conti- 
nentale ” fiottait” sur la croute oceanique et ces ’’radeaux” continentaux 
se deplagaient sur le plancher oceanique. D’apres Wegener, la force qui 
permettait le deplacement des continents etait dirigee vers l’equateur et 
causee par une legere difference entre la resultante des forces de gravite 
et de la force d’Archimede. 


1.0.4 Reponse a l’hypothese de Wegener. 

Les geophysiciens, en particulier Sir Harold Jeffreys ont immediatement 
attaque les faiblesses du mecanisme propose par Wegener. Les mesures 
sismiques montraient que la croute oceanique est plus dure que la croute 
continentale (les modules elastiques sont plus eleves) et les continents ne 
pouvaient done pas simplement ecarter la croute oceanique. Par ailleurs, 
la resultante des forces suggerees par Wegener est ridiculement faible. 

En depit de ces objections, l’hypothese de Wegener a ete serieusement 
consideree pendant une dizaine d’annees. Ce n’est qu’apres le sympo- 
sium de l’A.A.P.G., en 1928, que la majorite des geologues a ignore la 
theorie de Wegener. Ceci est du en partie aux erreurs dans les mesures 
astronomiques et la mauvaise chronometrie de Wegener. 

Apres la mort de Wegener en 1929, une minorite de chercheurs a 
continue a defendre sa theorie. II faut surtout citer: 

• Le geologue sud-africain duToit a travaille sur la similarity de la strati- 
graphie de part et d’autre de l’Atlantique. Argand a etudie la tec- 
tonique de l’Asie et avait conclu que seule la convergence de l’lnde 
et de l’Asie pouvait produire les grandes nappes himalayennes. A 
l’epoque, ces arguments ont eu peu d’impact car la majorite des 
geologues pensaient que d’enormes deplacements horizontaux pou- 
vaient etre induits par des mouvements verticaux. 

• Le geophysicien Arthur Holmes qui a propose que le deplacement 
des continents et l’ouverture de bassins oceaniques etaient dus a des 
courants de convection dans le manteau terrestre. Jeffreys pensait que 
les chaines de montagne etaient causees par le refroidissement et la 
contraction thermique de la Terre. Holmes pensait que le manteau 
sous-continental etait chauffe par la radioactivite des roches continen- 
tales et que cela causait la rupture des continents et l’ouverture de 
bassins oceaniques. 

• Les travaux experimentaux de Griggs (1940) sur la rheologie des roches 
ont montre que la rheologie dependait de l’echelle de temps consideree. 
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Une roche pouvant etre elastique pour des phenomenes sismiques et 
visqueuse a l’echelle de quelques milliers d’annees. 

Le debat sur la derive des continents fait partie d’un debat plus general 
sur le mecanisme de refroidissement de la Terre. La question avait deja 
ete soulevee au 19 e siecle lorsque Lord Kelvin avait calcule l’age de 
la terre a partir du gradient thermique. Le calcul de Kelvin ne tenait 
pas compte des sources de chaleur interne (la radioactivite) et surtout 
supposait que la Terre se refroidissait uniquement par conduction. La 
decouverte de la radioactivite n’a pas resolu la question soulevee par 
Kelvin car le probleme principal est celui du mecanisme de refroidisse- 
ment. Jeffreys et la majorite des geopliysiciens ont continue a soutenir, 
coniine Kelvin, que la Terre se refroidissait uniquement par conduction. 
Par contre, c’est parce qu’il a compris que la conduction n’etait pas un 
mecanisme efHcace de refroidissement de la Terre que Holmes a propose 
que la manteau terrestre ne pouvait se refroidir que par convection. In- 
cidemment, la convection permet une evolution cyclique tant qu’il y a 
des sources de chaleur, tandis que la conduction donne une direction au 
temps. 


1.0.5 Exploration des fonds marins. Paleomagnetisme. 

Les techniques developpees pendant la guerre pour lutter contre les sous- 
marins ont servi a explorer les fonds marins apres 1945. Par ailleurs, 
pour des raisons militaires, de nombreux programmes d’exploration et de 
cartographie des fonds marins se sont developpes. Parmi les precurseurs 
de l’apres guerre, il faut signaler Vening-Meinesz qui fit des mesures de 
gravite dans les fosses oceaniques, Bullard qui fit les premieres mesures 
de flux de chaleur, et M. Ewing qui developpa la sismique marine. 

Le paleomagnetisme fut developpe par Blackett et Runcorn et leurs 
etudiants, en particulier Ted Irving au Canada. Leurs travaux furent 
longtemps ignores par la communaute scientifique. 

C’est la demonstration des inversions du champ magnetique et la rnise 
au point de la clrronometrie precise de ces inversions (Cox, Doell, Dal- 
rymple) qui a permis a Morley, au Canada, et a Vine et Matthews, en 
Angleterre, de proposer l’hypotliese d’ expansion des fonds oceaniques. 

Le concept de plaque est du en grande partie du a J. Tuzo Wilson 
(de l’Universite de Toronto) qui a reconnu la nature et le role des failles 
transformantes. Le concept a ete generalise et introduit en 1967 a la 
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conference de l’AGU par W. Jason Morgan, Dan P. McKenzie et Robert 
L. Parker. 


1.1 Remarques 

II faut noter que revolution des idees a ete tres rapide des que des 
mesures tres precises ont permis de verifier l’hypothese d’expansion des 
fonds oceaniques. Entre 1962 et 1967, la grande majorite des geophysiciens 
occidentaux a change d’opinion et a adopte le nouveau paradigme. II faut 
aussi noter que de nombreux chercheurs de tres grande reputation sont 
restes jusqu’a leur mort opposes a la tectonique des plaques (Jeffreys, 
Ewing, Belloussov). 

Remarquons que la tectonique des plaques, telle que formulee, il y 
a 30 ans, n’avait plus grand chose a voir avec la theorie de Wegener. 
Depuis la theorie a encore beaucoup evolue. Un premier point a ete 
la reconnaissance progressive que la tectonique des plaques avait opere 
pendant la plus grande partie de l’histoire de la Terre. C’est ainsi que la 
geologie des continents s’explique par la succession de cycles de Wilson 
d’ouverture et fermeture de bassins oceaniques. Par ailleurs, les progres 
technologiques, surtout ceux des mesures satellitaires permettent main- 
tenant de rnieux mesurer les deformations tectoniques. On a ainsi pu 
etablir que les deformations ne se limitent pas toujours aux frontieres de 
plaques. Dans le cas de collisions continentales, l’interieur des plaques 
peut-etre tres profondement deforme. 

L’avenement de la tectonique des plaques a aussi marque le declin, 
sinon la disparition, de la ’’geologie de terrain” dans la recherche en 
sciences de la Terre. Ce sont les donnees geophysiques qui ont etabli la 
theorie. Meme dans le cas des failles transformantes, ce sont les donnees 
sismiques qui ont etabli leur nature et leur role. Aujourd’hui, la geodesie 
spatiale et la seismologie sont les outils privilegies pour l’etude de la 
tectonique active. 

L’ evolution des sciences de la terre peut etre decrite de differentes 
fagons. L’interpretation de Kuhn doit etre corrigee car il n’y avait pas 
de paradigme avant l’emergence de la tectonique des plaques, sauf peut- 
etre en Russie ou l’academicien V.V. Belloussov avait decrete que toute 
la tectonique est due a des deplacements verticaux. La tres grande ma- 
jorite des geologues et geophysiciens sovietiques s’est done raliee a la 
tectonique des plaques avec 20 ans de retard. Mais, la science dans 
l’URSS stalinienne at post-stalinienne avait une logique qui echappe a 
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toute analyse logique. Le concept de programme scientifique de Lakatos 
s’applique bien pour expliquer l’attitude et les choix de la communaute 
scientifique. On peut considerer que la contribution de Wegener etait 
la formulation d’un programme de recherche. Ce programme etait peu 
’’rentable” en 1928. Pour Pauli, le pire defaut d’une theorie est de n’etre 
rnerne pas fausse, cad qu’elle est impossible a tester. C’ etait un peu le 
cas de la theorie de Wegener dans les annees 1920s. C’est parce que le 
programme d’exploration de la structure interne de la Terre produisait 
des resultats nouveaux que la majorite des geophysiciens l’ont choisi en 
1928. L’avenement de la tectonique des plaques marque autant le declin 
de ce programme que le succes du programme concurent. 
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2.1 Structure sismique de la Terre 

La structure interne de la Terre (Fig 2.2) a ete definie par la sismologie. 
En premiere approximation, la terre est une sphere symetrique et ses 
proprietes (vitesses sismiques, etat physique et densite) ne dependent 
que de la profondeur. L’etat physique (solide, liquide) est determine par 
la presence d’ondes de cisaillement (ondes S) qui ne se propagent pas 
dans les liquides. 

Du point de vue sismique, la structure interne de la Terre est divisee en 
trois grandes regions (croute, manteau, noyau) separees par des discon- 
tinuites de vitesse (Jeffreys, 1976; Bullen, 1984). II existe une structure 
plus fine a l’interieur de chacune de ces regions. 

• La croute dont l’epaisseur varie entre 5km sous les oceans et 40km 
sous les continents. L’epaisseur de la croute oceanique est a peu pres 
constante (5-6km). La croute continentale est souvent epaissie sous les 
chaines de montagne, jusqu’a 80km dans l’Himalaya. 

• Le manteau, qui est solide jusqu’a 2900km de profondeur. Deux zones 
de transition a 450knr et 650km separent le manteau superieur du 
manteau inferieur. A la base du manteau, la couche D” est probable- 
ment une zone oil le gradient de temperature est eleve. Les etudes de 
tomographie sismique montrent que cette couche est tres heterogene. 

• Le noyau externe est liquide et ne transmet pas d’ondes de cisaille- 
ment. II est compose a 90% de Fe mais contient d’autres elements. 

• Le noyau interne est partiellement solide. Son etat est solide a ete 
deduit de l’augmentation de la vitesse des ondes de compression. Ce 
n’est que tres recemment que l’on a identifie l’onde PKJKP (onde P 
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Figure 2.1 Profils moyens de vitesses sismiques et densite a l’interieur 
de la la terre. Les discontinuites marquent les grandes enveloppes. 
L’absence d’onde de cisaillement dans le noyau externe implique qu’il 
est liquide 



convertie en onde de cisaillement qui aurait traverse le noyau interne 
avant d’etre reconvertie en P)(Cao et al., 2005). 

La discontinuite de Mohorovic et la discontinuite manteau-noyau cor- 
respondent a des changements de composition. Les transitions entre le 
manteau superieur et inferieur et entre les noyaux internes et externes 
sont interpretees comme des changements de phase. 


2.2 Isostasie 

Les premieres mesures des deviations de la verticale a proximite de 
chaines de montagnes furent faites en Amerique du Sud par l’astronome 
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Figure 2.2 Structure interne de la Terre 



Pierre Bouguer. Bouguer constata que l’attraction gravimetrique de la 
chaine Andine etait beaucoup plus faible que ce qu’il avait calcule. Les 
conclusions de Bouguer furent confirmees par les mesures faites par Ever- 
est en Inde. II observa egalement que l’attraction gravitationnelle de la 
chaine himalayenne etait beaucoup plus faible que prevue. Pratt et Airy 
en deduirent leurs modeles d’isostasie. Leurs modeles impliquent que 
la masse de la chaine de montagne est compensee en profondeur par 
un deficit de masse. Cette hypothese est confirmee par la correlation 
negative entre les anomalies de Bouguer et la topographie (voir Figures 
2.4 et 2.4 et appendice A). Selon Pratt, c’est parce que la croute est 
moins dense que les chaines de montagnes sont elevees. Selon Airy, la 
croute de densite moyenne constante et d’epaisseur variable, repose sur 
un manteau plus dense (comme un iceberg). Dans les deux modeles, la 
masse totale d’une colonne de roche entre la surface de la Terre et le 
niveau de compensation (auquel toutes les differences disparaissent) est 
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Figure 2.3 Variations de g et de la pression depuis le centre jusqu’a 
la surface de la Terre. On peut noter que g est peu pres constant a 
travers le manteau. 



constante. II en resulte que la pression est constante au niveau de com- 
pensation. Cette situation correspond a la condition d’equilibre dans un 
fluide (hydrostatique) . La pression a la profondeur z comp est: 

T^comp 

P(zcomp) = / p{z)gdz (2.1) 

J Zq 

oil z est la profondeur, zq est la surface de la Terre, g est l’acceleration 
de la gravite, et p est la densite. 

Le fait que l’isostasie soit observee dans les bassins sedimentaires ou 
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Figure 2.4 Topo^aphie du continent nord Americain. 
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dans les orogenes quasi independament du niveau d’erosion implique 
qu’il y a en permanence reajustement en reponse au chargement de la 
surface de la Terre. Ceci a ete confirme par l’etude du rebond post- 
glaciaire. 

Les mecanismes de Pratt et d’Airy sont des mecanismes locaux. En 
fait, l’isostasie est un phenomene qui n’est valide que regionalement 
(a 1’echelle de plusieurs centaines de km). Ainsi, c’est l’ensemble des 
Rocheuses qui sont compensees, rnais chaque vallee ou montagne n’est 
pas compensee individuellement. Le mecanisme de compensation regionale 
fait appel a la flexure d’une lithosphere elastique cjui supporte les charges 
de courtes longueurs d’onde. 
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Figure 2.5 Anomalie de Bouguer du continent nord Americain. Notez 
que l’anomalie de Bouguer est d’autant plus negative que Paltitude 
est elevee, suggerant une densite plus faible et/ou un croute plus 
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2.3 Structure rheologique 

II est egalement possible de definir dans le manteau superieur une lithosphere 
et une asthenosphere. Ces definitions ont ete proposees d’abord pour ex- 
pliquer le concept d’isostasie. La presence d’une region aisement deformable 
en profondeur etait necessaire pour permettre les ajustements isosta- 
tiques. La lithosphere est done l’enveloppe externe de la Terre qui resiste 
aux deformations; 1’asthenosphere est une couche plus malleable. La 
lithosphere ne coincide pas necessairement avec les plaques bien que les 
concepts soient proches. De meme, les definitions basees sur les pro- 
prietes mecaniques ou thermiques ne coincident pas exactement. La 
definition de 1’ asthenosphere est restee assez floue. Les sismologues avait 
identifie la base de la lithosphere avec un region du manteau a vitesse 
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sismique reduite (LVZ = low velocity zone). Cette definition est main- 
tenant demodee. Une definition rheologique de la lithosphere est basee 
sur les mecanismes de deformation des roches entre une region du man- 
teau qui ne se deforme pas a l’echelle geologique et et une zone qui 
peut se deformer (rheologie visqueuse). On utilise souvent un isotherme 
(1300°C) pour definir la base de la lithosphere. On a aussi introduit une 
epaisseur elastique effective de la lithosphere qui n’a rien a voir avec les 
definitions precedentes. 


Un solide elastique ne peut se deformer de fagon permanente. La deformation 
est proportionnelle a la contrainte et le solide retrouve son etat initial 
en l’absence de contrainte. Par contre un fluide se deforme sous la con- 
trainte et il ne retrouvera pas son etat initial quand la contrainte n’est 
plus appliquee. 

2.3.2 Deformation continue des roches 

Dans un fluide, la deformation e augmente tant que la contrainte a est 
appliquee. Le cas le plus simple est le cas du fluide newtonien pour lequel: 


Le taux de deformation est done proportionnel a la contrainte. La con- 
stante de proportionnalite v est la viscosite du fluide. En fait, les experiences 
montrent que si une contrainte est appliquee pendant longtemps a un 
solide, il ne reviendra pas exactement a son etat initial apres disparition 
des contraintes. Ceci est aussi le cas pour les roches. Les experiences 
entreprises par Griggs sur la deformation de roches soumises a des con- 
traintes de longues durees ont ete faites systematiquement sur la plu- 
part des roches. Ces experiences ont montre que les roches et mineraux 
peuvent subir une deformation permanente et continue. Le taux de 
deformation e depend de la composition, de la temperature, et de la 
contrainte a ; en premiere approximation, l’effet de la pression est sou- 
vent neglige. Les lois de deformation sont done du type : 


2.3.1 Rheologie 


a = ve 


( 2 . 2 ) 


e = Aa n exp (-H/RT) (2.3) 



(2.4) 
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ou A, n, et l’enthalpie d’activation H sont des constantes dependant 
de la composition. R = 8.3144J mol -1 K -1 est la constante des gaz, 

T est la temperature thermodynamique. Dans la croute n ~ 3 tan- 
dis que dans le manteau n = 1 (fluide newtonien). Cette loi est done 
non lineaire dans la croute ou la contrainte n’est pas proportionnelle 
au taux de deformation. A faible contrainte, le taux de deformation est 
quasi nul, rnais ce taux augmente beaucoup plus rapidement que la con- 
trainte. Ainsi, pour n = 3, la deformation est multipliee par 8 alors 
que la contrainte double. Reciproquement, la contrainte necessaire pour 
maintenir un taux de deformation n’augmente pas aussi rapidement que 
celui-ci. Pour le manteau, les donnees experimentales et d’autres donnees 
geophysiques sont compatibles avec une loi lineaire. Dans ce cas 

a = vi (2.5) 

oil v est la viscosite. Les contraintes se mesurent en Pa, les taux de 
deformation en s^ 1 , la viscosite en Pa s. 

La figure 2.6 represente en fonction de la temperature la contrainte 
necessaire pour maintenir un certain taux de deformation. Cette con- 
trainte decroit tres rapidement avec la temperature. Ainsi, e’est la temperature 
qui determine si les roches se deferment continuement (regime ductile) 
ou par fracture (regime cassant). Lorsque la contrainte pour maintenir 
un taux de deformation continue est inferieure a la contrainte de cisaille- 
rnent sur un plan de fracture, le regime passe de cassant a ductile. La 
profondeur de cette transition depend de la temperature. 


2.3.3 Loi de Byerlee 

La seule deformation permanente pour un solide est la fracture (regime 
cassant, brittle). Un materiel se fracture lorsque la contrainte de ci- 
saillement sur un plan de fracture excede la friction le long du plan. 
Les experiences ont montre que cette lirnite ne depend que de la pres- 
sion (contrainte normale) et est independante de la composition des 
roches. La loi de Byerlee donne la valeur de la contrainte de cisaillement 
necessaire a s pour qu’il y ait mouvement sur une fracture: 

(j s = a(a n - Pf) (2.6) 

ou a n est la contrainte normale, Pf est la pression de fluide et 

a ~ 0.8 a n < 200 MPa 
a « 0.6 <j n > 200 MPa 
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Figure 2.6 Contrainte necessaire pour maintenir un taux de 
deformation donne en fonction de la temperature (Celsius) pour trois 
roches differentes: granite, granulite, et dunite. Ces roches correspon- 
dent a la croute superieure, inferieure, et au manteau. La temperature 
augmente et la contrainte diminue avec la profondeur. 


stress for 1 e— 1 6 1/s strain rate 



Le gradient de pression lithostatique est a peu pres 30MPakm -1 . La 
pression de fluide est souvent egale a la pression hydrostatique (lOMPa km -1 ). 
La contrainte de cisaillement necessaire pour causer le deplacement le 
long d’un plan de fracture devient considerable dans la croute inferieure. 

La plupart des tremblements de terre se produisent a une profondeur de 
moins de 20km. A cette profondeur, la contrainte de rupture est deja 
tres superieure aux contraintes tectoniques moyennes (20 MPa). II faut 
done que la croute soit localement affaiblie ou que la contrainte puisse 
etre concentree localement. 
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Figure 2.7 Contrainte necesaire pour deformer les roches de la 
lithosphere en fonction de la profondeur. Si la contrainte ductile 
est plus grande que cello pour fracturer une roche, le regime de 
deformation est cassant. A faible profondeur, les roches se deferment 
de maniere cassante; lorsque la profondeur augmente, on entre dans 
le regime ductile. Lorsque la composition change a mi-croute at dans 
le manteau, on retrouve le regime cassant. Notez que dans le regime 
cassant, la contrainte necessaire est plus importante en compression 
qu’en extension. 


strength profile 



stress (MPa) 
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3.1 Fonds oceaniques 

La geologie et la geophysique marines ont mene a la theorie de la tec- 
tonique des plaques, alors que la geologie et la geophysique dans les con- 
tinents avaient atteint une impasse. Ainsi que l’a fait remarquer J. Tuzo 
Wilson, ce n’est pas en regardant le pont d’un bateau mais la mer que 
Ton observe le mouvement du bateau. La theorie de l’expansion des fonds 
oceaniques ( sea floor spreading) fut une etape importante sur la voie de 
la tectonique des plaques. La geophysique marine s’est developpee apres 
la guerre 1939-1945. Le programme d’exploration systematique des fonds 
marins a mis en evidence des differences fondamentales entre oceans et 
continents. 


3.1.1 Morphologie des fonds oceaniques 

La bathymetrie des fonds marins est maintenant bien connue. Les elements 
essentiels sont les dorsales mid-oceaniques avec une vallee mediane (rift 
central), et leurs flancs qui descendent progressivement vers les bassins 
oceaniques. Les dorsales forment une chaine quasi continue sur 35,000km. 
Contrairement, aux chaines orogeniques dans les continents, dies sont 
formees exclusivement de roclies ignees et ne contiennent pas de roches 
sedimentaires ou metamorphiques. Les fosses profonds entourent le Paci- 
fique. La plupart des fosses profonds sont paralleles a des arcs vol- 
caniques insulaires (iles en arc). La concavite de ces arcs est toujours 
dirigees vers l’ocean. La profondeur des bassins oceaniques est « 4.5km 
alors que celles des dorsales est «2km; la dorsale emerge parfois, en Is- 
lande par exemple. Les fosses ont une profondeur de ss 10km. Le plateau 
continental peu profond fait partie du continent. C’est le talus qui definit 
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la transition entre les continents et les bassins oceaniques. A l’exception 
des dorsales, la bathymetrie des fonds oceaniques est reguliere. II existe 
quelques anomalies bathymetriques isolees sur les fonds oceaniques: les 
’’seamounts” et les ’’guyots” 1 . 

Cette structure globale a ete progressivement raffinee, surtout grace 
aux donnees satellitaires. Ces donnees revelent toutes les irregularites de 
la bathymetrie des fonds oceaniques. Ainsi, le nornbre de volcans sous- 
rnarins en particulier dans le Pacifique sud est beaucoup plus eleve qu’on 
le pensait. Par ailleurs, ces donnees revelent des bombements a grande 
echellc du plancher oceanique ( swells et superswells ) . 

3.1.2 Structure sismique de la croute oceanique 

Entre 1955 et 1965, un programme d’etude de la croute de l’ocean Atlan- 
tique par sismique refraction a ete developpe par M. Ewing et l’equipe 
du Lamont-Doherty Observatory. Le principe de la methode est rap- 
pele en Appendice. Les donnees de sismique refraction ont rapidement 
rnontre que la croute oceanique etait tres uniforme contrairement a la 
croute continentale. Son epaisseur (5-6km) ne varie pratiquement pas 
de la dorsale au talus. Les etudes sismiques ont permis de distinguer 3 
couches principales dans la croute oceanique. II faut remarquer que seule 
la couche superficielle, qui n’est pas presente a la dorsale, s’epaissit en 
s’eloignant de la dorsale. Cette epaisseur reste faible (< 1,000m). La fig- 
ure 3.1 rnontre une coupe composite de l’ocean Atlantique de la dorsale 
au plateau continental americain. 

• couche 1: vp ~ 1.5 — 2.5km s -1 

• couche 2: vp ~ 3.5 — 6.2km s -1 

• couche 3: vp ~ 6.5 — 7.5km s -1 

• manteau: vp > 8.1km s^ 1 


1 Le terme guyot a ete introduit par Hess. II s’agit de volcans sous-marin dont la 
partie superieure est plate parce qu’erodee par les vagues. 


exosup.com 


page facebook 


3.1 Fonds oceaniques 


25 


Figure 3.1 Structure sismique de la croute oceanique 


E 



Alors que les vitesses sismiques et l’epaisseur totale de la croute con- 
tinentale sont tres variables d’une region a l’autre, la croute oceanique 
est tres uniforme. La largeur de la zone de transition ocean-continent 
est tres variable (25 a >500km). 


3.1.3 Gravite 

Les anomalies d’air libre 2 de grandes longueurs d’onde (>200km) sont 
faibles sur les dorsales et les bassins oceaniques. Ceci implique qu’il y a 
compensation isostatique a cette echelle. Par contre, des anomalies de 
courtes longueurs d’onde sont bien marquees et refletent directement la 
bathymetrie des fonds oceaniques (i.e. il n’y a done pas compensation iso- 
statique pour les courtes et moyennes longueurs d’onde). Les anomalies 
positives correspondent a la masse supplementaire au dessus du planclier 
oceanique. Les plus fortes anomalies d’air libre au monde proviennent 
des fosses profonds qui montrent un couple d’anomalies, 1’ime negative 
vers le fosse, l’autre positive vers l’arc insulaire. 

Pour les courtes longueurs d’onde, la surface du geoide suit tres bien 
la bathymetrie. Les donnees satellitaires permettent ainsi d’interpoler la 
bathymetrie entre les profils recueillis par les vaisseaux oceanographiques. 


3.1.4 Flux de chaleur 

Bullard a fait les premieres mesures de flux de chaleur en mer en 1948. 

Les resultats de ces mesures ont ete tres surprenants. La croute continen- 

2 En mer, on utilise plus souvent les anomalies d’air libre que de Bouguer en partie 
parce que la bathymetrie n’est pas toujours connue avec assez precision pour les 
corrections. La correction de Bouguer se calcule en remplagant l’eau de mer par 
des roches de densite 2800 Mg m -3 . 
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tale est epaisse et radioactive. On attendait done un flux de chaleur beau- 
coup rnoins eleve dans les oceans que dans les continents. Les premieres 
nresures indiquaient un flux de chaleur moyen egal pour les oceans et 
continents. L’ensemble des nresures disponibles maintenant montre un 
flux oceanique superieur au flux continental. Le flux moyen global est 
estime a 87mW m~ 2 ; le flux continental moyen est de 65mW m~ 2 et 
le flux oceanique de lOlmW m -2 . Contrairement, aux autres champs 
geophysiques, le flux de chaleur varie beaucoup et de maniere systematique 
depuis les dorsales (> 250nrW m -2 ) jusqu’au bassins oceaniques (~ 
50mW m -2 ). 


3.1.5 Geologie marine 

On peut attribuer l’idee d’expansion des fonds oceaniques a Hess [1962]. 
Hess avait propose, uniquement sur la base des quelques donnees geologiques 
et de bathymetrie, que la croute oceanique etait plus jeune que la croute 
continentale et se formait aux dorsales. Hess suggerait notamment que 
les ’’guyots” etaient d’anciens volcans qui s’etaient formes aux flancs 
des dorsales et avaient ete erodes par Faction des vagues. Ils avaient 
ete transporters avec le plancher oceanique qui s’ecarte de la dorsale. 

Les anomalies magnetiques marines ont permis de tester de maniere 
systematique l’hypothese de Hess. 


3.2 Les inversions du champ magnetique terrestre. 

Deux decouvertes ont permis de verifier la theorie de l’expansion des 
fonds oceaniques: (1) la demonstration des inversions du champ magnetique 
terrestre, (2) la decouverte des bandes paralleles d’anomalies magnetiques 
marines et leur symetrie apparente par rapport aux dorsales. 


3.2.1 Anomalie magnetique 

En premiere approximation, le champ magnetique terrestre est celui d’un 
dipole, dont l’axe est legerement incline par rapport a Faxe de rotation 
de la Terre. Le champ magnetique total observe a la surface de la terre 
comprend une composante interne (qui provient du noyau), une com- 
posante crustale et une composante externe qui varie rapidement (Is — > 
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24h) dans le temps. * * 3 Les mesures du champ magnetique doivent etre cor- 
rigees pour eliminer les variations externes et le champ du noyau. Ces 
corrections permettent de definir les anomalies magnetiques, dont les 
sources sont les roches magnetisees de la croute et du manteau superieur. 
II faut rappeler que l’anomalie magnetique, qui depend de l’inclinaison 
du champ magnetique, n’est en general pas symetrique par rapport a la 
source. Lorsque le champ magnetique est vertical (au pole), l’anomalie 
magnetique est symetrique. Aux hautes latitudes (> 50 °), l’anomalie 
magnetique est presque symetrique. Aux basses latitudes, les anomalies 
sont asymetriques et beaucoup plus difficiles a interpreter. Toutefois, 
il est possible de traiter les anomalies magnetiques et de calculer ce 
qu’elles seraient si l’inclinaison etait 90 ° (reduction au pole). Ces anoma- 
lies reduites au pole sont symetriques. 


3.2.2 Proprietes magnetiques des roches et mineraux 

Certains mineraux ont des proprietes magnetiques. Ces proprietes sont 
decrites en appendice. La magnetisation M peut etre soit induite 4 , soit 
remanente. Cette magnetisation remanente est celle qui reste lorsqu’il 
n’y a plus de champ magnetique ambiant. En laboratoire, le phenomene 
d’hysterese magnetique montre l’acquisition d’une magnetisation remanente 
(en general en presence d’un fort champ ambiant). II y a plusieurs 
mecanismes d’acquisition de magnetisation remanente dans les mineraux. 

Pour la geodynamique, le mecanisme le plus important est la thermorcmancnce . 

La temperature de Curie (Tc) est la temperature au dela de laquelle 

toutes les proprietes magnetiques dispar aissent. La temperature de Curie 
est une propriete physique qui depend de la composition des mineraux. 

(Elle varie entre 450 °C et 700 °C). En presence d’un champ magnetique 
ambiant, un mineral qui se refroidit sous la temperature de Curie peut 
acquerir une magnetisation remanente. La magnetisation ne changera 
plus apres que la temperature descend sous la temperature de blocage, 

(« 50 °sous la temperature de Curie). La magnetisation remanente 
est done determinee par le champ magnetique ambiant pendant que 
la temperature est comprise entre la temperature de blocage et celle de 
Curie. II s’agit la de la composante primaire de la magnetisation. II est 
possible qu’une composante secondaire soit imposee par la longue ex- 

3 Le champ interne varie aussi dans le temps (a une echelle annuelle) et les 
corrections tiennent compte de ces variations. 

4 Dans ce cas, la magnetisation est proportionnelle au champ ambiant H: 

M = \H ou la susceptibility magnetique \ est une propriete physique des roches 
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position de la roche a un champ magnetique different. Ce phenomene 
appele viscoremanence est important en paleomagnetisme. 


3.2.3 Magnetisme inverse dans les roches 

En 1904, Brunhes fut le premier a remarquer que la magnetisation de 
certaines roches etait opposee au champ magnetique terrestre actuel. 
Cette observation n’implique pas que le champ magnetique s’est inverse. 
II se pourrait que ce soit la magnetisation de la roche qui s’inverse. Le 
debat a dure d’autant plus longtemps que les physiciens ont propose 
des mecanismes qui permettraient une inversion de la magnetisation. De 
plus, au moins un exemple fut rapidement trouve oil l’on pu etablir que 
c’est la magnetisation qui s’etait inversee. 


3.2.4 Chronometrie des inversions 

Pour demontrer que c’est bien le champ magnetique qui s’est inverse, 
il faut demontrer la synchroneite des polarites de la magnetisation de 
toutes les roches a la surface de la terre. Cette demonstration dut atten- 
dre la mise au point d’un outil geochronometrique suffisamment precis 
pour dater des roches volcaniques jeunes. Ce n’est qu’apres le developpement 
de la methode K-Ar que Ton put enfin demontrer la synchroneite des 
magnetisations et que Ton put etablir une chronometrie precise des inver- 
sions du champ magnetique terrestre. Cette echelle de temps est devenu 
un outil tres important (magnetostratigraphie) . 

Cette chronometrie fut etablie entre 1960 et 1970 par les etudes de 
l’age et de la polarite de la magnetisation dans des coulees recentes de 
lave. Les travaux les plus importants sont ceux de Allan Cox et ses 
collaborateurs bases au U.S.G.S. et ceux du groupe de McDougall en 
Australie. La decouverte de l’evement de Jaramillo, qui a permis de 
replacer correctement le debut de la periode actuelle de polarite normale, 
a ete la cle de la mise en place de cette chronometrie (Figure 3.2). 

La chronologie des polarites du champ magnetique distingue entre les 
epoques (de longues periodes, ~lMa) marquees par une polarite dorni- 
nante au cours desquelles interviennent de courts (<50,000a) evenements 
de polarite inverse. Ainsi, l’epoque normale presente est celle de Brun- 
hes qui a ete precedee par l’epoque de Matuyama. Un evenement in- 
verse (Laschamp) semble s’etre produit il y a 40,000a. L’evement de 
Jaramillo est un evenement de polarite normale durant l’epoque inverse 
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de Matuyama. II n’y a aucune regularity ou periodicite dans les inver- 
sions du champ magnetique. II sernble y avoir par le passe de tres longues 
epoques sans inversion (120— >85Ma). Depuis 75Ma, la frequence des in- 
versions sernble plus elevee qu’auparavant. 
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Figure 3.2 Chronometrie recente des polarites du champ magnetique. 


Time Epoch Polorily Polorily Polarity 

(Ma) chrons subchrons 




Brunhes 


c 



V 

o 

- 0 72 

1 




to 



CL 

Motuyomo 

2 





-2 47 



Gouss 


to 

c 



o 

o 

- 341 


CL 


A 


Gilbert 




5 5 

« 

c 

to 

U 

o 

- 5 44 


5 




Joromillo 

G'lSO 

Oiduvai 

Reunion 


Koena 

Mommoth 

Cochiti 

Nunivok 

SidutjoM 

Thvero 


3.2.5 Mecanisme des inversions 

II est maintenant bien etabli que le champ magnetique terrestre est 
produit par des courants de convection dans le noyau liquide [ Jacobs , 
1987]. Le mecanisme d’une dynamo auto-excitee avait ete propose par 
Bullard et par Elsasser. II suppose la pre-existence d’un faible champ 
magnetique. En presence de ce champ, les courants electriques dans le 
noyau amplifient le champ qui interagit avec le courant. Ce mecanisme 
avec retroaction est non-lineaire et semble ’’chaotique” 5 . Le champ 

5 Par chaotique, on entend un phenomene dont revolution est determinee par un 
systeme d’equations, mais dont le comportement est tellement sensible aux 
conditions initiales qu’il est pratiquement impossible de faire des previsions a 
long terme. Les tremblements de terre, le climat sont d’autres exemples de 
systemes chaotiques. 
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3.3 Anomalies magnetiques marines: Hypothese de Vine-Matthews-MorlSf 

magnetique pourrait done s’inverser spontanement. Des experiences de 
laboratoire ont montre que deux dynamos couplees pouvaient produire 
un champ qui s’inverse spontanement. Glatzmaier et Roberts [1995] ont 
recemment reussi a modeliser numeriquement une geodynamo qui s’inverse. 

A l’echelle geologique, l’inversion est un phenomene tres rapide (~l,000a). 
Le comportement du champ magnetique durant une inversion reste tres 
mal documente. Certaines mesures sur des coulees de lave qui auraient 
passe la Tq durant une inversion indiquent que le paleopole varie rapide- 
ment. Ces mesures suggerent que l’axe du dipole se deplace en quelques 
jours (et non que le champ magnetique terrestre s’annule avant de changer 
de signe). Cette interpretation reste controversee. Merrill et McFadden 
[1999] presentent une revue des donnees recentes a ce sujet. 


3.3 Anomalies magnetiques marines: Hypothese de 
Vine-Matthews-Morley. 

Les bandes paralleles d’anomalies magnetiques marines ont ete remarquees 
pour la premiere fois a proximite de la dorsale de Juan de Fuca, sur la 
cote ouest du Canada par Mason & Raff [1961]. Ce patron apparut se 
trouver le long de toutes les dorsales. A cause de 1’inclinaison du champ 
magnetique, la symetrie des anomalies apparait beaucoup mieux aux 
hautes latitudes (>45°) et lorsque les dorsales sont orientees NS. Leur 
amplitude est de l’ordre de 500nT. 

Apres Larry Morley, de la Commission Geologique du Canada, dont 
le papier n’a jamais ete publie, Vine et Matthews [1963] ont suggere que 
la croute oceanique se forme aux dorsales, acquiert sa magnetisation 
lorsqu’elle s’eloigne de la dorsale et se refroidit sous la temperature Tc- 
La croute oceanique enregistre done la polarite du champ magnetique 
au moment ou elle passe sous la temperature de Curie. La chronometrie 
des polarites du champ magnetique pernret done de determiner Page 
des differentes bandes de croute oceanique et le taux de formation de 
plancher oceanique. 

L’hypothese de Vine-Matthews-Morley a ete testee avec des donnees 
provenant de toutes les dorsales au cours des annees qui ont suivi. En- 
tre 1963 et 1967, la precision de la chronometrie s’est enormemenent 
amelioree. En particulier, l’identification de l’evenement Jaramillo, a 
perrnis de mieux modeliser quantitativement les anomalies magnetiques 
marines a proximite de toutes les dorsales (Fig 3.5) et de convaincre 
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Figure 3.3 Anomalies magnetiques a proximite de la dorsale de Reyk- 
janes (Islande). 



l’ensemble de la communaute scientifique de la validite de l’hypothese 
d’expansion des fonds oceaniques. 
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Figure 3.4 Hypothese de Vine Matthews Morley. Magnetisation de 
la croute oceanique. 



r 


3.4 Applications 

Par la suite, on a etendu la chronometrie des inversions du champ 
magnetique jusqua 200Ma. II est impossible de remonter plus loin dans 
le temps parce que le plancher oceanique est detruit et qu’il n’y a plus 
d’enregistrement continu du champ magnetique terrestre. 6 

II est done possible de dater directement le plancher oceanique en 
identifiant les anomalies magnetiques. La signature de chaque sequence 
d’anomalies magnetiques marines est suffisamment caracteristique pour 
que l’on puisse innnediatement l’identifier dans la chronometrie des in- 
versions (pour autant que le taux d’expansion des fonds oceaniques soit 
reste constant pendant la periode consideree). Ainsi, chaque bande par- 
allele d’anomalies magnetiques correspond a une isochrone. Les cartes 
de Page de la croute oceanique sont ainsi etablies a partir des anomalies 
magnetiques. Pour que cela soit possible, il faut supposer que le plancher 
oceanique ne se deforme pas, c.a.d. implicitement introduire le concept 
de plaque. 

Le ”taux d’expansion” aux dorsales peut etre determine directement 

6 Les mesures paleomagnetiques dans les continents ne permettent pas un 

echantillonnage suffisant pour construire cette chronologie. Toutefois, elles ont 
demontre que le champ magnetique terrestre a existe et s’est inverse depuis 
l’Archeen. 
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Figure 3.5 Exemple de modele expliquant les anomalies magnetiques 
marines. Les donnees (ligne continue) sont comparee avec l’anomalie 
calculee pour le modele de magnetisation totale de la croute 
oceanique. 



en mesurant Page du plancher oceanique en fonction de la distance a la 
dorsale. C’est done l’inverse de la pente de l’age mesure en fonction de 
la distance a la dorsale. C’est done le taux de production de plancher 
oceanique par flanc de dorsale. Ces taux d’expansion varient beaucoup 
d’une dorsale a l’autre. Les taux les plus rapides sont observes sur YEast 
Pacific Rise (>5cm/a) et les plus lents sont mesures a la dorsale de 
Lonronossov dans l’ocean arctique (<lcm/a). Les taux d’expansion dans 
l’Atlantique sont relativement faibles (l-2cm/a). Le taux d’expansion de 
la dorsale de Juan de Fuca a l’ouest de Vancouver est de ~ 3cm/a. 

Avec un taux d’expansion moyen de 2cm/a, les 35,000km de dorsales 
produisent done k 2x 3510 6 x 210 -2 = 1.4 x 10 6 nr 2 de nouveau plancher 
oceanique par an, ou encore, pour 6 km d’epaisseur de croute, les dorsales 
produisent ss 8.5 km 3 de nouvclle croute oceanique par an. 

Cette chronometrie permet ainsi de reconstruire l’histoire des oceans: 
pour reconstruire l’ocean Atlantique ou l’ocean Indien a 45Ma, il ’’suffit” 
de faire disparaitre tout le plancher oceanique plus jeune que 45Ma. Un 
exemple recent est donne par le papier de Roest et Srivastava [1989] qui 
reconstruisent a partir des anomalies magnetiques marines Pouverture 
de la mer du Labrador. 

Les forages du Deep Sea Drilling Program (DSDP auquel a succede 
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ODP) ont permis d’echantillonner la croute oceanique et de dater ces 
echantillons. Les ages obtenus ont confirme les ages predits a partir des 
anomalies magnetiques marines. 

L’hypothese d’expansion des fonds oceaniques perrnet aussi de mieux 
comprendre la nature de la croute oceanique. La couche sismique 1 corre- 
spond aux sediments et s’epaissit avec l’age. La couche 2 est une couche 
de roches volcaniques alterees par la circulation hydrothermale et la 
couche 3 est une couche de gabbros. Cette interpretation a ete confirmee 
par la decouvertes de lambeaux de croute et de manteau oceaniques 
transportes sur de la croute continentale. Ces roches oceaniques torment 
les sequences ophiolitiques. De telles suites se trouvent dans la plupart 
des chaines de montagnes jeunes (incluant les Appalaches a Terre Neuve 
et Thetford-Mines). Les ophiolites les plus complets et les mieux etudies 
sont ceux d’Oman et de Chypre. II taut noter que les ophiolites provien- 
nent probablement des bassins marginaux et non de veritables bassins 
oceaniques. 


3.5 Problemes en suspens 

L'hypothese d’expansion des fonds oceaniques explique bien le cadre 
tectonique des dorsales. L’origine des magmas qui alimentent le volcan- 
isme des dorsales (profondeur des sources, temperature dans le manteau) 
n’est pas entierement elucidee. Les tairx d’expansion varient enormement 
d’une dorsale a l’autre et la cause de ces variations n’est pas toujours 
bien comprise. 

Les taux d’expansion sont stables a l’echelle de quelques dizaines de 
Ma, mais ils varient a plus longue echelle. Par contre, il n’est pas evident 
que le tairx de production global de nouveau plancher oceanique ait varie 
a cette echelle de lOOMa. 

Du point du champ magnetique terrestre, on s’interesse de plus en 
plus aux variations d’intensite du champ magnetique terrestre qui ont 
causes des variations de la magnetisation de la croute oceanique. Ces 
variations de la magnetisation ont ete mesurees par des magnetometres 
sous marins a proximite de dorsale a taux d’expansion rapide. 
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Figure 3.6 Chronometrie des polarites du champ magnetique depuis 
165Ma 
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3.6 Notes supplementaires 

3.6 Notes supplementaires 

3.6.1 1: La methode sismique refraction 

Deux types d’ondes de deformation (ondes elastiques) peuvent se propager 
dans un solide: les ondes P (ondes de compression associees a un change- 
ment de volume) et les ondes S (ondes de cisaillement associees a un 
changement de forme). Les ondes P sont toujours plus rapides que les 
ondes S. Dans la croute, v p ss 6km s^ 1 et v s ss 3.5km s , dans le man- 
teau, v p > 8km s -1 et v s > 4.7km s -1 . 

Dans un milieu homogene (a vitesse constante), les ondes se propa- 
gent en ligne droite. A l’interface entre deux milieux homogenes, il y 
a reflexion et refraction des ondes. (II peut y avoir conversion partielle 
d’ondes P en S et reciproquement). Les lois de Snell determinent les 
directions de propagation des differentes ondes: 

Vpl Vp2 V 8 i V s2 -j\ 

sinipi sinz P 2 sinz s i sinz S 2 

Tangle i est Tangle entre le vecteur vitesse et la normale a Tinterface. 

Si V 2 > vi, quand Tangle d’incidence est tel que sin* > U 2 /U 1 , il ne 
peut pas y avoir d’onde refractee (sin *2 > 1). Dans ce cas, il y a reflexion 
totale. L’angle i c tel que sini c = vi/v^ est appele Tangle critic^ue. A 
Tangle critique, Tonde refractee suit Tinterface, (*2 = 7t/ 2). Une partie de 
cette onde est refractee de nouveau dans le premier milieu. La methode 
sismique refraction consiste essenticllement a observer et interpreter les 
temps d’arrivee de ces ondes refractees. 

Considerant une couche horizontale d’epaisseur h et vitesse vi sur un 
milieu de vitesse V 2 (V 2 > v\). Soit O une source d’ondes sismiques a la 
surface. A distance x de O, Tonde directe arrive au temps: 

t d = ~ (3.2) 

Vl 

et pour x > 2 htani c /v 2 , il y aura une onde refractee qui arrivera au 
temps: 


2 h x — 2htaxii r 

(3.3) 

t r — . + 

V\ COS l c V2 

x 2 h , uisin* CN 

= - + 1-- ~) 

(3.4) 

V2 V\ COS l c V2 

x 2hcosi c 

(3.5) 

= 1 
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Figure 3.7 Trajet de l’onde refractee 


O 


X 



Notez que des que, x > x c , ou: 

x c 2hcosi c x c 

— + < — 3.6 

V 2 Vi V\ 

l’onde refractee arrive avant l’onde directe. 

La methode sismique refraction consistait a observer le temps d’arrivee 
de la premiere onde en fonction de la distance. Pour x < x c , le temps 
d’arrivee est une droite de pente l/iq; pour x > x c , le temps d’arrivee 
est une droite de pente l/v 2 - On peut done directement determiner v± 
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et V 2 - En prolongeant la droite donnant le temps d’arrivee de l’onde 
refractee on determine l’ordonnee a l’origine to: 


to — 


2 h cos i c 


Vi 


a partir de laquelle on peut deduire l’epaisseur: 

Vl t 0 


h = 


2\A - vl/v \ 


(3.7) 


(3.8) 


Dans le cas de plusieurs couches horizontales dont la vitesse aug- 
rnente a chaque interface, le temps d’arrivee de l’onde refractee le long 
de l’interface K est donne par: 


t = 


K- i 0 , 

x v 2hsk cos Z]$ 

vk + hi vk 


(3.9) 


hk et Vk sont l’epaisseur et la vitesse de la couche k et sin ik = Vk/vx- 
En pratique, on peut done determiner les vitesses dans chaque strate 
en mesurant la pente des temps d’arrivee en fonction de la distance. 
L’intersection de chacune de ces droites avec l’axe x = 0, permet de 
calculer de fagon recursive l’epaisseur de chaque strate. 

Les methodes recentes d’interpretation consistent a interpreter non 
seulement le temps d’arrivee mais aussi l’amplitude de toutes les ondes 
enregistrees. 


3.6.2 Notes supplementaires: Flux de chaleur 

Le flux de chaleur mesure la quantite d’energie qui traverse une sur- 
face par unite de temps. La determination du flux de chaleur permet de 
determiner le bilan d’energie de la terre. Le flux se mesure en J m~ 2 s -1 = 
En pratique, on ne mesure que le flux conductif qui est proportionnel au 
gradient de temperature: 

dT 

, = (3.10) 

oil I\ est la conductivite thermique (W m _1 K -1 ) et le gradient de 
temperature est mesure en K.m _1 . Dans les publications anterieures a 
1970, le flux est mesure en h.f.u. (1 h.f.u. = /real cm -2 s _1 = 41.8mW nr - 

)■ 


W m 2 . 
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3.6.3 Notes supplementaires: Anomalies du champ 

magnet ique 


3.6.4 Champ d’un dipole. 


II n’existe pas de charge magnetique isolee (monopole) et les champs 
magnetiques sont deceits comrne causes par un couple de charges op- 
posees (dipole). Le moment du dipole m est un vecteur produit de la 
charge par la distance qui les separe. Le potentiel magnetique U est 
defini par: 

17 = m • V(-) (3.11) 

r 

Le champ magnetique est le gradient du potentiel. A distance R du 
dipole, le champ magnetique varie en fonction de Tangle 6 que fait le 
vecteur position avec l’axe du dipole. 


H r 


H e 


2m cos 9 
R 3 

in sin 9 
Z R 3 


(3.12) 

(3.13) 


pour les composantes radiale (i.e. verticale) et tangentielle (i.e. horizon- 
tale) du champ. En premiere approximation, le champ magnetique ter- 
restre est celui d’un dipole geocentrique, legerement incline par rapport 
a l’axe de rotation de la terre. 

Dans le systeme d’unite international (SI), l’unite de mesure du champ 
magnetique est A m -1 ; le moment d’un dipole est mesure en A m 2 Dans 
le vide, l’induction magnetique est proportionnelle au champ magnetique: 
B = /ioH oil po = 47t10 7 H m 1 est la permeabilite magnetique du vide. 
En general, B = po(H + M) , oil M est la magnetisation du milieu. En 
pratique, e’est l’induction qui est la quantite physique mesuree. L’unite 
de mesure de l’induction magnetique est le T (Tesla). La magnetisation 
(densite de moment dipolaire) a les memes unites que le champ. L’unite 
pratique en geophysique est le nT (nanoTesla, parfois encore appele 7 ). 
Au pole, l’intensite de l’induction magnetique terrestre est de l’ordre de 
60,000nT; a l’equateur, elle est de l’ordre de 30,000nT. 

Le vecteur champ magnetique est defini par son intensity (champ to- 
tal) et sa direction definie a partir de Tinclinaison (angle du champ 
magnetique avec l’horizontale) , et de la declinaison (angle entre la direc- 
tion du nord geographique et celle du champ magnetique horizontal, nord 
magnetique). L’inclinaison est liee a la latitude magnetique X = tt/2 — 6: 


tan/ = Hr/Hb = 2cot0 = 2tanA (3-14) 
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3.6.5 Proprietes magnetiques des roches et des 

mineraux. 

• Magnetisation remanente et induite. 

• Diamagnetisme, paramagnetisme, ferro, anti-ferromagnetisnie, ferri- 
magnetisme. 

• Mecanismes d’acquisition de la magnetisation remanente (thermique, 
visqueuse, chimique). 

La magnetisation des materiaux peut etre soit induite par un champ 
magnetique ambiant, soit remanente. La plupart des anomalies magnetiques 
dans la croute continentale sont causees par la magnetisation induite. 

La magnetisation induite est proportionnelle au champ: 

M = X H (3.15) 

B = /i 0 (l + x)H (3.16) 

ou la susceptibility magnetique x est une propriety physique des materiaux 
et n’a pas de dimension en unite SI. Parmi les mecanismes de magnetisation 
induite, il y a le paramagnetisme (dans la direction du champ ambiant) 
et le diamagnetisme (dans la direction opposee au champ ambiant). 

Les mecanismes les plus importants de magnetisation remanente sont 
le ferromagnetisme (les moments magnetiques des atomes sont paralleles) 
et le ferrimagnetisme (les moments ne sont antiparalleles mais leur sonirne 
n’est pas nulle). A la temperature de Curie, les proprietes magnetiques 
dispar aissent. Quand un mineral se refroidit et passe la temperature 
de Curie en presence d’un champ magnetique, il peut acquerir une 
magnetisation dans la direction du champ ambiant. La magnetisation 
ne varie plus des que la temperature est tombee sous la temperature de 
blocage. (Cette magnetisation thermoremanente (TRM) est le mecanisme 
le plus courant d’acquisition de magnetisation). La temperature de Curie, 
qui est de l’ordre de 600 °C, depend de la nature des mineraux. La 
temperature de blocage est inferieure de quelques dizaines de °C a la 
temperature de Curie. 


3.7 Notes supplementaires: Chronometrie K-Ar. 

La plupart des geochronometres isotopiques (U-Pb, Rb-Sr) n’etaient 
pas utilisables pour les roches volcaniques jeunes 7 . La chronometrie K- 

7 Ceci n’est plus vrai aujourd’hui et les spectrometres de masse modernes 
permettent une tres grande precision meme si les concentrations sont infimes. 
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Expansion des fonds oceaniques 


Ar utilise la disintegration du K, present dans pratiquenrent toutes les 
roches. Le 40 AT se desintegre soit en 40 Ar soit en 40 Ca. Le Ca est trop 
abondant pour etre utile. L 'Ar est un gaz inerte et est done en principe 
absent du reseau cristallin au moment de la formation du mineral. La 
constante de disintegration a iti diternrinie: A = 5.544 x 10 -lo a -1 
avec 89.5% vers le 40 Ca et 10.5% vers 40 Ar. La demi-vie est done 
t = 1.25 x 10 9 a. A partir du moment ou le cristal se ferme (t = 0) 
pour le K et l’Ar, on a: 

m K(t) = 40 K( 0) exp (—At) 

40 Ar(t ) = .105 4o A% 0) (1 - exp(-At)) 

^^ = .105(exp(At)-l) (3.17) 

qui pernret de dater le mineral. 

Techniquenrent, la mise au point fut difficile a cause de la contami- 
nation par l’Ar atmosphirique. Elle est corrigie en mesurant 36 Ar dans 
l’ichantillon qui est d’origine atmosphirique. 


exosup.com 


page facebook 


Figure 3.8 Carte des ages du plancher oceanique d’apres Muller et al. 
(2008) 



Age (Myr) 


CO 


exosup.com 


page facebook 


3. 1 Notes supplementaires: Chronometrie K-Ar. 


4 

Les plaques: Cinematique 


L’hypothese de Vine-Matthews-Morley d’expansion des fonds oceaniques 
ne concerne que la formation de plancher oceanique aux dorsales. La 
dorsale ou centre d’expansion est un premier type de frontiere ou il y a 
divergence entre plaques. L’expansion des fonds oceaniques n’impliquait 
pas necessairement la tectonique des plaques. Ainsi, d’autres hypotheses 
ont egalement ete envisagees. Certains geologues (Carey, Heezen) et des 
physiciens ont considere tres serieusement que la Terre est en expan- 
sion. Pour en venir a la tectonique des plaques, il a fallut d’une part 
reconnaitre le role des autres frontieres entre plaques: les zones de sub- 
duction et les failles transformantes, d’autre part demontrer l’absence 
de deformations a l’interieur des plaques. Pour les oceans, les anomalies 
magnetiques marines qui restent rectilignes et paralleles quel que soit 
leur age, demontrent clairement que les plaques oceaniques ne sont pas 
deformees. 


4.1 Failles transformantes 

Les failles qui joignent des segments de dorsales medioceaniques appa- 
raissent comrne de grandes failles de decrochement. Cette interpretation 
se heurte a une difficulte parce que le mouvement de ces failles s’evanouit 
sans qu’il y ait trace de deformation a leurs extremites. 

C’est Wilson (1965) qui a identifie le probleme geometrique et pro- 
pose que les deplacements entre ’’plaques” le long de failles transfor- 
mantes etaient limites a la region comprise entre deux autres frontieres 
de plaque. (Fig 4.1). Les observations simiques ont confirme cette in- 
terpretation. La sismicite est concentree sur la partie de la faille trans- 
formante entre les deux dorsales. On verra egalement que Ton peut 
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Figure 4.1 Failles transformantes entre deux dorsales ou entre deux 
zones de subduction. Le mouvement le long de la faille doit etre limite 
au segment T\ Tb- La polarite de la subduction est indiquee par la 
concavite de l’arc insulaire. L’arc est concave du cote de la plaque 
qui chevauche. 



determiner le mouvement le long de la faille a partir des observations 
sismiques qui ont ainsi confirme l’hypothese de Wilson. 

Wilson a egalement reconnu les differents types de failles transfor- 
mantes entre deux dorsales, entre deux arcs insulaires, ou entre dorsale 
et arc insulaire et a montre comment evoluerait chacun de ces systemes. 
La faille est dite transformante parce que la nature du cisaillement le long 
de la faille est transforme soudainement. Cette hypothese est confirmee 
par la concentration de l’activite sismique le long de la faille transfor- 
mante entre T\ et sur la figure 4.1 et par les mecanismes au foyer des 
tremblements de terre qui ont confirme la nature du mouvement le long 
de la faille transformante. 
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4.2 Zones de subduction 


L ’importance et role des fosses profonds avaient ete reconnus grace 
aux etudes gravimetriques de Vening-Meinesz dans les annees 1940. 
En effet c’est a proximite des zones de subduction que Vening-Meinesz 
avait decouvert les anomalies d’air libre les plus fortes au monde. II en 
avait deduit que ces anomalies refletent une absence d’equilibre due aux 
courants de convection dans le manteau. La concentration de l’activite 
sisnrique et volcanique autour du Pacifique avait egalement ete observee 
depuis longtemps. Les arcs insulaires sont associes aux zones d’activite 
sisnrique profonde, zones dites de Wadati-Benioff. La profondeur de 
l’activite sismique augmente avec la distance au fosse profond. La pro- 
fondeur maximum des trenrblements de terre est de 600knr. Le pendage 
des zones de Wadati-Benioff varie d’une zone a l’autre et change egalement 
a l’interieur d’une zone; il est compris entre 25 °et plus de 60 Avant la 
tectonique des plaques, ^interpretation la plus courante etait que cette 
sismicite etait due aux contraintes dues aux differences de densite et de 
regime thernrique entre oceans et continents. Dans le papier de Wilson 
(1965), les arcs insulaires sont inrplicitement reconnus comme zones de 
convergence entre plaques et de destruction d’une des deux plaques. 
Tuzo Wilson avait insiste sur la concavite des arcs insulaires et fait 
l’analogie avec la surface d’une balle de ping-pong coupee et poussee 
vers l’interieur. 


Les frontieres de plaque, oil se concentre toute l’activite tectonique, sont 
definies par l’activite sismique. Les anomalies nragnetiques marines qui 
restent paralleles montrent que le plancher oceanique, meme ancien, 
n’est pas deforme entre la dorsale et les fosses profonds. La concentration 
de 95% de l’activite sismique dans aux frontieres de plaque montre aussi 
que c’est la qu’est concentree la plus grande partie des deformations 
actuels. Les mecanismes aux foyers des tremblements de terre permet- 
tent de determiner la direction du mouvement relatif des plaques. 

Ainsi, il y a trois types de frontieres entre plaques (Figure 4.2: frontieres 
divergentes (centres d’expansion ou dorsales) oil se forme la lithosphere 
oceanique, frontieres convergentes (zones de subduction) oil la lithosphere 
oceanique est detruite, et failles transformantes oil le mouvement relatif 
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Figure 4.2 Les trois types de frontieres entre plaques. 



de deux plaques se fait parallelement a leur frontiere et ou il n’y a ni 
formation ni destruction de lithosphere. 

La tectonique des plaques a eu pour effet de focaliser l’interet des 
tectoniciens sur les frontieres de plaque et d’exclure tous les autres 
phenomenes en supposant que les plaques sont rigides et ne se deferment 
pas. A l’origine, la tectonique des plaques expliquait exclusivement revolution 
de la lithosphere oceanique ou ellc s’applique aisement et avec succes. 

On reconnaissait bien qu’il y avait differentes possibility de convergence 
(ocean-ocean; ocean-continent; continent-continent). Ce n’est que pro- 
gressivement que la theorie a evolue pour inclure les deformations plus 
profondes a l’interieur des plaques ainsi que la possibilite de frontieres 
plus diffuses (surtout dans les continents). II est clair maintenant que, 
dans le cas de la collision de l’lnde avec l’Asie, les deformations sont tres 
profondes et concernent tout le continent asiatique. Par exemple, il sem- 
ble maintenant que la plaque australo-indienne consiste en deux ou trois 
plaques avec frontieres diffuses (Royer et ah, 1997). La Mediterrannee a 
ete divisee en une douzaine de microplaques. Le mecanisme de la colli- 
sion continentale reste assez mal compris. Le role de la subduction de 
croute continentale reste un probleme controverse. 
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Figure 4.3 L’activite sismique entre 1975 et 1995. L’activite sismique 
est limitee aux frontieres de plaques. C’est d’ailleurs ainsi que sont 
definies ces frontieres. 
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Figure 4.4 Les frontieres entre les plaques principales dans les oceans. 


A l’interieur des continents les frontieres de plaques sont en general 
plus diffuses que dans les oceans. On a egalement identifie de nombreuses 
microplaques en particulier dans les zones de collision. Ainsi, la tec- 
tonique et l’activite sismique en Mediterrannee demandent une douzaine 
de microplaques. L’Inde et l’Australie font nraintenant partie de plaques 
distinctes avec une frontiere diffuse. Le rift est-africain separe la plaque 
africaine de la plaque Somalienne. 

La compilation de Bird (2003) recense un total de 52 plaques et micro- 
plaques (voir Figure 4.5). 
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Figure 4.5 Les front ieres entre les plaques principales d’apres Bird 
(2003). Notez que dans ce modele le nombre de plaques et micro- 
plaques est passe a 52. Les regions ou la deformation est diffuse sont 
hachurees. 


00 4a 20* a 20* 40* eo* 



o 

CM 


s 

8 


o 

CM 


o 

T 


o 

ID 


O 

00 


o 

(O 

9 

S 

8 

3 

S 

? 

8 


exosup.com 


page facebook 



4-4 Cinematique des plaques. 

4.4 Cinematique des plaques. 


51 


Une des applications directes de la tectonique des plaques concerne 
Pactivite aux frontieres des plaques qui a pu etre expliquee grace a des 
contraintes geometriques tres simples. 

4.4.1 Description du mouvement dans l’espace des 


S’il est possible de determiner la vitesse relative de deux plaques, il est 
impossible de definir leur vitesse absolue. Pour cela, il faudrait disposer 
d’un systeme de reference fixe par rapport auquel on pourrait mesurer 
la vitesse. 

Il est utile de representer le mouvement des plaques dans l’espace des 
vitesses, ou la vitesse est un vecteur. On peut arbitrairement clioisir 
de mesurer toutes les vitesses relativement a celle de la plaque A (la 
vitesse de A est nulle). Le vecteur avb represente la vitesse de de B 
relativement a A. Si on considere plusieurs plaques, il faut que les vitesses 
s’additionnent: 


Pour les dorsales, la composante de la vitesse perpendiculaire a la 
dorsale peut etre calculee a partir du taux d’expansion mesure par les 
anomalies magnetiques marines. Toutefois, la vitesse n’est pas necessairement 
normale a la dorsale. La seule contrainte sur la direction du vecteur 
vitesse provient des failles transformantes. La vitesse doit etre parallele 
a la faille transformante puisqu’il n’y a ni creation, ni destruction de 
lithosphere. Dans le cas de dorsales ou de zones de subduction, la vitesse 
n’est pas necessairement perpendiculaire a la frontiere des plaques. Dans 
la plupart des cas, la subduction est oblique. Par contre, la plupart des 
failles transformantes entre dorsales sont a peu pres perpendiculaires 
aux dorsales, mais il y a des exceptions. 

Jusqu’a tres recemment, seuls les taux d’expansion aux dorsales pou- 
vaient etre mesures a partir des anomalies magnetiques marines. On en 
deduisait la vitesse dans la direction parallele aux failles transformantes 
et la vitesse de convergence pour la subduction. On peut maintenant 
mesurer directement la vitesse relative instantanee des plaques a l’aide de 
la geodesie spatiale (GPS, VLBI). Il semble cependant que cette vitesse 
instantanee soit differente de la vitesse moyenne a l’echelle geologique. 


vitesses 


AVb + BVC = AVC 
AVB + BVC + CVA = 0 


(4.1) 

(4.2) 
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Figure 4.6 Exemple de frontieres de plaques et de leur espace de 
vitesses. Les unites de distance et de vitesse sont arbitraires. 

espace des vitesses 3 dorsales 




espace des vitesses 


1 dorsale & 2 subductions 



Ainsi, sur la faille de San Andreas il y a discordance entre les vitesses 
mesurees aujourd’hui par la geodesie et les vitesses geologiques. 


4.4.2 Mouvement sur une sphere. Poles d’Euler 

Les contraintes de geometrie plane ne sont qu’approximatives car il faut 
tenir compte de la sphericite de la terre. Ainsi, le mouvement relatif 
entre deux plaques rigides sur la surface d’une sphere est defini comrne 
une rotation par rapport a un pole fixe (theoreme d’Euler). Dans cette 
geometrie les failles transformantes sont des par alleles, la vitesse rela- 
tive entre deux plaques, le long d’une faille transformante, depend de la 
distance au pole de rotation, i.e. de la colatitude <j>: 

v t =nxRx sin(^) (4-3) 

ou Q est la vitesse angulaire de rotation. 
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Le pole de rotation est determine par l’intersection des arcs de grand 
cercle perpendiculaires aux failles transformantes. Les failles transfor- 
rnantes entre la plaque nord americaine et la plaque pacifique (San 
Andreas en Californie, Queen Charlotte en Colombie britanique) don- 
nent un pole de rotation situe pres de Terre Neuve. L’ouverture de 
l’Atlantique sud correspond a une rotation dont le pole se trouve a 
proximite de l’lslande. Morgan (1968) a ainsi rnontre que la vitesse 
d’ouverture de l’Atlantique sud (mesuree a partir des anomalies magnetiques 
marines) dependait bien de la distance au pole suivant l’equation 4.3. 

La demonstration de Morgan etait non seulement elegante, mais surtout 
cllc montrait qu’il etait possible de comparer les predictions theoriques 
avec les donnees geophysiques. 


4.4.3 Lois de conservation 

Bien que la composition de deux rotations n’est pas aussi simple que 
l’addition de vitesse dans la plan, la loi d’addition des rotations est 
simple. Soient trois plaques A, B, C et est la rotation de B relativement 
a A , avec la condition a^b = —b^a, il faut: 

a&b +b &c Ac &A = 0 (4.4) 


4.4.4 Stability des points triples 

Un point ou se rencontrent trois plaques distinctes est un point triple. 
Des conditions de compatibility entre les vitesses des plaques doivent 
etre verifiees pour qu’un point triple soit stable. La plupart des points 
triples sont instables. Pour qu’un point triple soit stable, il faut pouvoir 
determiner sa vitesse dans l’espace des vitesses. Cette vitesse doit etre 
celle d’un point frontiere entre plaques A et B 1 et C , et C et A. 

Pour qu’un point reste sur une dorsale, il faut que sa vitesse soit 
parallele a la dorsale et a mi chemin entre les vitesses des deux plaques. 
Pour qu’un point reste sur une faille transformante, il faut que sa vitesse 
reste parallele a la plaque transformante. Pour qu’un point reste sur une 
zone de subduction, il faut qu’il reste sur la plaque qui passe au dessus; 
il peut se deplacer parallelement a la frontiere. Sa vitesse est celle de la 
plaque superieure plus une composante arbitraire parallele a la zone de 
subduction. 

Ainsi, supposons trois dorsales et leur triangle de vitesse, le point 
triple doit se trouver a l’intersection des trois mediatrices. Parce que les 


exosup.com 


page facebook 


54 


Les plaques: Cinematique 


trois mediatrices s’intersectent toujours en un point, ce point definit la 
vitesse du point triple. Pour trois failles transformantes, les vitesses sont 
paralleles aux failles. Un point sur la frontiere entre B et A reste sur 
la frontiere si sa vitesse est sur la droite AB dans l’espace des vitesses, 
de meme pour les autres frontieres. II n’y a done pas dans l’espace des 
vitesses d’intersection entre ces trois droites. II y a done incompatibility?, 
et une triple jonction entre 3 failles transformantes est toujours instable. 
Pour une zone de subduction, un point sur la frontiere doit rester avec la 
plaque qui passe au dessus, il doit done avoir la vitesse de cette plaque, 
ou une vitesse parallele a la zone de subduction. 

II faut remarquer que le type de point triple le plus communement 
observe est celui qui est inconditionnellement stable (jonction de trois 
dorsales). 


4.5 Applications 

Les contraintes geometriques sur le mouvement des plaques permettent 
de mieux interpreter le contexte tectonique d’une region. Par exemple, 
les nrodeles devolution tectonique et le regime actuel de la cote Pacifique 
de l’Anrerique du Nord ont ete reveles par de simples applications de la 
tectonique des plaques. 

4.5.1 Tectonique de la cote Pacifique de l’Amerique du 

Nord 

L’age des anomalies magnetiques marines a 1’ouest de la cote californi- 
enne augnrente vers l’ouest. Ceci indique qu’une dorsale etait active en- 
tre la cote californienne et ces anomalies. L’age de la plus jeune de ces 
anomalies donne l’age maximum de la disparition de la dorsale qui a ete 
subductee sous la Californie. Entre la Californie et cette dorsale, se trou- 
vait la plaque de Farallon. Depuis la disparition de la dorsale, la faille de 
San- Andreas nrarque la frontiere entre la plaque Pacifique et la plaque 
nord americaine. La faille transfornrante de San Andreas et la faille trans- 
formante de Mendocino se joignent au Nord de SF. Le long de la cote 
au Nord de la Californie, et en Oregon, l’activite sismique est faible 
et il n’y a pas evidence directe qu’il y a subduction. Deux hypotheses 
sont ainsi possible: (1) Il existe une mini plaque entre l’Anrerique du 
Nord et la dorsale de Juan de Fuca, cette miniplaque est subductee sous 
l’Oregon (2) La dorsale de Juan de Fuca est la frontiere entre la plaque 
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Figure 4.7 Condition de stability de jonctions triples. Trois dorsales. 
La jonction triple est toujours stable 



Figure 4.8 Condition de stability de jonctions triples. Trois failles 
transformantes, la jonction triple n’est jamais stable. 



Figure 4.9 Condition de stability de jonctions triples. Trois zones de 
subduction. La jonction triple n’est jamais stable. 
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nord americaine et la plaque pacifique. Toutefois, cela est incompatible 
avec les vitesses. Pour eliminer cette contradiction il faut introduire une 
mini plaque entre la dorsale de Juan de Fuca et l’Amerique du Nord. 
La vitesse de cette plaque est determinee par le triangle de vitesse et 
implique qu’il y a subduction oblique sous l’Oregon. Le point triple au 
nord de la faille de San Andreas est stable. 
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Figure 4.10 Les frontieres de plaques reconnues sur la cote pacifique. 
L ’introduction d’une miniplaque F est necessaire pour que les vitesses 
soient compatibles. La resolution des vitesses montre qu’il y a sub- 
duction oblique de cette plaque sous l’Oregon. 
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L’analyse des vitesses a done ete utile pour etablir la nature des 
frontieres des plaques sur la cote ouest. En analysant la stability des 
deux points triples a chaque extremite de la faille de San Andreas, on 
peut montrer que ces deux points s’eloignent et done que la faille de 
San Andreas grandit. Aux taux actuels, cela impliquerait que la faille de 
San Andreas a commence son existence a 20 Ma. Ceci est aussi confirme 
par les anomalies magnetiques marines. Le plancher oceanique a l’ouest 
de la Californie est relativement jeune et son age augmente vers l’ouest, 
cad que la dorsale dont il provient doit se trouver a Test. La plaque 
qui se trouvait a l’est de cette dorsale, Plaque de Farallon , et la dorsale 
elle-meme ont ete subductees sous la Californie a 20 Ma. 
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Figure 4.11 La situation actuelle sur la cote ouest de l’Amerique du 
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Figure 4.12 Les anomalies magnetiques et l’age du plancher 
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Figure 4.13 Evolution de la cote ouest et de la faille de San Andreas 
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L’ analyse des mesures recentes du mouvement des plaques par VLBI 
(very long baseline interferometry) montre que le mouvement est encore 
plus complique et il necessaire d’introduire unc microplaque supplemental 
pour absorber une partie du mouvement entre la faille de San Andreas 
et la plaque Nord Americaine. 

Bien d’autres problemes de tectonique contemporaine peuvent ainsi 
etre analyse a l’aide des contraintes geometriques de la tectonique des 
plaques. C’est le cas des frontieres de plaques autour de la Nouvelle 
Zelande ou des regimes tectoniques de la Turquie et du Moyen-Orient. 
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Point Chauds-Cycle de Wilson 


5.1 Points Chauds 

L’activite volcanique a l’interieur des plaques est anormale du point de 
vue de la tectonique des plaques. Wilson (1963) avait remarque que les 
iles Hawaii, formees par des volcans actifs aujourd’hui, faisaient partie 
d’une chaine rectiligne d’iles volcaniques qui se prolonge vers le Nord 
Quest du Pacifique. 


Tuzo Wilson avait suggere que cette activite volcanique intraplaque 


Figure 5.1 La chaine de Hawaii dans locean Pacifique 
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Figure 5.2 La trace des points chauds 


FU*d 
Mol Spot 


est causee par la presence de panaches (plumes) chauds provenant du 
manteau inferieur et qui penetrent la lithosphere. Si ces points chauds 
sont fixes dans le manteau, ils fournissent un systeme de reference pour 
mesurer le mouvement des plaques par rapport au manteau. Ainsi, la 
chaine des volcans d’Hawaii a enregistre le mouvement de la plaque 
Pacifique par rapport a un point chaud fixe. A partir d’Hawaii, la trace 
du point chaud se prolonge vers le WNW puis tourne soudainement vers 
le Nord. Les iles volcaniques a hauteur du Kamchatka ont un age de 
80Ma. Les iles a hauteur de la cassure ont un age de 45 Ma. II semble done 
que la plaque Pacifique se dirigeait vers le Nord jusqu’a 45 Ma et qu’elle 
a change de direction. Ce changement correspond a une reorganisation 
globale du mouvement des plaques. Apparement, cette reorganisation 
s’est faite tres rapidement. 

Par la suite, ? a suggere que la convection dans le manteau se fait 
a travers ces panaches qui transportent la chaleur depuis le manteau 
inferieur jusqu’a la base de la lithosphere. 
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5.1.1 Points chauds actifs 

On a repertorie une centaine de points chauds actifs aujourd’hui. Les 
plus importants sont Hawaii, l’lslande, les iles de la Reunion, les Afars. 
A l’interieur de l’Amerique du Nord, Yellowstone est un point chaud. 
Les points chauds sont plus frequents dans les oceans que dans les con- 
tinents. L’Afrique oil il y a de nombreux points chauds fait exception. 
Parce que la croute oceanique est recyclee, il y a peu de traces de points 
chauds au dela de 200Ma. Les Monteregiennes au Quebec et en Nou- 
velle Angleterre sont la trace d’un point chaud. Les essaims de dykes 
du McKenzie (McKenzie dykes swarm) qui partent de l’embouchure de 
la riviere McKenzie et traversent le bouclier canadien jusqu’au sud de 
la baie d’Hudson sont probablement dus a Taction d’un point chaud a 
l,200Ma. 

La lithosphere continental est supposee etre trop epaisse pour per- 
mettre a un point chaud de percer, sauf si le continent est immobile, 
ce qui est le cas du continent africain. Il faut noter que le volcanisme 
du point chaud de Hawaii perce la plaque Pacifique bien que celle-ci se 
deplace tres rapidement. 


5.1.2 Fixite des points chauds 

L’hypotliese que les points chauds sont fixes dans le manteau suppose 
que leur mouvement relatif est nul. Cela peut etre verifie en comparant 
les traces des points chauds. Il semble que cela soit approximativement 
vrai jusqu’a 120 Ma. Dans ce cas, les points chauds fixes forment un 
systeme de reference pour mesurer le mouvement absolu des plaques. 
Minster & Jordan [1978] ont ainsi determine le mouvement absolu de 
toutes les plaques dans un systeme de reference ou les points chauds 
sont fixes. Argus & Gordon [1991] ont montre que ce syterne de reference 
des points chauds est different du systeme dans lequel la resultante de 
toutes les rotations est nulle. Ceci inrplique que l’ensemble des plaques se 
deplace legerement vers l’ouest par rapport au manteau. Ce deplacement 
est peut-etre cause par l’orientation dominante des zones de subduction 
qui tirent les plaques vers l’ouest. 


5.1.3 Application: Ouverture de l’Atlantique 

L’ouverture de l’ocean Atlantique a ete precedee par une intense activite 
de points chauds. Les traces des points chauds peuvent etre utilisees pour 
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reconstruire le mouvement des plaques. Les resultats sont compatibles 
avec les anomalies magnetiques marines [ Morgan , 1983]. 
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EXPLANATION 

Divergent plate boundaries— 
Where new crust is generated 
as the plates pull away from 
each other 

Convergent plate boundaries — 
Where crust is consumed in the 
Earth's interior as one plate 
dives under another. 

Transform plate boundaries— 
Where crust is neither produced 
nor destroyed as plates slide 
horizontally past each other. 

Plate boundary zones — Broad 
belts in which deformation is 
diffuse and boundaries are not 
well defined. 

Selected prominent hotspots 
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J. Tuzo Wilson avait propose que la formation des Appalaches etait 
le resultat d’un cycle d’ouverture et de fermeture du proto-ocean at- 
lantique, Iapetus (Wilson, 1963). C’etait la premiere veritable applica- 
tion de l’expansion des fonds oceaniques a la geologie. Cela impliquait 
aussi que la tectonique des plaques avait fonctionne pendant une grande 
partie de l’histoire de la Terre et que la formation de super-continents 
et l’ouverture de bassins oceaniques font partie de cycles qui se repro- 
duisent. C’est Kevin Burke qui a donne le nom de cycle de Wilson au 
cycle d’ouverture et fermeture d’un bassin oceanique. 

Le cycle de Wilson , qui place l’orogenese dans le cadre de la tectonique 
des plaques, comprend les etapes suivantes: 

• Formation d’un rift intracontinental. II y a extension de la croute con- 
tinentale, formation d’un graben et activite volcanique. Aujourd’hui 
le rift est-africain et le rift du Rio Grande en sont des exemples. 
L’ouverture d’un bassin oceanique ne suit pas toujours la formation du 
rift intracontinental. II existe de nombreux exemples de rifts avortes a 
l’interieur des continents (graben du St-Laurent, rift de Keeweenawan) . 

• Debut d’ouverture d’un ocean. A ce stade, il y a formation de croute 
oceanique a partir d’un rift central. La Mer Rouge est un exemple de 
bassin oceanique en formation. 

• Evolution d’un bassin oceanique. Le plancher oceanique continue de 
se former a partir de la dorsale. les marges continentales, dont la 
croute a ete amincie, s’enfoncent et permettent l’accumulation des 
sediments. Les marges continentales jouent le role des geosynclinaux. 
Cette situation correspond a celle de l’ocean Atlantique. 

• Initiation de la subduction et destruction du plancher oceanique. For- 
mation d’iles en arc par subduction ocean-ocean. Subduction ocean- 
continent et formation de marges continentales actives. L’ocean Paci- 
fique contient tous ces elements. 

• Fermeture d’un ocean avec collision continentale. L’exemple est la 
collision de l’Inde avec l’Asie qui forme l’Himalaya. 

Ce resume tres simplifie ne rendaient evidement pas compte de tous les 
phenomenes. Par ailleurs, certaines etapes du cycle peuvent manquer. 
Mais le cycle de Wilson met la plupart des phenomenes geologiques dans 
un contexte global. 
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Les frontieres de plaque, ou se concentre toute l’activite tectonique, sont 
definies par 1‘activite sismique. La nature de cette activite permet de 
determiner le mouvement relatif des plaques. Avec les ondes sismiques, 
il est possible de determiner le mouvement des failles et d’en deduire le 
mouvement relatif de deux plaques. 


6.1 Observations des tremblements de Terre 

Les ondes sismiques sont des ondes de deformation qui se propagent dans 
les solides. La deformation associee aus ondes P est un changement de 
volume. Le passage de l’onde cause une compression suivie de dilata- 
tion, ou une dilatation suivie de compression. La deformation associee 
aux ondes P est un changement de forme sans changement de volume 
(cisaillement). La direction du premier mouvement des ondes sismiques 
permet de determiner le plan de rupture et la direction du mouvement 
et des contraintes qui ont cause la rupture sur la faille. 

Les ondes sismiques produites par les tremblements de terre sont ob- 
servees a l’aide de sismographes qui enregistrent les accelerations du sol. 
Les observatoires enregistrent les trois composantes de l’acceleration (en 
general dans deux gammes de frequences: courtes periodes ss Is, longues 
periodes « 40s) La difference des temps d’arrivee entre les ondes P et S 
permet de determiner la distance du tremblement de Terre. Avec suff- 
isamment d’observations, on peut localiser le foyer. Le foyer designe le 
point precis ou a eu lieu le tremblement de terre. L ’epicentre designe sa 
projection a la surface. 
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Echelle de Richter 

Pour quantifier les tremblements de terre, on utilise la magnitude (echelle 
de Richter). II s’agit d’une echelle logarithmique qui mesure l’amplitude 
du mouvement vertical du sol enregistre lors du passage d’une onde de 
frequence donnee a une certaine distance de l’epicentre. II y a plusieurs 
definitions de la magnitude dependant de l’onde choisie. L’energie du 
tremblement de terre est liee a l’amplitude. Empiriquement, on a con- 
state que l’energie etait multipliee par 30 quand la magnitude augmen- 
tait d’une unite. 1 La relation entre l’energie E s et la magnitude M s est 
donnee par une relation de Gutenberg-Richter: 

Logi 0 E s = 1.5M s +4.8 (6.1) 

Dans cette relation, l’energie est calculee en Joules. Un tremblement de 
Terre de magnitude 8.5 libere une energie > 10 17 J. l’energie totale liberee 
par toute l’activite sismique mondiale est de l’ordre de > 10 18 J/a. 

Moment sismique 

Les sismologues utilisent maintenant la notion de moment d’un trem- 
blement de terre qui depend du deplacement de la faille et de la surface 
totale de la zone de rupture. 

Le moment sismique est un tenseur dont les composantes M tJ sont: 

Mij = p(Auinj + AujUi)A ( 6 . 2 ) 

ou /i est le module de cisaillement, Atq est la composant i du deplacement, 
nj est la normale au plan de la faille, et A est la surface de la zone de 
rupture. Le moment se mesure en N m. 


6.2 Sismicite et definition des frontieres de plaques 

La distribution geographique des tremblements de terre suit un pa- 
tron bien defini (Figure 4.3). Les frontieres de plaque peuvent ainsi etre 
definies a partir des epicentres des seismes. 

II y a plusieurs milliers de tremblements de terre enregistres chaque 
annee. En moyenne, il y a un tremblement de terre de magnitude > 8 par 
an. Un tremblement de terre de magnitude > 8 degage autant d’energie 
que l’ensemble de tous les tremblements de l’annee. Les tremblements 
de terre de magnitude >8.5 sont exceptionnels. 

1 Theoriquement, l’energie est proportionnelle au carre de l’amplitude et le facteur 
devrait etre 100 au lieu de 30. 
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A proximite des dorsales, la plus grande partie de Pactivite est liee 
aux failles transformantes. II y a aussi un pen d’activite sur les Hanes de 
la dorsale. Ces tremblements de terre sont superficiels, et leur magnitude 
est en general <6.5. 

Les tremblements de terre sur les failles transformantes qui ne sont 
pas directement attachees a une dorsale (San Andreas) peuvent etre plus 
violents. 

Les tremblements de terre les plus violents sont lies a la subduction 
et a la collision continentale. Ce sont aussi les seuls tremblements de 
terre profonds (< 100km pour la collision continentale, < 600knr pour 
la subduction). 

Les donnees instrumentales de l’activite sismique sont tres recentes 
(salOO ans). En Amerique du Nord, les donnees historiques ne couvrent 
que quelques centaines d’annees. Ces donnees sont tres insuffisantes pour 
en deduire des tendances a long terrne. La tectonique des plaques a ainsi 
attire l’attention sur des regions ou il y a un deficit de Pactivite sismique 
(Oregon, lie de Vancouver). Dans la plupart de ces regions, des etudes 
de paleosismicite ont confirme que le risque sismique a long terme est 
eleve. 


6.3 Mecanisme au foyer des tremblements de terre 


Rappel: Geometrie des failles. Inclinaison et declinaison. Rejet. Con- 
traintes principales. Types de failles. 

Rappel: Representation a l’aide de projection stereographiques du 
plan d’une faille et du mouvement. 

6.3.2 Determination des mecanismes au foyer 

Un tremblement de Terre est cause par un mouvement soudain le long 
d’un plan de faille. La friction sur le plan de faille empeche tout mouve- 
ment de se produire et les contraintes s’accumulent autour de la faille. 
Lorsque la contrainte depasse le seuil de friction (loi de Byerlee), il y a 
un tremblement de Terre avec liberation brutale de Penergie accumulee. 
On peut montrer que dans un solide, la direction de la contrainte de 
cisaillement maximum est a 45 °entre les contraintes principales max- 
imum et minimum. Il y a done deux directions perpendiculaires pour 
lesquelles la contrainte de cisaillement est maximum. 


6.3.1 Description des failles 
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Figure 6.1 Mouvement d’une faille. Premiers mouvements des ondes 
P. Les zones ombragees sont celles ou le premier mouvement de l’onde 
P est une compresion. 






S'CHMAl. FAVI-T 




Pour qu’il y ait mouvement sur une faille, il faut que la contrainte 
de cisaillement soit superieure a la friction. On peut representer la con- 
trainte de cisaillement par un couple de forces paralleles au plan de la 
faille. Pour qu’il n’y ait pas de rotation, ce couple doit etre contrebal- 
ance par un couple de direction opposee. Ce couple agit sur le plan 
perpendiculaire a la direction du cisaillement dans le plan de faille. Ce 
deuxieme plan est dit plan conjugue. Les axes principaux, T de tension 
maximum, et P de compression maximum, sont a 45 °entre les deux 
plans conjugues. Le plan de faille et le plan conjugue divisent la terre 
autour de la faille en 4 quadrants. Lorsqu’il y a rupture de la faille, le 
premier mouvement dans chacun de ces quadrants sera soit une com- 
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Figure 6.2 Mouvement d’une faille et premiers mouvements des ondes 
P 


pression, soit une dilatation. Ce premier mouvement est celui des ondes 
P. II faut noter que ce premier mouvement est une compression dans les 
quadrants qui sont en tension (comme lorsqu’un ressort tendu casse, il se 
contracte) . Dans les regions en compression, le premier mouvement sera 
une dilatation (comme un ressort qui est comprime se detend lorsqu’il 
n’est plus retenu). 

En determinant la distribution geographique de la direction du pre- 
mier mouvement des ondes P, il est possible de determiner le plan de 
faille et le plan conjugue ainsi que les axes principaux des contraintes. 
S’il n’y a aucune ambiguite pour les axes T et P, il est par contre impos- 
sible de reconnaitre a partir des donnees sismiques, le plan de faille du 
plan conjugue. Ces deux plans sont appeles plans nodaux. L’ambiguite 
peut etre resolue si la direction des failles est connue par la geologie. 
la determination des mecanismes au foyer est faite automatiquement a 
l’aide d’ordinateurs. Elle peut se faire facilement en utilisant les pro- 
jections stereographiques. On represente en projection, la demi-sphere 
inferieure et on y reporte les premiers mouvements des ondes P. En 
principe, la polarite du premier mouvement permet de determiner la 
trace des deux plans conjugues. On peut determiner les axes P et T de 
compression et tension maximum. Si la geologie permet de determiner 
le plan de la faille, on peut alors determiner la direction du mouvement. 

6.3.3 Application 

Les mecanismes au foyer ont demontre la nature du mouvement le long 
des failles transformantes. Si la discontinuite de la dorsale etait due a 
un decrochement, le mouvement de la faille serait un mouvement lateral 
gauche. Par contre, le mouvement de la faille transformante est un mou- 
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vement lateral droit. Dans le premier cas, les quadrants 1 et 3 enreg- 
istrent une compression et les quadrants 2 et 4, line dilatation. Sinon 
c’est l’inverse. Les etudes de la dorsale mid Atlantique ont clemontre 
que l’interpretation de Tuzo Wilson etait correcte. 
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Figure 6.3 Contraintes principales et mecanismes au foyer pour les 
differents types de failles. En general les mecanismes sont plus com- 
plexes avec une composante de decrochement superposee a la com- 
posante normale ou inverse. 
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Figure 6.4 Faille transformante. S’il s’agissait d’une faille de 
decrochement classique, les mecanismes au foyer seraient inverses 
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Figure 6.5 Un tremblement de terre le long de la dorsale du Chili. Si le 
plan de faille est EW, le mouvement du bloc Nord est vers l’Est. Pour 
une faille NS, le mouvement du bloc Est est vers le Nord. II y a une 
composante normale superposee a la composante de decrochement. 

WEST CHILE RISE 

Mw 6.2 

(USGS Rapid Moment-Tensor Solution) 



Date: 7 SEP 1 998 
Time: 0:39:31 .63 UTC 
Epicenter: -35.902 -97.274 
Depth: 13 km 
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Figure 6.6 Un tremblement de Terre le long de la faille de la Reine 
Charlotte, a l’ouest du Canada. La faille est orientee NWN, et le 
mecanisme indique un mouvement NW du bloc Pacifique 

QUEEN CHARLOTTE ISLANDS REGION 

Mw 6.1 

(USGS Rapid Momeni-Tensor Solution) 


Date: 30 AUG 1998 
Time: 11:33:34.10 UTC 
Epicenter: 51 .131 -130.434 
Depth: 3 km 


De part et d’autre de la vallee centrale de la dorsale, il y a des failles 
normales. Les mecanismes au foyer des tremblements de Terre qui provi- 
ennent des flancs de la dorsale montrent bien un mouvement correspon- 
dant a une extension le long de ces failles normales. 

La sismicite des zones de subduction est beaucoup plus complexe. 
Elle implique toute la zone de Wadati-Benioff jusqu’a 600 km de pro- 
fondeur. Elle implique de la deformation dans les iles en arc derriere la 
zone de subduction. Pres de la surface, les mecanismes au foyer de la 
subduction sont des mecanismes de failles inverses. Plus en profondeur 
(100-300km), les mecanismes au foyer montrent qu’il y tension dans la 
plaque en subduction. Les mecanismes changent pour de la compression 
pour des profondeurs > 400km. 

Dans File en arc, en arriere de la zone de subduction, les mecanismes au 
foyer indiquent des failles normales et un regime d’extension. 

6.3.4 Sismicite intraplaque 

Bien que la plupart des tremblements de Terre se produisent aux frontieres 
de plaques, il existe des zones sismiques loin de toute frontiere entre 
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Figure 6.7 Un tremblement de Terre pres de la fosse des Nouvelles- 
Hebrides. Le mouvement correspond a un chevauchement. Le plan de 
faille pourrait etre incline soit vers le SW, soit vers le NE. Dans cette 
region oil la plaque Australienne passe sous la plaque Pacifique, une 
faille inclinee vers le NE est plus vraissemblable. 

VANUATU ISLANDS 

Mw 6.3 

(USGS Rapid Momeni-Tensor Solution) 



Date: 21 SEP 1998 
Time: 12:9:39.19 UTC 
Epicenter: -13.438 1 66.298 
Depth: 34 km 


plaque. Ainsi, es tremblements de terre les plus violents historique- 
ment rapportes en Amerique du Nord sont ceux de 1811 et 1812, a 
New Madrid, dans la vallee du Mississippi. La vallee du Mississippi est 
le centre d’une region d’activite sismique qui continue aujourd’hui. Au 
Canada, il y a de 1’activite sismique dans la vallee du Saint Laurent, 
principalement autour de l’astrobleme de La Malbaie. II y a egalement 
un peu d’activite sur la marge Atlantique. 

II est peu vraissemblable que ces crises de sismicite intraplaque soient 
importantes a 1’echelle geologique. Dans du Saint-Laurent, le debut de 
l’activite presente n’est pas connu, mais il est etabli qu’elle dure depuis 
au moins lOka. Il est probable que cette activite est sporadique et qu’elle 
pourrait se deplacer dans des regions aujourd’hui calmes. 
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Figure 6.8 La sismicite des Etats Unis. Notez l’absence de sismicite 
de la zone de subduction de Cascadia en Oregon et la sismicite in- 
traplaque de l’ouest des USA. 

Seismicity of United States: 1975 - lit* 



Figure 6.9 La sismicite du Canada. 

Seismicity of Canada: 1977 - 1997 JT ^ 



6.4 Notes supplementaires: Tenseur des contraintes. 

En general, la force, F qui agit a travers une surface dans un milieu con- 
tinu, depend de l’orientation de cette surface. L’orientation d’une surface 
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est definie par la direction de la normale a cette surface, c.a.d. par les 
composantes du vecteur unite n normal a la surface; les composantes de 
ce vecteur sont les cosinus directeurs de la normale (cos de Tangle entre 
la normale et chacun des axes d’un systeme cartesien de coordonnees). 
Si la force depend lineairement de Torientation de la normale, elle est 
donnee par: 

F x = T xx n x + T xy n y + T xz n z (6.3) 

Fy — Ty X Tl X “t“ TyyJly Ty Z Tl Z (6-4) 

F z = T zx n x + T zy n y + T zz n z (6.5) 

Par convention, la normale est prise vers Texterieur d’un volume. Les 
contraintes sont done definies positives vers Texterieur (signe contraire 
de la pression). La contrainte est une force par unite de surface, c.a.d. 
une pression dont l’unite de mesure est le N m~ 2 ou Pa. Le gradient de 
pression lithostatique dans la croute est de Tordre de 30 MPa. m -1 . 

En utilisant une notation avec indices (i,j = 1,2,3 pour x,y,z), la com- 
posante Fi du vecteur est donne par: 

Fi = TijUj ( 6 . 6 ) 

avec sommation sur Tindice qui est repete. Cette relation definit la force 
par unite de surface en fonction du tenseur des contraintes T. (On intro- 
duit le concept de tenseur parce qu’il s’agit d’une grandeur physique done 
independante du systeme de coordonnees. Par contre, les composantes de 
ce tenseur dependent du systeme de coordonnees). On peut montrer que 
cette definition conduit a une relation d’equilibre entre toutes les forces 
agissant sur la surface d’un cube ou d’un triedre, en considerant que la 
force agissant sur une surface perpendiculaire a l’axe des x a pour com- 
posantes T xx , T yxi T zx , la force agissant sur une surface perpendiculaire a 
l’axe des y a pour composantes T xy , T yy , T zy , et la force agissant sur une 
surface perpendiculaire a l’axe des z a pour composantes T xz ,T yz ,T zz . 

La condition d’equilibre du moment des forces agissant sur un cube 
implique que le tenseur des contraintes est aussi symetrique, cad: 

T xy = T yx ; T xz = T zx ; T yz =T zy (6.7) 

6.4.1 Contraintes principales 

Toute matrice symetrique d’ordre 3 admet 3 valeurs propres et vecteurs 
propres, c.a.d. 3 solutions a Tequation: 

Av = Av (6.8) 
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Pour chaque valeur propre A il existe un vecteur propre Vi tel que 
l’equation est satisfaite. On peut demontrer que: 

1. Si les valeurs propres sont distinctes, les vecteurs propres sont orthog- 
onaux. 

2. Si deux valeurs propres sont egales, toute combinaison lineaire des 
vecteurs propres correspondant est vecteur propre. 

3. Si les trois valeurs propres sont egales, tout vecteur est vecteur propre. 
Dans ce cas, la matrice A est multiple de la matrice unite et ses 
composantes sont les memes dans tous les sytemes de coordonnees. 

Dans le systeme de coordonnees propres (les axes sont definis par les 
directions de trois vecteurs propres orthogonaux), le tenseur des con- 
traintes est diagonal: 

/ <71 0 0 \ 

T = 0 < 7-2 0 (6.9) 

\ 0 0 <73 / 

Les valeurs propres du tenseur des contraintes (par convention <j\ > oy > cr 3 ) 
sont dites contraintes principales et les axes definissent les directions 
principales (directions selon lesquelles les contraintes sont paralleles aux 
normales). 

Dans le cas ou les trois contraintes principales sont egales, la contrainte 
est toujours normale a la surface et independante de la direction. Le 
champ de contrainte est dit lithostatique (la contrainte est egale a la 
pression due au poids des roches). 

Les contraintes sur un plan parallele a la direction de oy et faisant 
un angle 6 avec la direction de ay peut se decomposer en contraintes 
normale a n et tangentielle oy. La contrainte tangentielle correspond au 
cisaillement. Les contraintes normale et tangentielle sont egales a: 

a n = (Ji cos 1 2 6 + 03 sin 2 6 (6.10) 

a t = (<7i — (73) cos0sin0 = i(cri — er 3 ) sin(2d) (6-11) 

Le cisaillement est done maximum le long de plans a ± 45 °entre les 
directions des contraintes maximum et minimum. En principe, e’est done 
le long de ces plans que la fracturation se produit. 

Ceci permet de determiner le type de faille en fonction de l’orientation 
des contraintes principales. Trois cas sont a distinguer: 

1. ( 7 1 et <r 3 horizontales: faille de decrochement. 

2. (j 1 horizontale et cr 3 verticale: regime tensionnel: faille normale. 
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3. a i verticale et <73 horizontale: regime compressif: faille inverse. 

La convention rle signe est que la contrainte est positive vers l’exterieur. 
En general, la pression lithostatique est dirigee vers l’interieur et les 
contraintes sont negatives, la contrainte principale maximum est done 
la plus petite en valeur absolue. 


6.4.2 Effet de la friction 

En pratique les failles normales et inverses font rarement un angle de 45 
°avec la verticale. Ceci est du au fait que la friction le long de la faille 
n’a pas ete consideree. La loi d’Amonton suppose que la friction est pro- 
portionnelle a la pression lithostatique. La constante de proportionalite 
est le coefficient de friction. La loi de Byerlee est basee sur des mesures 
en laboratoire des variations de la contrainte de cisaillement necessaire 
pour surmonter la friction et declencher la rupture en fonction de la 
pression effective. 

a s =a(a n -p w ) (6-12) 


oil le coefficient de friction a est 0.8 pour une contrainte effective a n — p w < 
200 MPa et 0.6 pour a n — p w > 200MPa. La contrainte effective a n — p w 
est la difference entre la pression lithostatique et la pression des fluides 
dans les pores de la roche, souvent supposee egale a la pression hydro- 
statique. 

En tenant compte de la friction, la contrainte est maximum pour un 
angle 0 tel que: 

tan(29) = ±— (6.13) 

a 


oil le signe est positif pour une faille normale et negatif pour une faille 
de chevauchement. La difference entre la plus grande et la plus petite 
contrainte principale, <x 1 — 03, doit etre au moins egale a: 


cri - a 3 > 


2 0.(7 n 

'Jo? + 1 ± a 


(6.14) 


Le signe au denominateur depend de la direction de la contrainte princi- 
pale maximum. II est positif si la contrainte maximum est horizontale et 
negatif si elle est verticale. En pratique, cela implique que la difference 
des contraintes principales doit etre beaucoup plus grande pour sur- 
monter la friction pour une faille de chevauchement que pour une faille 
normale. (En supposant que la pression de fluide est hydrostatique, la 
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difference est de l’ordre de 2 fois la contrainte lithostatique en chevauche- 
rnent et de 0.7 fois la contrainte lithostatique pour une faille normale). 
La lithosphere est done plus resistante a des contraintes compressives 
que tensiles. 

6.4.3 Changement de systeme d’axes de coordonnees 

Dans un systeme cartesien de coordonnees, un vecteur v est exprime par 
ses composantes v-i en fonction des vecteurs unites ei dans la direction 
des axes. 

v = i^ei (6.15) 

Pour un systeme cartesien, o; ■ ej = Sij ( 5ij = 1; i = j;5ij = 0 ;i ^ j). 
Dans un autre systeme cartesien dont les vecteurs unites sont el, le 
vecteur v est donne par: 

v = v^i = u'el = v'j (el ■ ei)e ; = v^R^e ; (6.16) 

La matrice R est une matrice de rotation qui permet de transformer les 
composantes d’un vecteur d’un systeme cartesien a un autre. Elle est 
orthogonale: 

RijRkj = $ik (6-17) 

R t R = RR t = I (6.18) 

OU r t est la transposee de R et I denote la matrice unite. Pour un 
tenseur cartesien de composantes Tjj, le changement de systeme de co- 
ordonnees est donne par: 

Tij = Rki.RijT kl (6.19) 

Apres changement du systeme de coordonnees, l’equation aux valeurs 
propres, Av = Av, peut s’ecrire: 

A'v' = A'Rv = Av' = ARv (6.20) 

et done: 

A = R t A'R (6.21) 

A' = RAR t (6.22) 
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7.1 Introduction 

C’est un nouveau type de magnetometre congu par le physicien Black- 
ett qui a perrnis le developpement du paleomagnetisme dans les annees 
1950. Blackett, qui voulait expliquer le champ magnetique de la Terre 
et du soleil, avait propose que tout corps en rotation produit un champ 
magnetique. II avait mis au point un magnetometre pour tester sa theorie 
qui a rapidement ete abandonnnee. Mais ce nouveau magnetometre al- 
lait etre utilise par Keith Runcorn et ses etudiants, notament Ted Irv- 
ing, qui allaient rnettre au point la methode paleomagnetique. II est 
surprenant que le paleomagnetisme n’ait pas joue un role plus impor- 
tant dans l’histoire de la tectonique des plaques. Les premieres etudes 
(1955) qui demontraient le mouvement des continents relativement au 
pole magnetique ont ete ignorees par l’ensemble de la communaute 
geoscientifique. Ce n’est qu’apres la demonstration de l’expansion des 
fonds oceaniques que les etudes paleomagnetiques ont ete utilisees de 
fagon tout a fait systematique pour reconstruire les positions des conti- 
nents dans le passe. 

Contrairement aux anomalies magnetiques marines qui ne couvrent 
que 200 Ma, les donnees paleomagnetiques qui proviennent des conti- 
nents couvrent la plus grande partie de l’histoire de la Terre et permet- 
tent done de retracer des evenements qui remontent jusqu’a l’Archeen. 


7.2 Champ magnetique terrestre 

Le champ magnetique terrestre varie en grandeur et direction a la surface 
de la Terre. Localement, on peut mesurer la direction de la composante 
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Figure 7.1 Les elements du champ magnetique. 



horizontale du champ magnetique (ou declinaison) et l’angle du champ 
magnetique avec l’horizontale (ou inclinaison) . 

En premiere approximation, le champ magnetique terrestre est celui 
d’un dipole, dont l’axe est maintenant incline de 11 °par rapport a l’axe 
de rotation. 
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Figure 7.2 Le champ magnetique terrestre est celui d’un dipole. 



Pour un dipole, le champ magnetique a la surface de la terre varie du 
pole nord au pole sud en fonction de la colatitude 9: 

H r = cos 6 (7.1) 

a 6 

M 

Hg = — sind (7.2) 

a 6 

H r est la composante radiale (verticale) et Hg est la composante tan- 
gentielle (horizontale), M est le moment dipolaire, et a est le rayon 
terrestre. Pour un dipole, le champ ne depend done pas de la longitude, 
mais uniquement de la latitude. L’inclinaison et la declinaison dependent 
de la position par rapport au pole. Le rapport de la composante verti- 
cale a la composante horizontale du champ magnetique est determine 
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par la colatitude (distance au pole) et en mesurant leur rapport (qui 
est la tangente de l’inclinaison, tan/), il est possible de determiner la 
colatitude, 9 , ou la latitude, A, a l’aide des formules: 


. i , tan / 

0 = cotan ( — - — ) 

(7.3) 

, i , tan / 

A = tan ( 2 ) 

(7.4) 


L’inclinaison est positive vers le bas: la latitude est done bien definie 
positive dans /hemisphere nord. 

En fait, le champ magnetique terrestre n’est pas exactement celui d’un 
dipole. Le champ magnetique peut etre decompose en une composante 
dipolaire et des composantes multipolaires (qui varient sur de plus cour- 
tes distances). II faut aussi distinguer le pole geomagnetique, le point 
ou l’axe du dipole traverse la surface de la terre, du pole magnetique, 
le point ou le champ magnetique est vertical. II n’y a qu’un seul pole 
geomagnetique; par contre il y a deux poles nord magnetiques: l’un au 
Canada, l’autre en Siberie. Aujourd’hui, le pole geomagnetique est in- 
cline de 11 “par rapport au pole nord geographique. 

En principe, cad pour autant que le champ soit dipolaire, /equation 7.3 
permet de calculer la distance au pole geomagnetique. 

Outre les inversions de polarite, le champ magnetique varie a courte 
echellc de temps (100-1000 ans). Ces variations affectent surtout les com- 
posantes non dipolaires qui semblent se deplacer vers l’ouest ( westward 
drift). La vitesse apparente de ce drift est 0.2 deg /a. Par ailleurs, cer- 
taines mesures sismiques avaient suggere que la vitesse de rotation du 
noyau terrestre est differente de celle du manteau [Song & Richards , 
1993]. Cette rotation differentielle pourrait etre la cause du westward 
drift. Ce resultat a ete cependant remis en question par une etude 
recente des oscillations propres; cette etude a conclu que le taux de 
rotation differentielle etait faible ou nul [Laske & Masters, 1999]. Par 
ailleurs, /amplitude du champ dipolaire decroit depuis plusieurs cen- 
taines d’annees et l’axe du dipole semble tourner autour du pole nord 
geographique. Cette derniere observation suggere qu’en moyenne l’axe 
du dipole pourrait coincider avec l’axe de rotation de la terre. 


7.3 Proprieties magnetiques des roches 

Les proprietes magnetiques des roches ont ete decrites dans le chapitre 3. 
La magnetisation thermoremanente est acquise lorsque les mineraux se 
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Figure 7.3 Les variations recentes de la position du pole magnetique 
an nord du Canada 
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refroidissent et que leur temperature devient inferieure a la temperature 
de Curie ( Tc ~ 400 — 700 °C). La magnetisation ne change plus quand 
la temperature est inferieure a la temperature de blocage, (~ 50 °sous 

Tc)- 

Ce sont surtout les roches mafiques, riches en mineraux magnetiques, 
qui fournissent les meilleurs echantillons pour les etudes paleomagnetiques. 

En principe, les meilleurs echantillons proviennent de roches qui se sont 
refroidies rapidement et qui enregistrent done une direction bien definie 
du champ magnetique terrestre. Les meilleurs echantillons pour le paleomagnetisme 
proviennent done de dykes (parce qu’ils se refroidissent rapidement) de 
composition mafique, tels des dykes de diabase. Le temps de refroidisse- 
ment depend de la taille (quelques semaines pour un dyke epais d’un 
metre, quelques millions d’annees pour un pluton a quelques kilometres 
de profondeur. Plusieurs echantillons provenant de tels plutons permet- 
tent d’obtenir une meilleur moyenne (representant mieux le dipole moyen 
dont l’axe coincide avec l’axe de rotation). Par contre, on y trouve sou- 
vent des polarites differentes du champ magnetique terrestre. Ces inver- 
sions de polarite sont faciles a detecter et ne posent pas de problems 
majeur pour ^interpretation. 


7.4 Determination des paleopoles 

L’objectif principal d’une etude paleomagnetique est la determination 
de l’age d’un echantillon et du paleopole dont la position est determines 
en supposant que l’echantillon est reste fixe. Cette maniere de proceder 
est necessaire car il n’est pas possible de determiner la longitude et done 
une position geographique unique de l’echantillon par rapport a un pole 
fixe. 

La methode paleomagnetique est basee sur plusieurs hypotheses. 

1. La direction de l’axe du dipole qui represente le mieux le champ 
magnetique est restee en moyenne parallele a l’axe de rotation de 
la Terre. Ceci n’exclut pas les inversions de la polarite du champ 
magnetique sans changement de l’axe du dipole. Pour des echantillons 
recents (< 20Ma), on a pu verifier que la moyenne de tous les paleopoles 
coincidait avec le pole geographique. Pour des echantillons plus vieux, 
l’hypothese est incompatible avec l’immobilite des continents. 

2. Une moyenne temporelle (sur plusieurs ka) du champ magnetique 
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terrestre ne contient que la composante dipolaire et les composantes 
multipolaires s’annulent. 

3. Les roches magnetisees ont conserve la direction du champ magnetique 
terrestre au moment ou elles ont passe la temperature de bloquage. 
II existe d’autres mecanismes de magnetisation que la magnetisation 
thermoremanente (magnetisation visqueuse, alteration chimique, etc.). 
Ces magnetisations secondaires peuvent etre eliminees en laboratoire 
par demagnetisation progressive de l’echantillon. Cela se fait soit en 
chauffant l’echantillon par paliers jusqu’a la temperature de Curie, 
soit en plagant l’echantillon dans un champ magnetique alternatif 
dont l’intensite augmente par paliers jusqua demagnetisation complete. 
Pour chaque palier, on mesure la direction de la magnetisation. En 
principe, la magnetisation primitive sera la plus resiliente et done la 
derniere a disparaitre. 

4. Les echantillons n’ont subi aucune rotation tectonique apres avoir 
acquis leur magnetisation. Si une rotation peut etre mesuree sur le 
terrain, une correction sera faite. 

5. II faut aussi noter que si la temperature de fermeture des mineraux 
utilises pour dater 1’echantillon est different de la temperature de 
Curie, la geochronometrie ne permet pas de determiner Page de la 
magnetisation. Pour les dykes qui se refroidissent rapidement, ces 
deux ages seront proches, pour autant qu’il n’y ait pas eu d’evenement 
thermique apres la mise en place du dyke. 

En pratique, pour determiner un paleopole, il faut collecter plusieurs 
echantillons orientes de la merne roche. Plusieurs echantillons fournissent 
une moyenne temporelle qui elimine les composantes multipolaires. Si 
une rotation tectonique est evidente sur le terrain, on peut en tenir 
compte et faire une correction pour chaque echantillon. En laboratoire, 
le magnetometre permet de mesurer sur ces echantillons orientes la 
declinaison (qui donne la direction d’un arc de grand cercle passant 
par le pole) et Pinclinaison (qui donne la distance au pole le long d’un 
arc de grand cercle). A partir de ces mesures, il est facile de determiner 
la position du paleopole par rapport a l’echantillon. Ceci peut se faire 
par le calcul (voir appendice) ou graphiquement a l’aide des projections 
stereographiques [Cox & Hart , 1986, p315-316]. Pour chaque echantillon, 
on determine un pole geomagnetique virtuel. Le paleopole est obtenu 
comme moyenne des poles virtuels. Les differences mineures entre les 
poles virtuels sont dues surtout a des differences dans le temps de pas- 
sage de la temperature de bloquage et de la direction du dipole au mo- 
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ment du passage par la temperature de bloquage. Des differences irnpor- 
tantes indiquent souvent que la roche a ete deformee apres avoir acquis 
sa magnetisation. 

II faut noter que Ton peut determiner le paleopole de faqon unique, 
rnais qu’il est impossible de determiner uniquement la position de l’echantillon 
par rapport a un pole fixe. II est impossible de determiner la longitude. 

C’est une des raisons pour laquelle la methode paleomagnetique n’avait 
pas pu demontrer de maniere convaincante l’ouverture de l’Atlantique. 


7.5 Applications du paleomagnetisme 

A partir des nombreuses mesures paleomagnetiques, le mouvement du 
paleopole en fonction de Page a pu etre calcule pour les differents con- 
tinents depuis LArcheen jusqu’au present. Les resultats globaux ont ete 
compiles par Gordon (1995) et plus recemment par Torsvik et al. (2008). 

Une etude paleomagnetique consiste soit a comparer le paleopole d’un 
echantillon date avec celui du continent, soit a comparer les paleopoles 
de differents echantillons de meme age. Le paleomagnetisme a permis de 
demontrer entre autres: 

1. Les variations de la distance relative entre sites ou des echantillons de 
meme age ont ete recueillis. Ainsi, le mouvement recent de Linde par 
rapport a l’Asie est facilement demontre. Des echantillons de meme 
age provenant de Linde et de l’Asie ont des paleopoles differents. 
La distance entre ces paleopoles correspond a la distance minimum 
parcourue par Linde relativement a l’Asie. En supposant l’Asie fixe, 
Linde peut etre replacee a sa paleolatitude par rapport au paleopole 
de l’Asie. (A cause de Lambiguite sur la paleolongitude, il y a plusieurs 
positions possibles du continent indien). 

2. La rotation tectonique de fragments de continents. (Rotation de LItalie 
ou de la peninsule iberique par rapport a l’Europe). Les declinaisons 
d’echantillons de meme age provenant de la peninsule iberique ou des 
Appenins sont differentes de celles du reste de LEurope pour un meme 
age. La figure 7.5 illustre un exemple de rotation. 

3. Le paleomagnetisme fournit une contrainte essentielle pour recon- 
struire les paleocontinents. Par exemple, les reconstructions des con- 
tinents autour de l’ocean Atlantique avant l’ouverture d’lapetus, ou 
les reconstructions du supercontinent Rodinia doivent etre compati- 
bles avec toutes les donnees paleomagnetiques. 
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Figure 7.4 Paleopoles de l’Amerique du Nord pendant l’orogenie 
Grenvillienne entre 1100 et 800 Ma. Ces paleopoles impliquent un 
mouvement du continent de plus de 90°ou 10,000km. Dans ce cas, 
le paleomagnetisme pourrait servir a dater des echantillons avec une 
precision de l’ordre de 25Ma. 
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Figure 7.5 Paleopoles d’echantillons contemporains provenant du 
nord et du sud de la region des Afars. Ces paleopoles ont la meme 
inclinaison mais les differentes declinaisons impliquent une rotation 
du bloc des Afars pax rapport a l’Afrique. 
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Figure 7.6 Paleopoles d’echantillons de 520 Ma provenant de la 
Cordillere sud americaine en Argentine. Si l’Amerique du Sud est 
attachee a l’Afrique, les paleopoles sont incompatibles avec tous les 
autres poles du meme age. Pour reconcilier tous ces paleopoles, il 
faudrait que les terrains sud americains aient ete proches de Lauren- 
tia (continent nord-americain). 




4. L’origine allochtone de certains terrains (Terrains suspects dans la 
Cordillere Canadienne, les Appalaches, ou en Italie). 

5. La separation passee entre des blocs de croute continentale aujourd’hui 
soudes (Orogene Trans Hudsonien) et l’existence probable de micro- 
continents a l’Archeen et au PaleoProterozoique. 
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6. Une autre application elegante a ete la demonstration de deformations 
intracontinentales. Certains essaims de dykes, qui sont longs de plusieurs 
milliers de km, traversent des regions qui ont ete deformees apres leur 
emplacement. Ainsi, on a pu demontrer que le soulevement de Ka- 
puskasing a l’ouest de l’Abitibi, est posterieure a l’intrusion des dykes 
de Matachewan et calculer la deformation tectonique [Halls et al., 

1994]. Les cartes d’anomalies magnetiques demontrent assez claire- 
ment la deformation. II est possible de ’’reconstruire” la geometrie 
des differentes ceintures geologiques avant deformation en redressant 
les dykes. Dans ce cas, la reconstruction doit satisfaire les donnees 
paleomagnetiques (les paleopoles des differents dykes ainsi reconstru- 
its doivent coincider). 

7. En general, la methode paleomagnetique consiste a determiner le 
paleopole d’un echantillon dont on a mesure l’age independament. 

Elle peut aussi etre utile pour dater des echantillons dont on a determine 
le paleopole et en supposant que celui-ci coincide avec le paleopole 
regional. La precision d’une telle methode est mediocre et la methode 
n’est possible que pour les periodes pendant lesquelles le continent 
subissait des deplacements NS rapides. Toutefois, cette methode per- 
met de dater a peu de frais des roches (surtout recentes) ou de fournir 
un age approximatif quand les methodes usuelles de chronometrie 
sont inadequates. 

Bien que d’autres donnees (paleontologiques, paleoclimatiques, etc.) 
fournissent des indications sur la paleoposition des continents, le paleomagnetisme 
fournit la seule donnee qui permette de determiner avec precision la lat- 
itude passee d’un continent. II est maintenant possible de reconstruire 
revolution des differents blocs continentaux depuis au moins la fin de 
l’Archeen. 


7.6 Appendice: Determination du paleopole 

Les mesures paleomagnetiques donnent l’inclinaison I et la declinaison 
D de la magnetisation d’un echantillon. La colatitude 9 est determinee 
par l’inclinaison / 


tan / = 2cotan0 


(7.5) 
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Si A s et <fis sont les coordonnees du site les coordonnees du pole X p et 
(f> p sont: 

sin X p = sin A s cos 9 + cos A s sin 9 cos D 

<t>p = <t>s + /3 si cos 9 > sin A s sin X p 
<f>p = <j> s + 180 — p si cos 9 < sin A s sin X p 
sin j3 = sin 9 sin D / cos X p 
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Applications de la geodesie spatiale a la 
tectonique des plaques 


Le developpement de la tectonique des plaques a ete possible parce que 
les anomalies magnetiques marines ont perrnis de mesurer la formation 
de plancher oceanique a l’echelle de temps geologique. Le developpement 
de la geodesie spatiale au cours des 20 dernieres annees a perrnis de 
mesurer les deformations tectoniques en ’’temps reel”. De plus, les tech- 
niques de geodesie spatiales ont perrnis de quantifier toutes les deformations 
tectoniques. II est maintenant possible de mesurer directement des taux 
de convergence et les deformations a l’interieur des plaques. Cette revolution 
technologique a done eu de tres importantes implications conceptuelles. 

Grace aux satellites, il est maintenant possible d’obtenir rapidement 
et sur l’ensemble du globe des donnees geophysiques et geodesiques avec 
une resolution et une precision jamais atteintes auparavant. Les donnees 
sismiques avaient perrnis de demontrer la direction des mouvements et 
contraintes aux frontieres de plaque. Les methodes de geodesie spatiale 
permettent de quantifier ces mouvements. II est possible maintenant de 
mesurer en temps reel les deformations tectoniques et la vitesse relative 
des plaques. Ces techniques se sont developpees suite au deployment de 
satellites de positionement et d’observation, grace a la technologie spa- 
tiale qui permet de recueillir rapidement des donnees precises et surtout 
grace aux ordinateurs qui ont perrnis de les traiter. Le livre de Cazenave 
and Feigl (1994) contient une bonne revue des techniques nouvelles de 
geodesie spatiale. 

Les determinations de distances demandent des mesures precises du 
temps de trajet d’ondes electromagnetiques (visible, IR, radar, etc). 

II est souvent suffisant de mesurer des differences de temps de par- 
cours. Pour avoir une precision de l’ordre de 3cm, il faut mesurer des 
differences de temps avec une precision < 10” lo s. Ceci ne peut se faire 
par des methodes conventionnelles mais est possible avec les techniques 
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d’interferometrie. Pour mesurer les vitesses et accelerations de satellites, 
on utilise l’effet Doppler qui perrnet de mesurer de fagon quasi continue 
la vitesse dans la ligne de visee. 

Pour donner une idee des progres realises au cours des 20 dernieres 
annees, la precision des mesures a ete amelioree d’un facteur de 100. 
Ainsi, la position des satellites est maintenant determinee avec une 
precision de 1cm. Leur vitesse est connue avec une precision de l’ordre 
de 0.1 mrn/s. Le geoide est mesure avec une precision de 1cm et une 
resolution spatiale d’une dizaine de km. II faut noter que pour toutes 
ces mesures, les corrections sont souvent superieures de deux ou trois 
ordres a la precision des mesures. Ainsi, il est ncessaire de tenir compte 
du parcours des ondes dans Patmosphere et des variations de cette at- 
mosphere. Ces corrections pour le positionement sont de l’ordre de lm. 
De meme, pour mesurer l’acceleration de la gravite il faut tenir compte 
du freinage atmospherique et du freinage par le rayonnement solaire. La 
precision des mesures corrigees perrnet done de determiner des variations 
centimetriques de la distance entre stations au sol, ou des variations cen- 
timetriques du geoide causees par les courants marins ou des phenomenes 
tel El Nino. 

Deux satellites mis sur la sur la meme orbite permettent de mesurer 
de fagon continue les variations dans le temps du champ de gravite. Il 
s’agit de la mission GRACE Gravity Recovery and Climate Experiment) 
qui a debute en Mars 2002. Ces satellites permettent de mesurer les vari- 
ations du champ de gravite terrestre avec un precision de ss lz/Gal alors 
que les gravimetres conventionnels au sol ont rarement une precision 
< 10/rGal II deviendra possible de mesurer les changements de la gravite 
produits par un tremblement de Terre, par la fonte de glaciers, les 
deplacements de la nappe phreatique, ou par des changements dans les 
courants oceaniques superficiels ou profonds. 


8.1 Telemetrie satellite. SLR (Satellite Laser 

Ranging) 

A l’aide d’ ondes laser reflechies par les satellites, les stations de pour- 
suite au sol doivent determiner la position exacte de satellites. Au cours 
du temps, il a ete possible de montrer les changements de distance entre 
stations au sol dus au mouvement des plaques. 

La vitesse est mesuree instantanement par effet Doppler (change- 
ment de frequence de l’onde reflechie). Cela perrnet de determiner aussi 
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l’acceleration et done de determiner le champ de gravite global. Les 
mesures absolues de gravite au sol etaient peu nombreuses (les gravimetres 
ne permettant de mesurer que les variations de g ). Les mesures a partir 
de satellites ont permis de determiner les valeurs absolues de g avec une 
resolution d’une centaine de km. 


8.2 GPS-Global Positioning System 

Les techniques GPS utilisent les signaux emis par un ensemble d’une 
trentaine de satellites de navigation. Les horloges des differents satellites 
et celle du recepteur GPS ne sont pas synchronisees avec suffisament de 
precision. II faut done au moins quatre satellites pour pouvoir calculer 
les differences entre leurs horloges ainsi que la distance au recepteur. 
Chaque satellite emet deux signaux avec des longueurs d’onde differentes 
dans les bandes centimetriques («25cm). Le signal est egalement code 
suivant une sequence qui donne l’heure de l’horloge interne du satellite. 
La difference de phase entre les deux signaux est proportionnelle a la 
distance parcourue. Le codage du signal peut etre modifie pour que 
les mesures precises ne soient pas possibles pour les utilisateurs non 
autorises. Pour ameliorer la precision localement on peut utiliser un 
systeme differentiel avec deux GPS. 


Des mesures GPS repetees et de grande precision permettent de demontrer 
les taux de deformations dans les cliaines de montagnes actives (Alpes 
et Himalayas). Elies ont aussi mis en evidence des taux importants de 
deformation a l’interieur de l’Asie Argus and Gordon (1996). Les mesures 
GPS en Oregon, a proximite de la zone de la zone de subduction de Cas- 
cadia, ont demontre la rotation de l’Oregon par rapport a l’Amerique du 
Nord [McAffrey et al., 2000]. Pour l’Oregon, les variations du taux de 
deformation sont en partie dues au verrou de la zone de subduction et 
en partie a l’effet de l’extension dans la Province du Basin and Range. 
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Figure 8.1 Les deformations le long de la zone de subduction de Cas- 
cadia en Oregon mesurees par GPS. 



8.3 Interferometrie a large base: VLBI (Very Long 
Baseline Interferometry) 

Pour mesurer les positions relatives d’un certain nombre de stations 
fixes, on utilise les signaux emis par les quasars. 1 Ces signaux tres car- 

1 La technique a ete developpee par les astrophysiciens dans le but d’etudier les 
objets quasi-st ellaires. Par la suite, les geodesistes ont pu profiter des 
observatoires et des mesures astrophysiques. 
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Figure 8.2 Les deformations de la Californie mesurees par VLBI et 
GPS. 


Crustal Deformation Velocity wrt North American Plate 



acteristiques sont enregistres en des sites differents. En correlant ces 
signaux, on peut determiner la difference de temps de parcours avec 
une precision de l’ordre de la ps (10~ 12 s) et done la projection de la 
distance entre les stations sur la direction du quasar avec une precision 
de l’ordre du mm. En repetant les mesures avec plusieurs quasars, on 
determine les positions relatives de toutes les stations fixes. La precision 
des mesures repetees sur plusieurs annees est suffisament grande pour 
que l’on puisse maintenant determiner les mouvements des plaques tec- 
toniques. Le mouvements des plaques calcules par VLBI coincident bien 
avec les mouvements calcules par les methodes conventionnelles (anoma- 
lies magnetiques marines, etc) Gordon and Stein (1992). 

Un des resultats les plus interessants a ete la demonstration que la zone 
de deformation entre plaques est souvent diffuse. Ainsi les deformations 
entre la plaque americaine et pacifique en Californie impliquent non 
seulement la faille de San Andreas mais toute la Sierra Nevada. 
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8.4 Altimetrie satellite 

L’altimetrie satellite a ete utilisee pour mesurer precisement la surface 
du geo’ide en mesurant le temps d’aller retour d’un signal radar ernis par 
le satellite et reflechi par la surface de la mer. Ces mesures ont permis 
de determiner le champ de gravite sur tous les oceans 2 et souvent d’en 
deduire la bathymetrie (Smith and Sandwell, 1994). 

S’il y avait parfaite compensation isostatique, les anomalies d’air li- 
bre seraient nulles et on ne verrait pas la bathymetrie sur les cartes 
gravimetriques. C’est parce quc les courtes longueurs d’onde de la bathymetrie 
(<200 km) ne sont pas compensees que le geo'ide reflete la bathymetrie. 

En fait, les anomalies de geo'ide ne servent qu’a interpoler entre les 
mesures de bathymetrie faites par les batiments oceanographiques. 

La carte des fonds marins qui en resulte a une resolution qu’il serait 
impossible d’obtenir avec des mesures traditionelles. Par exemple, on 
voit tres clairement les zones de subduction sur la carte des fonds marins 
de la zone de Tonga Kermadec. En fait, cette carte est presque suff- 
isante en elle meme pour tracer toute l’activite tectonique des plaques 
oceaniques. 

Des mesures repetees ont permis de suivre les variations temporelles du 
geo'ide qui peuvent etre dues aux courants marins ou a des phenomenes 
environnementaux (variations de temperature, vents). Ainsi, les satel- 
lites tel TOPEX/POSEIDON ont maintenant pour mission de suivre le 
developpement des phenomenes du type El Nino. 


8.5 Interferometrie radar. (SAR- Synthetic Aperture 

Radar) 

On peut construire des interferogrammes entre signaux radar enregistres 
par les satellites a des temps differents. Apres une correction geometrique 
pour replacer les satellites au meme point, les interferences sont dues au 
deplacement de la surface du sol entre deux passages du satellite. On a 
pu ainsi mesurer avec grande precision les mouvements du sol apres des 
tremblements de terre, le mouvement de glaciers, ou les deformations 
de volcans. Ceci permet une mesure quasi instantannee de deplacements 
verticaux qui sont quasi impossibles a detecter meme avec des releves 
geodesiques precis. La figure 3 montre un exemple recent en Californie 

2 A ne pas confondre avec les mesures de g par poursuite de l’orbite des satellites 
qui ont une plus faible resolution 
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(Sandwell et al., 2000). La precision des mesures permet non seulement 
de mesurer la deformation co-seismique, mais aussi les reajustements 
post-seismiques [ Pollitz et al., 2001]. 

Ces techniques d’interferometrie ont une application immediate pour 
Evaluation des risques geologiques. Par exemple, l’interferometrie per- 
met de detecter les variations de la surface du sol causees par les re- 
montees de magma dans les chambres magmatiques [Lu et al., 2000] 
et eventuellement de prevoir les eruptions volcaniques. Jusqu’a present, 
cette methode n’a pas ete utile pour detecter les mouvements preseismiques 
de la surface du sol et prevoir les tremblements de Terre. 

Des interferences entre images radar prises avec des angles differents 
peuvent etre utilisees pour calculer de fagon precise la topographie. On 
a ainsi etabli un modele digital de la topographie a l’echelle globale. 


exosup.com 


page facebook 


102 Applications de la geodesie spatiale a la tectonique des plaques 


Figure 8.3 Les franges d’interferences causees par la deformation 
apres le seisme de la mine Hector en Californie (M=7.1). Les valeurs 
annotees sur les lignes de contour correspondent au deplacement en 
mm. 



8.6 Application en planetologie 

Des techniques similaires ont permis d’imager la surface des planetes et 
de mesurer le champ de gravite en suivant la trajectoire des satellites. 
Pour Venus dont la surface est cacliee par une epaisse atmosphere, les 
premieres mesures radar faites a partir d’observatoires terrestres man- 
quaient de resolution. Des satellites d’observation ont permis des mesures 
tres precises dans l’infrarouge. La surface de MARS a ete mesuree avec 
une precision et une resolution remarquables lors de la mission MOLA 
(Mars Orbiter Laser Altimetry) conduite par NASA [Smith et al . , 1999]. 
L’etude de l’orbite des satellites a permis de determiner les champs de 
gravite de la Lune, ainsi que de Mars et Venus. Aujourd’hui, la resolution 
est bonne pour la Lune et Mars, mais plus faible pour Venus. 
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Outre la prediction des risques naturels, les mesures satellites servent a 
detecter les changements environnementaux naturels ou anthropiques. 
Les mesures d’interferometrie permette de determiner les mouvements 
des glaciers. Les mesures d’interferometrie ont permis de demontrer la 
subsidence due a l’extraction du petrole dans les puits du Texas ou de 
Californie. Les mesures d’altimetrie satellite permettent de determiner 
les variations des courants marins et les phenomenes transitoires de type 
El-Nino. 3 Les possibility offertes par le systeme GRACE sont encore 
plus etonnante, on pourra detecter l’effet sur le champ de gravite global 
de variations de la masse de glace au Greenland ou en Antartique et leur 
redistribution dans l’eau des oceans. 


L’effet Doppler est un changement de la frequence d’une onde emise par 
une source un mouvement par rapport a un observateur. La frequence 
to = 2 k/T ou T est la periode du signal, i.e. pour un signal sinusoidal, 
la temps entre 2 pics successifs. Si la source se deplace avec une vitesse 
v vers l’observateur, le temps entre deux pics successifs sera reduit de 
vT/c (c est la vitesse de l’onde). La periode pour l’observateur sera done 
T x ( 1 — u/c) et la frequence deviendra w/ ( 1 — v/c), cad sera plus elevee si 
v > 0 et plus faible si v < 0, cad la source s’eloigne. Le changement vers 
l’aigu du klaxon d’un vehicule qui s’approche est une exemple. C’est ce 
meme principe qui est utilise par les radars qui mesurent la vitesse des 
voitures. Le deplacement vers le rouge (basses frequences) des signaux 
de galaxies lointaines montrent qu’elles s’eloignent et que l’univers est 
en expansion. 


La superposition de deux ondes produit des phenomenes d’interferences, 
partout ou les deux ondes sont en phase, leur superposition produit un 
signal renforce, quand elles sont en opposition de phase, cad decalees 
dans le temps d’une derni periode elles s’annulent. L’interferometrie peut 

3 A ce sujet, il y a d’autres informations utiles, notament la tempearture de 
surface de la mer, qui sont aussi mesurees par les satellites. 
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etre utilisee, soit pour mesurer la difference de distance a la source en- 
tre deux instruments qui mesurent le meme signal (c’est le principe de 
l’interferometrie a tres longue base), soit pour mesurer la difference de la 
distance aller-retour de deux signaux emis par la meme et reflechis (c’est 
le principe de dl’interferometrie satellite). En principe, la precision de la 
rnesure de la distance est une fraction de la longueur d’onde. Ainsi pour 
les signaux radar, les frequences {v) sont comprises entre 200 MHz et 
20000 MHz ce qui implique des longueurs d’ondes (A) de 1.5m a 1.5cm 
(A = v/c, avec c = 3 x 10 8 m s -1 ). Les radiotelescopes utilisent a peu 
pres les memes gamrnes de longueur d’onde. 

Plus d’information au site du Jet Propulsion Laboratory 
http : / / Southport . j pi . nasa. gov / index, html 

http:/ /southport. jpl.nasa.gov/scienceapps/dixon/report2. html 

8.8.3 Mesure de la distance par GPS 

Pour mesurer la distance a un satellite, le recepteur GPS mesure le 
dephasage de deux signaux de frequences differentes. 
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Figure 8.4 La diffrence de phase entre deux ondes de frequences 
differentes permet de mesurer la distance a la source. 



8.8.4 Satellite Laser Ranging 

Laser ranging is the most accurate technique available for observing the 
orbits of artificial satellites. Although it can only be used with satellites 
carrying retro-reflectors, these devices are entirely passive, long-lasting, 
light and small, and so are relatively easy to fit to satellites whose precise 
positions need to be known. Sometimes the results may be needed for 
one specific purpose, but most satellite ranging observations can be used 
in many different, but inter-related, investigations. 

The technique uses time-of-flight measurements with short pulses of 
light to determine the range to a satellite as it moves in orbit around 
the Earth. The light pulses come from a pulsed laser and are beamed 
at the satellite. The retro-reflectors work like ‘cat’s eyes’ and reflect the 
light they receive back to the receiving telescope. From the pulse of 10 1 ' 
photons emitted towards the satellite only about one is recorded by the 
receiver in every 5 to 10 shots. 

The principle for determining distance is the same as that of radar 
except that light, rather than radio waves, is used. Because a narrow 
pass-band filter is used to cut out the light from the sky, the equip- 
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ment at Royal Greenwich Observatory can work both by day and night, 
provided that the skies are clear. Its single-shot precision of about 2 
cm is typical of the ten or so “third generation” systems that are now 
participating in a worldwide campaign of coordinated observations. A 
single satellite pass, lasting up to 40 minutes, may yield more than 1000 
returns. Each return corresponds to the detection of only one photon. 
Ranging accuracy is ultimately limited by the correction for the slower 
speed of light through the atmosphere. The correction amounts to about 
2m at the zenith (directly overhead) and is computed from meteorologi- 
cal data; the estimated uncertainty of the model is about 1 cm. Checking 
and setting up of the system is done by ranging to local targets, and the 
observations are dated to within 1 millionth of a second so that they 
may be properly combined in analyzes with others obtained elsewhere. 


8.8.5 La mission des satellites GRACE 

The Gravity Recovery and Climate Experiment (GRACE) mission has 
been competitively selected by NASA for development under a new Of- 
fice of Mission to Planet Earth program called Earth System Science 
Pathfinders (ESSP). The two satellites fly at an altitude of between 
300 and 450 kilometers in nearly the same orbital plane, that is, one 
’’chasing” the other, with a separation of about 200 kilometers. As the 
satellites fly over a feature on the Earth, such as an ice sheet or a moun- 
tain range, their separation changes because each satellite is affected 
differently by the gravitational attraction of the feature. 

The satellites use a microwave link to measure both the exact sepa- 
ration distance and its rate of change to an accuracy of better than 1 
micrometer per second, allowing a precise ’’snapshot” of the gravity field 
to be measured about every two weeks. This precise accuracy allows sci- 
entists to use the GRACE mission to ” weigh” various parts of the Earth 
system, learning about the distribution and changes in ocean mass, the 
growth or shrinking of the polar ice sheets, the amount of water in un- 
derground aquifers and other issues profoundly affecting climate change. 
The GRACE data are also an excellent complement to ocean height mea- 
surements provided by the radar altimeter on the TOPEX/POSEIDON 
satellite. 

GRACE will provide a framework for studying the gravitational sig- 
natures of gigantic, continent-sized underground water reservoirs, or 
aquifers. It also will provide a never-before-available perspective on global 
ocean circulation and the time variability of Earth’s overall external 
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shape, or geoid. This fundamental data set could enable great improve- 
ments in existing ocean radar altimetry data sets, and retrospective im- 
provements of seasonal to inter-annual climate change estimates. 

Plus d’information sur le site de NASA 
http://science.nasa.gov/headlines/y2001/ast30oct_l.htm 
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Figure 8.5 Free air gravity map based on model GGM02 including 
data from the GRACE satellite mission (Tap ley et al., 2005). Note 
that the major gravity anomalies are in narrow belts. 
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Figure 8.7 Map of the ages of the sea floor from Muller et al. (2008) 
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Figure 8.8 Global elevation map-Bathymetry of the sea floor. Note 
the correlation of continental elevation with crustal thickness and sea 
floor bathymetry with age. 



m 


-6000 -4000 -2000 0 2000 4000 6000 

elevation 


page facebook 


exosup.com 


8.8 Methodes 



9 


Evolution de la lithosphere oceanique. 
Mecanismes de la tectonique des plaques. 


9.1 Modeles devolution thermique de la lithosphere 

oceanique 

II est maintenant etabli que le flux de chaleur moyen dans les oceans 
(lOOmW m -2 ) est nettement plus eleve que celui des continents (56mW m -2 ). 
De plus, le flux de chaleur (Appendice A) et la bathymetrie varient 
de maniere systematique dans les oceans. La profondeur du plancher 
oceanique augmente et le flux de chaleur diminue avec la distance a 
la dorsale, c.a.d. avec l’age du plancher oceanique. Turcotte & Oxburgh 
[1967] avaient deja montre comment la convection dans le manteau pou- 
vaient expliquer les variations de flux de chaleur et la bathymetrie des 
fonds oceaniques. Par la suite, differents modeles de refroidissement 
d’une plaque par conduction thermique (en Appendice) ont ete proposes 
et expliquent bien les variations de bathymetrie et de flux de chaleur 
avec Page du plancher oceanique [ Sclater et Francheteau , 1970; Parsons 
& Sclater, 1977]. 


9.1.1 Variations du flux de chaleur avec Page 


Le modele fait appel au refroidissement de la lithosphere dont la temperature 
initiale est supposee egale a la temperature de fusion des basaltes T m « 
1300 °C et qui se refroidit par la surface au contact de l’eau des oceans a 
temperature « 0°C. Le flux de chaleur se calcule en resolvant Pequation 
de la chaleur dont la solution fait apparaitre des terrnes en \ft. On ob- 
tient ainsi les variations du flux de chaleur q en fonction de Page t du 
plancher oceanique. 


q(t) = 


kTrr 


\f-Knt 


(9.1) 
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ou k est la conductivity thermique (~ 3W m -1 K” 1 ), et k est la diffu- 
sivite thermique (k ~ 10“ 6 m 2 s _1 ). Pour t= 50 Ma, q ~ 55m W m~ 2 . 

En fait les observations du flux de chaleur sont tres bruitees et done 
difhciles a comparer au modele. II est clair que le flux ne peut pas etre 
oo a la dorsale. A proximite des dorsales, le flux varie beaucoup sur 
de tres petites distances. Ces variations sont dues a la circulation hy- 
drothermale [ Stein & Stein, 1994; Lowell et al., 1995]. En general, le 
flux moyen observe est inferieur a ce que predit le modele parce qu’on 
ne mesure que le flux conductif et que la plus grande partie de la chaleur 
est transportee par la circulation lrydrothermale. Le flux observe tend 
vers le modele a une certaine distance de la dorsale (pour un certain 
age du plancher oceanique). Cet age varie d’un bassin oceanique a un 
autre sans qu’on sache precisement pourquoi. L’hypothese la plus plau- 
sible est que le regime de contraintes controle l’ouverture des fractures a 
travers lesquelles se fait la circulation lrydrothermale. Pres des dorsales, 
le regime est extensionnel et les fractures restent ouvertes. A distance 
des dorsales, le regime devient compressif et referme les fractures. 

Les donnees de flux de chaleur s’ecartent egalement du modele a tres 
grande distance de la dorsale, i.e pour un age de > 80 Ma. Le flux tend 
vers une valeur a peu pres constante (48mW.m~ 2 ). L’explication la plus 
plausible de cette observation est que la lithosphere oceanique ne peut 
s’epaissir indefiniment. En effet la lithosphere qui se refroidit devient 
plus dense que le manteau. Cette instability va induire de la convection 
a petite echelle sous la lithosphere oceanique et maintenir l’epaisseur de 
la lithosphere et le flux a la base a peu pres constants. II faut noter que 
si la base de la lithosphere est definie comrne isotherme (ss 1300°C), il 
existe une relation entre le flux et l’epaisseur de la lithosphere. II est 
malheureusement tres difficile de mettre en evidence cette convection a 
petite echelle. (Sclater et al., 1980; McKenzie et al., 1980). 


9.1.2 Variation de la bathymetrie des oceans avec l’age 

Le modele de refroidissement de la lithosphere oceanique permet aussi 
de predire la bathymetrie. En effet, la densite moyenne de la lithosphere 
augmente a mesure qu’elle se refroidit (contraction thermique). L’equilibre 
isostatique implique que la surface de la plaque doit done s’enfoncer. En 
fait, on peut montrer que la subsidence thermique des fonds oceaniques 
est proportionnelle a la quantite de chaleur perdue depuis la formation 
du plancher oceanique. Cette quantite peut se calculer facilement en 
integrant le flux. On obtient que la bathymetrie h(t ) augmente propor- 
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Figure 9.1 Variations du flux conductif de chaleur en fonction de 
l’age du plancher oceanique. Les donnees representent le flux moyen 
pour chaque ocean. Le courbe theorique est calculee pour un modele 
de plaque. 



tionnellement a la racine carree de Page: 

h{ t) = ho + ^M (9.2) 

oil a est le coefficient d’expansion thermique (voir Appendice C). Le co- 
efficient d’expansion thermique des roches est de l’ordre de 3 x 10~ 5 A' _1 . 
Un changement de temperature de 330K donne done lieu a une expansion 
(contraction) de 1%. Cet effet est amplifie legerement par l’ajustement 
isostatique sous la charge supplementaire de la colonne d’eau de mer, ce 
qui donne: 

h(t) = h 0 + —m tt-\— ( 9 - 3 ) 

Pm Pw 2 V 7T 

oil p m est la densite du nranteau et p w celle de l’eau de mer. 

Contrairement au flux de chaleur, la bathymetrie moyenne observee 
suit tres bien le modele a proximite des dorsales et la profondeur des 
fonds oceaniques plus jeunes c^ue 80 Ma augmente proportionellement 
a la racine carree de l’age. Au dela de 80 Ma, la profondeur des fonds 
oceaniques reste constante. Cela implique un regime d’equilibre ther- 
mique (la lithosphere ne se refroidit pas) et le flux a la surface doit etre 
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Figure 9.2 Variations de la bathymetrie en fonction de la racine 
carree de Page. La courbe theorique est calculee pour un modele de 
plaque. 
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equilibre par le flux maintenu par la convection a petite echelle a la base 
de la lithosphere. 


9.1.3 Effet des points chauds 

Lorsque la lithosphere passe sur un point chaud clle est rechauffee. Ce 
rechauffement produit un ’’rajeunissement thermique” de la lithosphere. 
Apres le passage de la plaque au dessus du point chaud, la bathymetrie 
sera moins profonde et le flux de chaleur plus eleve que pour une plaque 
du meme age [ Crough , 1978; 1979]. Un profll bathymetrique le long de 
la chaine Hawaii-Empereur dans le Pacifique montre une augmentation 
de la profondeur du plancher oceanique au Nord-Ouest de Hawaii. 
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9.1.4 Applications 

Le modele de refroidissement de la lithosphere oceanique a marque le 
premier succes de la modelisation geodynamique. Ce modele simple ex- 
plique bien des variations a grande echelle de grandeurs observables. Par 
la suite, il a ete applique avec plus ou moins de bonheur a de nombreuses 
autres situations, telles la subsidence des marges continentales passives 
et celle des bassins sedimentaires intracontinentaux [e.g. Sleep & Snell, 

1976; McKenzie , 1978]. Ces modeles different par leurs conditions ini- 
tiales, mais font tous appel au merne mecanisme de refroidissement et 
de contraction thermique. 


9.2 Forces agissant sur les plaques 

II est possible de determiner certains mouvement verticaux en con- 
siderant les conditions d’equilibre isostatique. L’equilibre isostatique n’implique 
pas du tout que les contraintes disparaissent dans la lithosphere. 

Des contraintes importantes sont produites par les variations de la to- 
pographic (bathymetrie) associees a des variations de densite. Ces varia- 
tions sont induites par des differences de composition ou de temperature. 

Les contraintes qui en resultent sont transmises par et vont agir sur la 
lithosphere. La plupart de ces contraintes proviennent des frontieres des 
plaques. II est possible d’evaluer les forces qui agissent sur les plaques 
[. Forsyth & Uyeda, 1975]. 

9.2.1 Poussee de la dorsale (ridge push) 

A cause de l’elevation de la dorsale par rapport au plancher oceanique, 
la pression lithostatique a une profondeur donnee est plus elevee sous la 
dorsale que sous les bassins oceaniques. A cause de la difference de den- 
site cette difference s’attenue avec la profondeur et disparait au niveau 
de compensation isostatique. Ces differences de pression induisent une 
poussee de la dorsale (pression elevee) vers les bassins oceaniques. En 
supposant une difference de bathymetrie de ~ 3 km, la poussee moyenne 
de la dorsale est de l’ordre de 30 MPa. 

9.2.2 Traction de la zone de subduction (slab pull) 

Les anomalies gravimetriques de grande longueur d’onde sont positives 
au dela des zones de subduction. Ces anomalies positives suggerent que 
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Figure 9.3 Difference de pression en fonction de la profondeur entre 
la dorsale et les bassins oceaniques. La difference de bathymetrie est 
3km. II y a equilibre de pression a 100km (base de la plaque). La 
pression est plus elevee sous la dorsale et cause une poussee vers les 
bassins. 


Ridge push for A z = 3km 



la plaque froide est plus dense que le manteau avoisinant. Si la plaque 
oceanique est agee et froide, sa temperature moyenne de l’ordre de 
600 °C est «600 K plus basse que celle du manteau. Cette difference 
de temperature implique que la densite moyenne de la plaque est ss 
0.05Mg m -3 (1.8%) plus elevee que celle du manteau. Considerant que 
la longueur de la plaque est 500 km, la contrainte moyenne induite est 
de l’ordre de 250 MPa {gApL). Cette contrainte a pour effet de tirer 
vers le bas la plaque subductee. Elle est au moins 10 fois plus elevee que 
la poussee des dorsales. 

II faut noter que des plaques jeunes peuvent aussi etre entrainees dans 
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les zones de subduction. C’est le cas notament de la plaque de Juan 
de Fuca a l’ouest du Canada et de la plaque des Cocos, a l’ouest de 
l’Amerique du Sud. Dans ces deux cas, l’inclinaison de la subduction est 
faible. 

Les profils bathymetriques montrent un bombement (1km) de la plaque 
qui retourne dans le manteau en amont de la zone de subduction. Ce 
bombement, sur une distance de 500 a 1000km, s’explique par la rigidite 
de la plaque froide. On pourrait faire l'analogie avec une poutre fixee 
d’un cote et sur laquelle s’exerce une charge verticale; la poutre est 
pliee par la charge. La distance caracteristique depend des proprietes 
mecaniques de la plaque et du couple de forces dues aux differences de 
densite. Les contraintes sont tensiles pres de la surface de la plaque. 

Souvent, il y a de l’extension dans la plaque qui chevauche, tres en aval 
de la zone de subduction (extension d’arriere arc) avec parfois formation 
de bassins d’arriere arc. L’explication la plus courante est que la plaque 
qui sombre dans le manteau rapidement suce et tire vers l’amont la 
zone de subduction et la plaque chevauchante. II y aura done un recul 
progressif de la zone de subduction. 


9.2.3 Forces sur les failles transformantes 

Ces contraintes sont tres difficiles a evaluer. On ne sait pas tres bien 
quelle est l’importance de la friction le long des failles transformantes. 
Pour les failles transformantes a proximite des dorsales, la lithosphere est 
chaude et peu epaisse. A faible profondeur, on atteint un regime ductile 
avec une viscosite effective et done une friction assez faible. Dans ce cas, 
la friction est vraissemblablement negligeable. Par contre, on ne sait pas 
ce qu’il en est pour les grandes failles transformantes comme la faille de 
San Andreas on la faille Alpine. Certaines donnees, notamment l’absence 
d’anomalie thermique causee par une friction importante, suggerent un 
tres faible couplage le long de la faille, qui serait peut-etre ” lubrifiee” . 


9.2.4 Friction a la base de la plaque (viscous drag) 

Cette force est produite par le mouvement relatif de la lithosphere par 
rapport a l’asthenosphere. Elle depend de la vitesse relative, de la vis- 
cosite de l’asthenosphere. Elle est probablement assez faible ss 10 MPa. 
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9.2 Forces agissant sur les plaques 

9.2.5 Contraintes intraplaques 

La figure 9.4 montre la direction des contraintes observees a l’interieur 
des plaques. Ces observations sont basees sur differ entes series de donnees: 

(1) Les mecanismes aux foyers des tremblements de terre, (2) des mesures 
in situ par hydrofracturation dans des forages, et (3) les donnees neotectoniques 
(mouvement recent de failles). 

Bien que la traction de la zone de subduction soit la plus impor- 
tante de toutes les forces qui agissent sur les plaques, la poussee de la 
dorsale joue un role important, parce qu’elle controle en grande par- 
tie les contraintes a l’interieur des plaques. 1 Ainsi les contraintes dans la 
plaque nord Americaine, a l’ouest des montagnes rocheuses sont dues a la 
poussee de la dorsale mid-atlantique. En premiere approximation, l’effet 
des autres forces est negligeable. Les contraintes induites a l’interieur 
des plaques ont ete modelisees [Solomon et al., 1975] et comparees avec 
les contraintes estimees par la geophysique. 


1 II faut aussi noter que les contraintes compressives sont mieux transmises que les 
contraintes tensiles. 
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Figure 9.4 Les contraintes estimees a l’interieur des plaques. La di- 
rection des contraintes est determinees par les mecanismes au foyer 
de tremblements de terre, l’activite neotectonique et des mesures in 
situ par hydrofracturation dans les forages. 



9.2.6 Stability des chaines de montagne 

Les calculs des contraintes dues a la hauteur des dorsales relativement 
aux bassins oceaniques montrent comment des differences de topogra- 
phie merne compensees isostatiquement, engendrent des contraintes im- 
portantes. Le meme type de calcul pourrait s’appliquer aux chaines de 
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montagnes. Les differences topographiques impliquent d’enormes con- 
traintes qui doivent etre equilibrees par les contraintes tectoniques. Une 
disparition des contraintes tectoniques pourrait ainsi declencher un col- 
lapse de la chaine par extension. Ce genre de collapse a ete observe dans 
de nombreuses chaines. Ainsi la Cordillere en Colombie Britanique sem- 
ble avoir subi une extension rapide a 50 Ma. Un tel collapse pourrait 
egalement expliquer l’extension presente dans le Basin and Range. Dans 
la Mediterrannee, il y a ainsi juxtaposition de zones d’extension et de 
zones de compression. 


9.3 Mecanismes de la Tectonique des plaques 

II est acquis que la tectonique des plaques est etroitement liee a la convec- 
tion thermique dans le manteau. Les premiers modeles avaient associe la 
tectonique des plaques a de grandes cellules de convection dans le man- 
teau [ Runcorn , 1964]. Les dorsales etaient associees aux courants chauds 
ascendants et les zones de subduction aux courants descendants. II est 
tout a fait certain que la convection dans le manteau opere de fagon 
beaucoup moins reguliere. La geometrie des dorsales et des zones de 
subduction est incompatible avec une convection organisee en quelques 
grandes cellules. Par ailleurs, les plaques ne sont pas passives rnais par- 
ticipent a la convection. L’existence des points chauds est bien etablie 
et Morgan [1972] a propose qu’ils jouent le role essentiel dans le re- 
froidissement du manteau par convection. II est maintenant clair que 
les points chauds ne contribuent qu’a <20% du flux de chaleur total. 
De nombreuses questions restent en suspens: la forme de la convection 
dans le manteau, le couplage entre les plaques et la convection, l’origine 
des points chauds, le recyclage dans le manteau des plaques oceaniques 
sub duct ees. 

La cause de la convection thermique est le refroidissement seculaire 
du manteau et du noyau. L’ energie (chaleur) du manteau et du noyau 
proviennent de l’accretion initiale de la Terre, de la differentiation du 
noyau, et de la radioactivite dans le manteau. Pour mesurer l’intensite 
de la convection, les physiciens utilisent le nombre de Rayleigh Ra (voir 
Appendice D). Le nombre (sans dimension) de Rayleigh mesure le rap- 
port des forces dues aux differences de densite aux forces de friction 
(viscosite), ou le rapport du temps de refroidissement par conduction 
thermique au temps par convection. La convection n’est possible que 
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si Ra depasse une valeur critique ( Ra c ) qui depend des conditions aux 
limites mais qui est de l’ordre de 1000. 

Pour calculer le nombre de Rayleigh du manteau et determiner si la 
convection est possible, il faut determiner la viscosite du manteau. Les 
mesures en laboratoire de la viscosite des roches dans les conditions de 
pression et temperature du manteau superieur sont tres difficiles. Par 
contre, il a ete possible depuis tres longtemps d’estimer la viscosite en 
mesurant le rebond post glaciaire. Le temps de retour a l’equilibre est 
proportionnel a la viscosite. Les premiers calculs de Haskell [1937] qui 
avait propose une viscosite de « 10 20 Pa s ont ete continues. Les etudes 
experimentales de la viscosite effective aux temperatures et pressions du 
manteau pour des roches comme la dunite ou de mineraux tels l’olivine 
donnent des valeurs comparables a cello de Haskell. 

Les meilleures estimations de Ra pour le manteau sont maintenant 
~ 10 6 — > 10 7 . Les experiences en laboratoire et des calculs numeriques 
montrent que dans ce cas, la convection n’est pas organisee en cellules 
mais se fait de fagon chaotique, avec des panaches qui se ferment et 
montent rapidement et des gouttes denses qui sombrent. 

9.3.1 Organisation de la convection dans le manteau 

Pendant longtemps, la question du role de la zone de transition revelee 
par la sismique entre 400 et 600 km dans le manteau a ete au centre 
des debats. Certains geophysiciens pensaient qu’il y a deux niveaux de 
convection dans le manteau et que la zone de transition etait une barriere 
a la convection. La majorite des geophysiciens pensent maintenant que 
la convection traverse tout le manteau. Malheureusement, les donnees 
disponibles restent contradictoires et ne permettent pas d’exclure les 
autres modeles. 
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Figure 9.5 Un modele de convection organisee dans le manteau. Ce 
modele est peu plausible. 



Les questions restent en suspens mais tout modele doit etre compatible 

avec un certain nornbre de donnees maintenant bien etablies. 

• Les geochimistes pensent que les differences entre basaltes de dorsales 
(MORB) et basaltes qui proviennent de points chauds exigent deux 
reservoirs distincts dans le manteau [O'Nions et ai, 1980]. Ceci im- 
plique qu’il n’y a pas de melange, ce qui semble incompatible avec la 
convection a travers tout le manteau 

• Certaines donnees sismiques a proximite des zones de subduction 
indiquent que les plaques penetrent profondement dans le manteau 
inferieur et certaines plaques pourraient meme etre suivie jusqu’a la 
base du manteau. Par contre, on a egalement vu par la sismique que 
certaines plaques semblent stagner au dessus de la zone de transition 
[van der Hilst et ai, 1997]. Les donnees sismiques sur la structure pro- 
fonde des zones de subduction sont resumees dans Particle de Fukua et 
al. [2001]. Un modle de penetration intermittente avec des avalanches 
dans le manteau inferieur a ete propose par [ Machetel et Weber 1991]. 

• Les experiences faites a haute pression permettent de verifier les pro- 
prieties physiques des mineraux dans les conditions de pression et 
temperature du manteau inferieur. Certaines interpretations de ces 
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experiences suggerent des differences de composition entre le nranteau 
superieur et inferieur. De telles differences sont egalement incompati- 
bles avec le melange que produirait la convection. 

• L’anisotropie de la vitesse des ondes sismiques dans le nranteau est 
due a l’orientation de l’olivine dans le champ de contraintes. Cette 
anisotropie est directenrent correlee avec la direction des courants de 
convection et permet done de determiner la direction des lignes de 
courant dans le manteau. 

• Depuis une dizaine d’annees, on peut mesurer des variations laterales 
des vitesses sismiques dans le manteau, par la technique de tomogra- 
phie sismique. L’interpretation la plus courante de la tomographie sis- 
mique est que les variations de vitesse sismique sont dues a des varia- 
tions de temperature. La vitesse des ondes sismiques, particulierement 
cello des ondes S, decroit lorsque la temperature augmente. La tomo- 
graphie sismique permet ainsi de determiner des variations laterales de 
temperature dans le manteau. Dans l’ensemble, les resultats de tomo- 
graphie sismique montrent que des anomalies de vitesse sismique sont 
correlees avec la geologie et la tectonique de surface. Ainsi, les cratons 
ont de profondes racines froides. Les zones de subduction apparaissent 
froides. Le manteau est chaud sous les dorsales, en particulier VEast 
Pacific Rise, et les points chauds. Certaines de ces anomalies ont une 
source profonde dans le manteau et semblent continuer jusqu’a la zone 
dite D” a la base du manteau. C’est le cas d’un point chaud sous le 
continent africain et de certaines zones de subduction (Grand et al., 
1997). 

• Les variations profondes de la densite se refletent mieux dans les 
anomalies du geoi'de. L’effet de la zone de transition sur les anomalies 
depend du contraste de viscosite. La correlation entre anomalies de 
tres grande longueurs d’onde du geoi'de et la vitesse des ondes P dans 
le manteau inferieur est compatible avec la convection a travers tout 
le manteau. 

• Depuis quelques annees, les geophysiciens envisagent un nouveau modele, 
qui suppose que le manteau inferieur serait compositionnellement different, 
plus dense, et se melange mal avec le reste du manteau. Les analyses 
recentes de donnees sismiques semblent supporter ce modele (Kellogg 

et al., 1999). La figure 9.7 montre le genre de structure de convection 
est induite dans un manteau avec un couche inferieure dense. 

• De nombreux travaux recents cherchent a etablir de fagon plus precise 
le lien entre la tectonique des plaques et la convection dans le man- 
teau (voir p. ex. les articles dans la monographie de Richards et al. 


exosup.com 


page facebook 


9.3 Mecanismes de la Tectonique des plaques 


125 


Figure 9.6 Un modele de convection dans le manteau avec points 
chauds et zones de subduction. Ce modele n’implique pas de rela- 
tion directe entre les dorsales et les structures dues a la convection 
profonde dans le manteau. 


Hoi spot 



[2000]). Les methodes numeriques de plus en plus rafinees permettent 
mettant d’introduire des plaques dans les calculs de convection dans 
le manteau. 
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Figure 9.7 Un modele du manteau avec un reservoir profond plus 
dense qui n’est pas entraine dans la convection. 



9.4 Flux de chaleur 

Le flux de chaleur mesure la quantity denergie qui traverse la surface 
de la terre par unite de temps. II se mesure en W m~ 2 . Sa mesure 
permet de faire le bilan d’energie de la Terre. Les variations locales et 
regionales peuvent egalement etre interpretees. Le flux moyen terrestre 
est ss 80mlT m~ 2 . 

En principe, la chaleur peut etre transportee par conduction ther- 
rnique. II y a transport denergie sans transport de matiere, ou par con- 
vection, l’energie est transportee par la matiere en mouvement. D’apres 
la loi de Fourier, le flux par conduction est proportionnel au gradient de 
temperature: 


q = -fcVT 

(9.4) 

dT 

(9.5) 

q z = -k— 

oz 


La constante de proportionnalite k est la conductivity thermique, une 
propriete physique caracteristique de chaque materiau. La conductivity 
thermique des roches est de l’ordre de 3 W to -1 K -1 . 

Pour mesurer le flux de chaleur, il faut done mesurer le gradient 
de temperature et la conductivity thermique des roches. Sur terre, les 
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mesures se font dans des forages. Dans les oceans, on fait les mesures a 
l’aide d’une sonde qui penetre (« 10 m) dans les sediments. On mesure 
le gradient de temperature et la densite des sediments qui permet de 
determiner leur conductivite. 


9.5 Modele de refroidissement de la lithosphere 

Le plancher oceanique se forme a la dorsale (x = 0) ou la temperature de 
la roche est approximativement la temperature de fusion T m (independament 
de la profondeur). Le plancher oceanique est entraine et s’eloigne de la 
dorsale. II se refroidit a partir de la surface qui, par contact avec Pocean, 
est maintenue a temperature T ss 0. 

En regime de quasi equilibre, la temperature en un point fixe de 
l’espace ne change pas au cours du temps (la seule variation est celle qui 
accompagne Pentrainement du plancher oceanique). La temperature est 
determin.ee par une equation aux derivees partielles similaire a P equation 
de la chaleur 


d 2 T dT 
K ~d^~ V ~d^ 


(9.6) 


ou k est la diffusivite thermique (m 2 s x ). La solution de cette equation 
est de la forme: 


T(x, z) = T m erf 
avec la fonction d’erreur 


—\ — 

2 V kx . 


2 r z 

erf (z) = —= / exp(— u 2 )dv 

Jo 


D’ou le flux de chaleur q est calcule: 


(9.7) 


(9.8) 


oil k est la conductivite thermique (W m -1 
du plancher oceanique. 



K 1 ) et t = x/v est 


(9.9) 

Page 
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9.6 Subsidence thermique 


Le coefficient d’expansion thermique a est la variation relative de volume 
due a une variation de temperature de 1 K: 


AV 

V 


= aAT 


^ = —aAT 


p(T) = Po (l - aT) 


Le coefficient d’expansion thermique a est exprime en iL -1 . L’equilibre 
isostatique sous la lithsophere oceanique implique que la pression d’une 
colonne de roche du bassin oceanique soit egal a la pression sous la 
dorsale. Cela requiert une difference de hauteur Ah: 



(9.10) 

(9.11) 

(9.12) 


oil C est la chaleur specifique de la lithosphere. En integrant la flux de 
chaleur on obtient done: 


akT m Vt 

2pC^TTK 



(9.13) 


9.7 Convection thermique 
9.7.1 Viscosite d’un fluide 

Un fluide en mouvement est sounds a des contraintes qui sont induites 
par les differences de vitesse a l’interieur du fluide. Ces contraintes ten- 
dent a attenuer les differences de vitesses. II est demontre experimentalement 
que pour de nombreux fluides, la contrainte est proportionnelle au gradi- 
ent de vitesse. La constante de proportionalite est une proprite physique 
du fluide, la viscosite. Si v est la vitesse du fluide, la contrainte T est: 

T = /zVv (9.14) 

oil /i est la viscosite. La viscosite est mesuree en Pa s. (L’unite cgs est 
la poise: 1 poise = 0.1 Pa s). La viscosite de l’eau est de O.OOlPa s. La 
viscosite du manteau terrestre est de l’ordre de 10 2O Pa s. La viscosite 
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du manteau terrestre a pu etre determinee grace a l’analyse du rebond 
post-glaciaire. Le temps necessaire pour que les regions couvertes par 
les glaciers remontent a leur niveau d’equilibre est proportionnel a la 
viscosite. 


9.7.2 Nombre de Rayleigh critique 

Pour maintenir la convection thermique dans un fluide, il faut que les 
forces dues aux differences de densite soient superieures aux forces dues 
a la viscosite. Le rapport des differentes forces agissant sur le fluide 
determine le nombre (sans dimension) de Rayleigh: 

Ra = gpaATV 

KfJ, 

g est l’acceleration due a la gravite, p est la densite, a est le coefficient 
d’expansion thermique du fluide, AT est la difference de temperature 
entre les surfaces superieure et inferieure du fluide, L est l’epaisseur de 
la couche de fluide, n est la diffusivite thermique, p est la viscosite. Une 
autre interpretation du nombre de Rayleigh est le rapport du temps car- 
acteristique pour transport de chaleur par conduction au temps necessaire 
pour transporter la chaleur par convection. Lorsque le nombre de Rayleigh 
est eleve, la convection est done plus efficace que la conduction pour re- 
froidir le fluide. La convection est possible des que le nombre de Rayleigh 
depasse une valeur critique (600 < Ra c < 2, 000) qui depend de la 
geometrie du fluide. 

Dans la convection de Rayleigh-Benard (fluide chauffe par le bas), au 
nombre de Rayleigh critique, la cellule de convection ne developpe qu’une 
seule longueur d’onde (proportionnelle a l’epaisseur du fluide. Lorsque 
Ra augmente, les celules de convection peuvent se developper avec un 
nombre croissant de longueurs d’onde. Pour Ra » Ra c , la convection 
est possible a toutes les longueurs d’onde. Cela donne done lieu a un style 
de convection chaotique. Le nombre de Rayleigh du manteau terrestre 
est probablement superieur a 1,000,000. Dans ce cas, la convection sera 
chaotique et non pas organisee en cellules de longueur d’onde fixe. 
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Figure 9.8 Les anomalies de vitesse determinee par tomographie sis- 
mique dans le manteau superieur et dans le manteau inferieur dans 
la zone D” . Les anomalies represente la variation de la vitesse des on- 
des S a une profondeur donnee par rapport a un modele de reference. 
Les couleurs chaudes correspondent aux zones a faible vitesse et les 
couleurs froides aux zones a vitesse elevee. 

(a) Dans le manteau superieur, les anomalies chaudes correspondent 
aux dorsales, au rift est africain, au Basin & Range. La signature 
froide des cratons Archeens est preservee jusqu’a plus de 250km. 

(b) Dans le manteau inferieur, les contrastes de vitesse sont impor- 
tants. Les zones chaudes dans le manteau inferieur pourraient ali- 
menter les points chauds (Rift africain et Pacifique). 
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Figure 9.10 Shear wave travel time anomalies for the 60 to 300knr 
slice. The travel times are calculated from the CUB2.0 model of 
Shapiro and Ritzwoller (2002), and the anomalies are relative to the 
world average. The value of 60km for the slice upper limit is selected 
to exclude the crust which is not very well resolved by the seismic 
model. 
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Figure 9.11 Shear wave velocity anomalies at 350km depth. The 
anomalies relative to the mean shear wave velocity at that depth 
are given in %. By convention, the color scale is inverted as negative 
velocity anomalies are interpreted as hot. The seismic model is that 
of Grand et al. (1997). 
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Figure 9.12 Shear wave velocity anomalies at 1810km depth. Note 
the reduced scale relative to Figure 9.11. 
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Figure 9.13 Shear wave velocity anomalies in the D” region at 
2800km depth. Anomalies are enhanced relative to lower mantle. 
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Rifts intracontinentaux. 

Marges continentales passives. Bassins 
sediment air es. 


10.1 Rifts intracontinentaux 

La rupture d’un continent et la naissance d’un bassin oceanique est tou- 
jours precedee d’un episode d’extension pendant lequel la lithosphere et 
la croute continentales sont amincies. Les rifts intracontinentaux sont 
des zones etroites (~ 100km) et alongees (> 1000km) d’extension a 
l’interieur d’un continent. Les rifts sont caracterises par une zone d’extension 
tres localisee autour d’une vallee centrale encadree par des failles nor- 
rnales. Les rifts actifs voient une forte activite sismique et volcanique. 

Le rift est africain, aujourd’hui actif, est l’archetype du rift intraconti- 
nental. Le rift du Rio Grande qui traverse le Nouveau Mexique et une 
partie du Colorado, le rift de Baikal dans l’est de la Siberie et le graben 
du Rhin sont des rifts contemporains, moins actifs que le rift est africain. 

Les rifts actifs ne peuvent etre isoles du contexte general de la tectonique 
des plaques, mais il n’y pas de relation universelle entre la tectonique 
regionale et la formation d’un rift. 

Le rift intracontental se distingue d’un bassin oceanique en formation 
(tel la Mer Rouge) par la nature de la croute qui reste continentale. 
Certains rifts evoluent pour former un bassin oceanique, mais il y a 
a l’interieur des continents de nombreuses traces laissees par des rifts 
avortes. 


10.1.1 Comparaison des rifts Est Africain et du Rio 

Grande 

Le rift de Rio Grande (RGR) oriente NS traverse le Nouveau Mexique 
et le Colorado. Au Sud, le RGR rejoint la province du Basin & Range 
(B&R) qui est une large (> 500km) zone d’extension. Le RGR est separe 
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Figure 10.1 Le rift est africain. 



du B&R par le Plateau du Colorado. La presence du Grand Canyon qui 
traverse le plateau montre que celui ci a ete recemment souleve sans subir 
d’autre deformation. Le Plateau du Colorado est ~ 1000m au dessus du 
B&R. II y a eu par contre du volcanisme depuis 30Ma aux marges du 
plateau. Les indicateurs de contrainte indiquent qu’il est en compression. 

Le rift est africain (EAR) est connecte au rift de la Mer Rouge et au 
rift du golfe d’Aden. II part du point triple des Afars au Nord, traverse 
l’Ethiopie, le Kenya et finit en Tanzanie. Au Kenya il se separe en 2 
branches autour du craton de Tanzanie. La branche de l’Ouest suit les 
grands lacs africains (Kivu, Tanganyika, Victoria). Le rift est aujourd’hui 
eleve (1500-2000m), mais le soulevement est recent et jusqu’a 10 Ma, le 
rift etait au niveau de la mer, tant en Ethiopie qu’au Kenya. 
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Figure 10.2 La seismicite de 1’Afrique 

Seismicity of Africa: 1977 - 1997 
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L’EAR et le RGR sont actifs depuis ca 50 Ma. Leur evolution est assez 
differente (voir Appendice). Aujourd’hui, le rift africain est nettement 
plus actif avec du volcanisme tres recent (<10ka) et une forte seismicite. 
Les mecanismes au foyer sont d’extension. 

Le rift de Baikal, situe en Siberie a 2000km au Nord Est de la chaine 
himalayenne, est plus jeune (< 20Ma). II est actif sismiquement. L’axe 
du rift coincide avec le Lac Baikal. 

Le graben du Rhin suit le cours du Rhin depuis les Alpes jusqu’a 
la Mer du Nord. L’activie sismique est moderee. Le volcanisme le plus 
recent date de 40ka. 

Les anomalies de gravite et les etudes par refraction sismique montrent 
que, sous tous les rifts, la croute est tres amincie (25-30 km). La vitesse 
des ondes sismiques est anormalement faible dans le manteau avec une 
vitesse des ondes P ( v p < 7.9 km s _1 ). Un faible vitesse suggere une 
temperature elevee. Du point de vus de la gravite, il y a superposition 
de: (1) une anomalie negative de tres courte longueur d’onde due aux 
sediments de la vallee du rift, (2) une anomalie positive de moyenne 
longueur d’onde due a Pamincissement de la croute, et (3) une anoma- 
lie negative de tres grande longueur d’onde due a la faible densite du 
manteau. 

Le flux de chaleur est tres variable dans la vallee centrale du rift a 
cause de la circulation hydrothermale dans les sediments. En moyenne, 
il est eleve (> 80mkEro -2 ). Le flux de chaleur reste souvent faible dans 
les regions non deformees qui entourent le rift. Ainsi, il est < 40mkEm -2 
dans le craton en Tanzanie, et ss 50mkETO -2 dans le plateau du Colorado. 

Dans le EAR on a observe une progression temporelle de la compo- 
sition des roches volcaniques. La petrologie indique que la source des 
magmas est de moins en moins profonde avec le temps. La remontee 
de la source des magmas est relativement rapide (de 150 a 40 km de 
profondeur en 20 Ma). 

10.1.2 Mecanisme de formation d’un graben 

Dans les rifts, les grabens sont encadres de maniere symetrique par des 
failles normales. Des qu’une epaule du rift est soulevee, les contraintes 
induites dans la croute, amplifiees par l’effet des sediments dans le demi 
graben, provoquent la formation d’un faille symetrique a une certaine 
distance (40-80 km) de la premiere faille. Le poids des sediments va 
contribuer a l’enfoncement du graben central. 

Les zones d’extension larges, telles le B&R, sont traversees par une 
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Figure 10.3 Les volcans actifs en Afrique 


faille de detachement a l’echellc de la lithosphere. De telles failles provo- 
quent une asymetrie des deformations de la croute superieure. On parle 
de cisaillement pur pour un rift symetrique et de cisaillement simple 
pour 1’ extension le long d’une faille de detachement a l’echelle lithospherique 
[Wernicke, 1985]. L’extension le long de la faille de detachement permet 
la remontee de croute profonde et formation de metamorphic core com- 
plexes. 

II faut bien noter que l’extension totale dans un rift intracontinental 
est tres limitee. Elle peut s’estimer grossierement a partir de l’amincissement 
de la croute qui est rarement superieur a 25%. Cela impliquerait une ex- 
tension du meme ordre sur une distance «50 km. Par contraste, l’amincissement 
de la croute dans le B&R est plus important et l’extension est de l’ordre 
de 50 a 100% sur des distances de plusieurs centaines de km. 

10.1.3 Modeles devolution 

Le mecanisme de formation d’un rift intracontinental est loin d’etre 
elucide. Dans les mecanismes actifs, le rift est cause par l’interaction 
de la lithosphere avec une asthenosphere anormale ou un point clraud. 

La pnetration de l’asthenosphere dans la lithosphere provoque du vol- 
canisme, le soulevement et l’extension de la croute. Dans les mecanismes 
passifs, ce sont les contraintes tectoniques qui provoquent l’extension et 
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la remontee adiabatique de Pasthenosphere. La fusion de Pasthenosphere 
et le volcanisme sont dus a la decompression adiabatique. 

Le contexte tectonique regional est tres different pour le RGR et pour 
PEAR et il n’est pas du tout evident que le meme mecanisme a joue dans 
les deux cas. Le RGR et Pextension du B&R sont lies a la subduction 
de la plaque de Farallon. Aujourd’hui, les contraintes d’extension pour- 
raient etre induites par le mouvement de la Plaque Pacifique le long de la 
faille de San Andreas. Pour PEAR, il est impossible d’expliquer les con- 
traintes extensionnelles par le mouvement des plaques. Au contraire, les 
dorsales qui entourent la plaque africaine devraient induire un regime de 
compression. Par ailleurs, PEAR part du point triple des Afars, un point 
chaud. La tomographie sismique semble indiquer que ce point chaud est 
enracine profondement dans le manteau. 

Le graben du Rhin est situe juste au Nord des Alpes et le rift de Baikal 
a ss 1500km au Nord de l’Himalaya. Leur formation est probablement liee 
aux collision entre l’Afrique et l’Europe et entre l’lnde et l’Asie. Dans la 
cas de l’Asie, les deformations dues a la collision sont tres profondes et 
affectent les contraintes malgre la distance de Baikal a PHimalaya. Dans 
la cas de tels rifts, on parle parfois d ’ impactogene. 

Il semble que le mode d’extension depend beaucoup du flux de chaleur. 
Les contraintes sont plus localisees et Pextension se limite a un rift 
si le flux de chaleur est faible (< 60m Wm~ 2 ). Pour un flux eleve (> 
80mPP?n -2 ), Pextension se fait dans une region tres large (e.g. B&R) 
[Buck, 1991], 


10.1.4 Rift oceanique 

Si Pextension du rift se poursuit a taux eleve, elle ne peut pas etre 
compensee par Pamincissement de la croute continentale et il va y avoir 
formation de plancher oceanique. Dans la Mer Rouge, une croute de 
type oceanique se forme depuis « 20 Ma. L’ouverture est tres lente et les 
anomalies magnetiques marines peu marquees. Les etudes recentes dans 
le Golfe d’Aden et la corne de PAfrique (Djibouti, Somalie) indiquent que 
Pouverture du Golfe d’Aden est posterieur a celle de la Mer Rouge. Cette 
ouverture a commence vers lOMa. Dans le Golfe d’Aden, Les anomalies 
magnetiques marines sont marquees et elles montrent la propagation du 
rift a partir de la dorsale de Carlsberg vers la pointe de PAfrique. La 
direction des failles transformantes montre que la vitesse de la plaque 
arabe par rapport a PAfrique est dirigee ENE. 
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10.2 Rifts avortes. Aulacogenes. 

10.2.1 Rifts avortes en Amerique du Nord 

De nombreux rifts intracontinentaux n’ont pas evolue pour donner nais- 
sance a un bassin oceanique. Ils laissent des traces dans les continents. 

La plupart des fleuves suivent d’anciens rifts (Saint Laurent, Missis- 
sipi, Rhin, Rhone). Les rifts avortes sont souvent associes a des zones 
de faiblesse ou est focalisee l’activite sismique a l’interieur des plaques 
(Reelfoot, St Laurent). Ce sont les anomalies de gravite et magnetiques 
qui revelent la presence de tels rifts avortes. De nombreux rifts fos- 
siles ont ainsi ete identifies en Amerique du Nord: Keeweenawan, Saint- 
Laurent, Ottawa-Bonnechere, Lac Champlain, Reelfoot, Richmond Gulf, 
McKenzie Delta. 

Le rift de Keeweenawan, sous les grands lacs d’Amerique est associe 
a la plus forte anomalie positive de gravite (200 mGal) sur le continent 
nord americain. De fortes anomalies magnetiques sont dues a la presence 
de roches basaltiques. Les etudes sismiques (profils GLIMPCE) montrent 
que ce rift avorte contient des coulees successives de laves d’une epaisseur 
totale > 15 km [Cannon et al., 1989]. 


10.2.2 Aulacogenes 

Les aulacogenes sont des bassins sedimentaires tres profonds (~ 10km) 
de forme allongee. Ils se forment souvent a partir d’un des trois rifts 
qui se developpent autour d’un point triple. De tels aulacogenes se trou- 
vent de part et d’autre de l’Atlantique. Ils coincident avec l’embouchure 
des grands fleuves et les provinces riches en hydrocarbures. A l’interieur 
du continent nord americain, d’autres aulacogenes (Bassin d’Anadarko, 
Southern Oklahoma aulacogen) se trouvent sur l’ancienne marge de 
l’Amerique du Nord. 


10.3 Marges continentales passives. Formation des 
bassins sedimentaires. 

10.3.1 Subsidence des marges continentales 

L ’evolution des marges continentales passives est determinee en partie 
par l’amincissement de la croute et les ajustements isostatiques durant 
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la phase de rifting, et en partie par le refroidissement, la contraction 
thermique et la subsidence pendant la phase de drifting. 

• Mecanismes devolution: subsidence thermique, flexure de la lithosphere. 

• Exemples: Marge continental a Test du Canada, Golfe de Biscaye. 

10.3.2 Bassins sedimentaires intracontinentaux 

La plupart des nrodeles proposes pour expliquer la subsidence des bassins 
intracratoniques supposent qu’ils sont le resultat du refroidissement et 
de la contraction thermique d’une croute anormalement chaude initiale- 
ment. Cette condition initiale est le resultat de l’activite d’un point 
chaud ou d’un episode d’extension. 

Le modele de McKenzie [1978] examine revolution d’un bassin apres 
un episode d’extension de la lithosphere. La condition initiale n’explique 
pas la symetrie circulaire des bassins sedimentaires. Cette symetrie est le 
resultat de la subsidence due a l’ajustement isostatique sous le poids des 
sediments. Le modele n’est pas universellement applicable. II n’y a pas 
de trace de rift sous le bassin de Williston ou la croute est epaissie. Le 
rift de Keeweenawan est beaucoup plus vieux que le bassin de Michigan 
et ne peut done pas avoir cause sa formation. Par contre, le modele 
semble bien s’appliquer au bassin de la Mer du Nord ainsi qu’au bassin 
de Paris. 

II faut remarquer que les bassins sedimentaires du continent nord 
americain (Wiiliston, Michigan, Illinois, Baie d’Hudson) ont debute leur 
histoire a peu pres a la meme periode (400-500 Ma) mais que la duree de 
la subsidence et la quantite de sediments accumules varient d’un bassin 
a l’autre. 

Les bassins dits d’avant-pays ( foreland basins), tel le bassin d’ Alberta, 
ont subi une evolution tres differente. La formation du bassin d’Alberta 
coincide avec celle des montagnes rocheuses qui ont ete transporters sur 
le continent Nord Americain. C’est la charge des montagnes rocheuses 
qui a provoque la flexure de la plaque nord-americaine et la subsidence du 
bassin. Cette subsidence est relativement bien expliquee par les modeles 
de plaque elastique [Beaumont, 1981]. 
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Figure 10.4 La topographie des zones d’extension au SW des USA 



Figure 10.5 Le point triple des Afars 
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10-4 Chronologie des rifts de VAfrique de I’Est et du Rio-GranddAb 

10.4 Chronologie des rifts de l’Afrique de l’Est et du 

Rio- Grande 

10.4.1 Rio-Grande 


Age (Ma) Rio-Grande 


Western American margin 


40-30 calalkaline volcanism 

32 rhyolitic volcanism 

30-20 Bimodal volcanism 

(rhyolites-basalts) 


Subduction of Farallon 
plate. 

Ridge off California 
is subducted. 

San- Andreas transform 
fault begins. 


15 

Basalts 


Growth of San Andreas. 

First uplift of Colorado Plateau 

7-4 

Uplift and graben formation 


Second uplift of Colorado 
Plateau. 

10.4.2 East Africa 

Age (Ma) 

Kenya 

Age (Ma) Ethiopie 

50-30 

basalts 

45 

Flood basalts 

21-14 

Flood basalts 

16 

basalts 

15-11 

Regional Uplift (500m) 

12 

rhyolites 

10 

Western shoulder uplift 

8 

Faulting 

7-5 

Alkaline basalts 



6-2 

Flood lavas 

8-0.5 

Flood basalts 
(2 episodes) 

4 

Main uplift (1,500m) 



3 

Dykes 



2 

Normal faulting (graben) 

1 

Plateau uplift 

1-present 

Volcanism-Faulting 
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Zones de convergence 1: Subdnction 


C’est dans les zones de convergence que se forme et s’accumule la plus 
grande partie de la croute continentale. La future croute continentale 
se forme essentiellement dans les bassins marginaux en aval des zones 
de subduction oceanique. Les sequences volcaniques et sedimentaires 
provenant de ces bassins marginaux, ainsi que les plateaux ooceaniques, 
sont arrimes aux continents le long des marges continentales actives telles 
les Cordilleres Nord et Sud Americaines. II y a aussi apport de croute 
continentale par l’activite magmatique de la subduction. 

Cette croissance de la croute continentale est peut-etre compensee 
en partie par la disparition de sediments emportes dans les zones de 
subduction ainsi que par la subduction possible de croute continentale 
dans les zones de collision continentale. 


11.1 Subduction oceanique 

Ce sont les anomalies de gravite et la sismicite qui ont attire l’attention 
sur l’importance des fosses profonds. La destruction du plancher oceanique 
dans les zones de subduction fut demontree par Isacks et al. (1968). De 
nombreuses questions restent non resolues: le mecanisme qui declenche 
la subduction, les changements de polarite et de direction de la subduc- 
tion, etc. Les zones de subduction ocean-ocean sont tres differentes des 
zones de subduction ocean-continent 

Les principales zones de subduction oceanique sont situees dans l’ouest 
du Pacifique: Tonga-Kermadec, Mer des Philippines, Japon, Indonesie. 

II y a aussi de la subduction dans l’Atlantique sud (South Sandwich) et 
dans les Caraibes. 
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Oceanic crust 


Continental 


Figure 11.1 Zone de subduction ocean-ocean 


Lithosphere 


Asthenosphere 


Lithosphere 


Oceanic-oceanic convergence 


11.1.1 Morphologie 

Les fosses profonds atteignent >10km sous le niveau de la mer. Des 
sediments sont souvent accumules dans un bassin d’avant arc (forearc 
basin) . Un arc volcanique est forme sur la plaque qui chevauche celle qui 
retourne dans le manteau. Un bassin extensionnel (back arc basin) se 
forme au dela de l’arc volcanique. 


11.1.2 Sismicite 

La sismicite se concentre le long d’une zone inclinee, dite de Wadati- 
Benioff. Les zones de Wadati-Benioff plongent avec une inclinaison generalement 
de 30 a 60 mais pouvant atteinndre 90°. Le pendage depend de la 
densite et de l’age de la plaque. Les mecanismes au foyer montrent des 
contraintes de tension jusqu’a «350km dans la plaque qui descend. La 
profondeur maximum des tremblements de terre est de 600knr mais de 
nombreuses etudes montrent que la plaque conserve son identite au dela 
de 600km. Geometrie de la plaque en subduction. 


11.1.3 Gravite 

Les zones de subduction sont caracterisees par un couple negatif-positif 
d’anomalies d’air libre. L’anomalie negative de courte longueur d’onde 
correspond au fosse profond et l’anomalie positive de grande longueur 
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Zones de convergence 1: Subduction 


d’onde est due a la densite de la plaque qui est subductee. Ces anoma- 
lies indiquent que les zones de subduction sont tres loin de l’equilibre 
isostatique. 


Les premieres experiences de tomographie sismique ont ete faites dans 
des zones de subduction. La vitesse des ondes qui se propagent dans la 
plaque qui descend est beaucoup plus elevee que celle des ondes qui tra- 
versent le manteau avoisinant. Cela suggere que la plaque reste froide et 
rigide. Cela est confirme par la plus faible attenuation des ondes dans la 
plaque que dans le manteau. La forte attenuation dans certaines regions 
du manteau suggere de la fusion partielle. Les modeles thermiques con- 
firment que la plaque qui retourne dans le manteau n’a pas le temps de 
se rechauffer. 


Le volcanisme dans bare insulaire est du essenticllement a la fusion par- 
tiellc dans le manteau. II est vraissemblable que la fusion ne se pro- 
duit pas dans la plaque qui descend mais dans le manteau au dessus 
de la plaque qui descend, suite a l’expulsion des volatiles presents. La 
temperature du solidus (temperature ou la fusion partielle apparait) du 
manteau baisse en presence de volatiles. 


Des bassins se forment par extension derriere bare volcanique dans la 
plaque qui passe par dessus. Plusieurs mecanismes ont ete proposes pour 
expliquer b extension dans ces bassins. La subduction pourrait engendrer 
une cellule de convection: bextension correspondrait a la branche mon- 
tante dans ces cellules. Alternativement, la plaque qui sombre dans le 
manteau s’enfonce rapidement et la zone de subduction recule: et ce recul 
produit des contraintes extensionnelles dans la plaque qui chevauche. 

(Mer des Philippinnes, mer du Japon, Tonga, Tasmanie, etc.) 


? ont etudie les conditions d’equilibre dans les prismes d’accretion. 


11.1.4 Structure sismique profonde 


11.1.5 Volcanisme 


11.1.6 Bassins marginaux 


11.1.7 Prismes d’accretion 
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11.1.8 Metamorphisme des zones de convergence 

La zone de subduction proprement dite produit un metamorphisme 
haute pression, basse temperature, tandis que Fare volcanique est as- 
socie a un metamorphisme de haute temperature. Les zones de subduc- 
tion produisent done une paire de ceintures metamorphiques distinctes 
(Japon, Californie, etc.) 


Le long de la cote ouest de l’Amerique, le plancher oceanique est detruit 
sous le continent. Le caractere de cette subduction varie de fa$on tres 
marquee, sans doute essentiellement a cause des differences d’age et de 
temperature de la plaque oceanique ainsi que des interactions avec des 
traces de points chauds et des fragments de dorsales inactives. 

C’est dans ces regions que de la croute continental juvenile, formee 
surtout dans les arcs volcaniques et les bassins marginaux, peut-etre 
accretee aux continents. II y a aussi addition de croute continentale par 
plutonisme et formation de grands batholites. 


II est important de noter la difference de style tectonique entre le nord et 
le sud de la Cordillere des Andes. Les mecanismes au foyer et l’activite 
neotectonique nrontrent que dans la partie Nord, le regime est essen- 
tiellement compressif, tandis que dans la partie sud, le regime est exten- 
sionnel. La zone de Wadati-Benioff sernble plus inclinee dans la partie 
sud que dans la partie nord. Cette difference de style de subduction 
s’explique peut-etre par la subduction dans la partie nord d’anciennes 
dorsales inactives provenant du point chaud des Galapagos. 

Les anomalies de gravite et les etudes sismiques indiquent que la to- 
pographic est largement compensee par une croute tres epaisse (70km). 
Get epaississement sernble du autant a la compression et au charriage en 
surface qu’a Fapport dans la croute inferieure de magmas du manteau. 

Dans la partie nord des Andes, au Perou, le batholite andin contient un 
volume considerable de roches plutoniques emplacees durant ss 80 Ma. 
Les magmas proviennent du manteau mais les signatures isotopiques 
(des rapport initiaux 87 Sr/ 86 Sr plus eleves que ceux du manteau appau- 
vri) suggerent une contamination crustale. La presence de 10 Be, de tres 


11.2 Marges continentales actives 


11.2.1 Cordillere des Andes 
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Figure 11.2 Zone de subduction ocean-continent 
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courte duree de vie (1.6Ma) et qui ne peut-etre forme que dans la haute 
atmosphere avant d’etre incorpore dans les sediments marins, montre 
que les sediments oceaniques sont en partie recycles dans les zones de 
subduction. 


11.2.2 Cordillere nord americaine 

Bien qu’il n’y a plus de subduction le long des failles de San Andreas, 
en Californie, et de la Reine Charlotte, en Colombie Britanique qui se 
sont formees suite a la disparition de la plaque de Farallon. la Cordillere 
Nord Americaine est tres semblable a la Cordillere des Andes. La croute 
en Amerique du Nord est moins epaisse que dans les Andes, a cause de 
l’extension recente du Basin & Range. Le batholite de la Sierra Nevada 
est l’equivalent du batholite andin. La presence de terrains allochtones 
est notablement plus importante dans la Cordillere Nord Americaine, 
en particulier au Canada. L’architecture de ces terrains et leur relation 
avec le socle du bouclier canadien a ete etudie par les transects South- 
ern Cordillera (Cook, F. (editor), 1995) et SNORCLE (?) de LITHO- 
PROBE. 


11.2.3 Bassins d’avant pays (Foreland basins) 

Pendant la formation des Montagnes Rocheuses, des couches de roches 
sedimentaires ont ete transportees vers Test sur la plaque nord americaine. 
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Figure 11.3 Structure profonde de la zone de subduction pres de l’ile 
de Vancouver 



Distance from trench axis (km) 

En meme temps, un bassin sedimentaire etait forme a l’est des Rocheuses, 
notament en Alberta. C’est sous le poids de la masse de roches trans- 
porters que la plaque s’est enfoncee permettant ainsi la formation d’un 
bassin (foreland basin) en avant de la chaine de montagnes Beaumont 
(1981). 


11.3 Megathrust earthquakes in the Cascadia 
subduction zone 

A large earthquake hazard is recognized for southwest British Columbia 
and adjacent United States as reflected in the seismic zoning. However, 
a poorly know threat has been from great thrust earthquakes (M>8) on 
the Cascadia subduction zone extending from southern British Columbia 
to northern California. Although common for most other subduction 
zones, no such events have occurred on this coast in the 200 year histori- 
cal period. A multidisciplinary program has now demonstrated that they 
have occurred in the past and that elastic strain is accumulating toward 
future fault rupture. Paleoseismic evidence includes abruptly submerged 
intertidal coastal marshes and turbidite layers from widespread deep sea 
landslides. The interval between past great earthquakes average about 
500 years but have been irregular; the last was 300 years ago. Repeated 
precise geodetic measurements show the shortening and buckling of the 
coastal region expected for a locked subduction thrust fault. The seis- 
mic zone is unusually narrow, limited to beneath the continental shelf, 
which reduces the hazard somewhat at inland cities. However, complete 
rupture, as suggested for the last great event, gives a magnitude of M«9. 
Hyndman and Wang (1995). 
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Figure 11.4 Seismogenic zone on the Pacific coast 
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11.4 Notes supplementaires: Flexure de la 
lithosphere. 

La lithosphere se comporte comme une plaque elastique. Cela a d’importantes 
consequences pour les mecanismes d’evolution des marges continentales 
passives, des bassins sedimentaires intracontinentaux et des zones de 
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Figure 11.5 Mechanism of stress accumulation in a subduction zone 

BETWEEN EVENTS 


EARTHQUAKE 
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subduction. La rigidite de la plaque lithospherique permet egalement 
d’expliquer la formation des bassins d’avant-pays. A grande echelle, on 
peut considerer que la lithosphere est analogue a une poutre ou a une 
barre. Si une charge est placee sur la poutre, celle-ci va s’inflechir. Cette 
inflexion peut s’observer a proximite des iles oceaniques (seamounts) 
et peut servir a determiner l’epaisseur de la lithosphere elastique. De 
meme, sous la charge de sediments deposes a la surface, la lithosphere 
flechit, ce qui tend a elargir les bassins sedimentaires. Dans le cas des 
seamounts et des bassins sedimentaires, la rigidite de la lithosphere aug- 
mente lorsque la lithosphere se refroidit. On peut egalement calculer le 
prohl d’equilibre d’une barre qui s’inflechit sous son propre poids (ana- 
logue a la lithosphere oceanique qui s’enfonce dans le manteau). Les cal- 
culs predisent que le rayon de courbure de la barre est grand (plusieurs 
centaines de km) et depend evidement de la rigidite de la barre (cad 
de la temperature). Ceci permet d’expliquer que le plancher oceanique 
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est souvent souleve de plusieurs centaines de m en amont des zones de 
subduction. 

En pratique, la rigidite de la plaque lithospherique se mesure par la 
rigidite flexurale, D: 


D = 


Eh 3 

12(1 -a 2 ) 


( 11 . 1 ) 


ou E est le module de Young de la lithosphere, a le coeficient de Poisson, 
et h l’epaisseur de la plaque elastique. La rigidite flexurale s’exprime en 
N.m , c.a.d. a les dimensions d’un moment d’une force. Elle mesure la 
resistance de la plaque aux moments des forces qui lui sont appliquees. 
Le module de Young, E, mesure la resistance a l’elongation d’une barre 
soumise a tension T : 


A l _T 
T ~ E 


( 11 . 2 ) 


Le coeficient de Poisson a mesure le rapport du racourcissement trans- 
verse a l’elongation dans une barre en tension: 


Ad 

~d 


A l T 

~ _cr T " ~ a E 


(11.3) 


Le coeficient de Poisson d’un solide varie entre 0 et 0.5. La rigidite 
flexurale de la lithosphere oceanique est de l’ordre de 10 24 — 10 25 1 V.to. 

L’echelle de longueur pour la deflection de la lithosphere est determinee 
par la longueur d’onde de la deflection, A: 


A 


2n 



(11.4) 


avec p m densite du manteau et g l’acceleration due a la gravite. Si D 
est petit, la longueur d’onde est courte et la compensation est locale 
(Airy); si D est tres grand, la longueur d’onde devient grande, il n’y a 
pas compensation. 
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Collision continentale 


12.1 Zones de convergence 2: Collision continentale 

Les Alpes, la chaine de Zagros en Iran et les Himalayas montrent differents 
stages de la colision continentale. Dans la cas des Alpes, la collision con- 
tinentale est inachevee et les continents africains et europeens ne con- 
vergent plus bien que les deux continents ne sont pas vraiment entres 
en contact. Toutefois, l’ltalie est un promontoire du continent africain 
qui a ete accrete a l’Europe. La ” suture” entre terrains europeens et 
africains se fait a la ligne insubrique a peu pres le long de la frontiere 
entre Suisse et Italie . Dans l’Himalaya, la collision de l’lnde avec l’Asie 
est profonde. L’lnde s’est enfoncee de plus de 2000km a l’interieur de 
l’Asie et la convergence continue aujourd’hui. 

La marge Sud de l’Eurasie a ete active pendant au rnoins 150 Ma avec 
la subduction de Thetys et l’accretion des terrains qui aujourd’hui for- 
rnent le plateau du Tibet. L’lnde s’est detachee du supercontinent a ca 
80Ma. Cet evenement precede l’emplacement des Deccan traps (65Ma), 
un enorme plateau de basaltes au sud ouest du bouc.lier indien. La con- 
vergence entre et l’lnde et l’Asie et la disparition de Thetys a ete tres 
rapide et le debut de la collision date de 50Ma. Depuis, la plaque in- 
dienne continue de penetrer en Asie mais a vitesse reduite de moitie. 

Les ophiolites sont exposes le long de la suture de l’lndus-Zangbo, juste 
au Nord de la haute chaine de l’Himalaya. Au Sud de la suture, les failles 
majeures (Main Frontal Thrust, Main Boundary thrust, et Main Central 
Thrust sont des failles de chevauchement inclinees vers le Nord. Au Nord 
de la suture apparaissent les grandes intrusions de granite (<30Ma). Le 
plateau du Tibet est compris entre les failles de Karakorum et d’Altyn- 
Tag, une grande faille de decrochement gauche. L’Indochine semble etre 
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Collision continentale 


Figure 12.1 Collision des plaques indiennes et arabes a la marge sud 
de l’Eurasie. 
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Figure 12.2 La collision himalayenne. Dans la partie Sud de la chaine 
himalayenne, le socle indien n’est pas deforme, tandis que le socle et 
l’ensemble du continent asiatique sont deformes en profondeur. 
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un bloc de l’Asie ecarte par l’intrusion de PInde et aujourd’hui le plateau 
du Tibet semble aussi etre expulse lateralement par la poussee de PInde. 
L’anisotropie sismique sous le plateau du Tibet est compatible avec cette 
extrusion. Les modeles reduits de la penetration d’un poingon dans une 
plaque reproduisent assez bien les styles de deformation dans le continent 
(Molnar and Tapponnier, 1975). L’epaississement de la croute fournit un 
autre source de contraintes qui pourraient expliquer les deformations en 
Asie (England and McKenzie, 1982). 
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Collision continentale 


Figure 12.3 Reconstruction des positions de la plaque indienne et 
arabe depuis 56Ma 
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Figure 12.4 Les failles principales dans l’Himalaya et le plateau du 
Tibet 



Le soulevement du plateau du Tibet a eu lieu assez recement. La to- 
mographie sismique semble indiquer un manteau chaud qui serait com- 
patible avec les modeles de type delamination du manteau lithospherique 
pour expliquer ce soulevement. Aujourd’hui, il y a de nombreux signes 
d’extension dans le plateau du Tibet. 

Les mesures geodesiques (mesures reptees de GPS et interferometrie 
a grande ligne de base) demontrent que des deformations importantes 
se produisent aujourd’hui en Chine (Zhang et ah, 2004, e.g.,). 

Les montagnes rocheuses aujourd’hui et certaines portions des Ap- 
palaches sont un exemple de ’’tectonique superficielle” ( thin skinned tec- 
tonics) oil les deformations sont confmees dans la partie chevauchante 
de la croute et laissent le socle intact. Les deformations profondes en 
Asie aujourd’hui et des orogenes tels l’orogene de Grenville demontrent 
une tectonique profonde ( thick skinned tectonics ) oil le socle est pro- 
fondement deforme. 
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Collision continentale 


Figure 12.5 Un des nombreux modeles proposes pour revolution du 
plateau tibetain 
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12.2 Himalaya. Tibetan plateau 

Some questions 

• 2 contradictory approaches in modeling stress field on large scale. 
Rigid blocks (e.g. (Molnar and Tapponnier, 1975, 1978)) vs contin- 
uum mechanics (e.g. (England and McKenzie, 1982; Houseman and 
England, 1986)). 

• Related to rheology of lithosphere (e.g. (Maggi et al., 2000)) and ther- 
mal control (e.g. (Huerta et al., 1998)). 
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• Link between uplift and exhumation on one side and structure of 
mantle revealed by seismology on the other ((Priestley et ah, 2006, 
2008)) 

• Link between erosion and uplift (e.g. (England and Molnar, 1990; 
Molnar et al., 2010)). Chicken and egg question? 
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Evolution tectonique des planetes internes 


13.1 Proprieties des planetes internes 

Les planetes internes sont Mercure, Venus, Terre, Mars. Elies ont une 
composition relativement semblable a celle de la Terre, contrairement 
aux planetes externes dont l’enveloppe gazeuse est beaucoup plus im- 
portante. Parmi les planetes internes, seule la Terre a un satellite, dont 
la masse est relativement importante. Mars a deux minuscules satellites. 
Toutes les grandes planetes externes ont des satellites. 

Les mesures astronomiques ont permis de determiner un certain nom- 
bre de proprietes de planetes internes (masse, densite, moment d’inertie). 
Parmi celles ci, la densite moyenne et le moment d’inertie sont impor- 
tants. La densite moyenne est un indicateur de la composition. Le mo- 
ment d’inertie indique si la planete est homogene ou stratifiee. 

Le moment d’inertie de Mars indique que Mars a un noyau dense. 
Le moment d’inertie de Venus est inconnu, mais il est vraissemblable 
que Venus a egalement un noyau. La Terre est la seule des planetes 
internes a avoir un champ magnetique. Ni Venus, ni Mars n’ont un champ 
magnetique aujourd’hui. Cependant, la croute martienne est magnetisee 
et il est vraissemblable que Mars a eu un champ magnetique dans le 
passe [ Connerney et al, 1999]. 


13.2 Atmosphere et Climat 

Le climat des planetes ne depend pas que de la distance au soleil et de 
la radiation solaire mais depend de fagon critique de la presence d’une 
atmosphere et de sa composition. La temperature d’equilibre avec le ray- 
onnement solaire est la temperature de corps noir. C’est la temperature 
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Table 13.1 Proprieties physiques generates 



Masse 

(kg) 

Rayon 

(km) 

Densite 
(g cm" 3 ) 

If 

g 

(ms" 2 ) 

Jour 

Annee 

Terre 

5.98 x 10 24 

6,371 

5.5 

0.335 

9.8 

1 

1 

Lune 

7.41 x 10 22 

1,738 

3.4 

0.4 

1.6 

28 


Venus 

4.87 x 10 24 

6,052 

5.3 

? & 

8.9 

243 

0.61 

Mars 

6.5 x 10 23 

3,397 

4.0 

0.365 

3.7 

1.025 

1.88 


fll s’agit du moment d’inertie reduit, cad divise par M x a 2 . (M est la masse 
et a le rayon). Pour une sphere homogene, le moment d’inertie reduit est 0.4 
J La periode de rotation de la lune est egale a la periode de son orbite autour 
de la Terre, de telle sorte que la lune presente toujours la meme face a la 

Terre 

& Parce que la surface de Venus n’est pas visible, il est difficile de bien 
mesurer sa rotation et le moment d’inertie est longtemps reste inconnu. Des 
mesures recentes de maree planetaire (i.e. deformations de la planete 
induites par le champ de gravite solaire) suggerent que Venus a un noyau 
comparable a celui de la Terre. 


Table 13.2 Conditions a la surface des planetes internes 



T surface 
(°C) 

T Corps noir 

(°C) 

P atmospherique 

variations 

temperature 

composition 

atmosphere 

Terre 

20 

-15 

0.1 MPa 

annuel, journalier 

N 2 , O 2 

Lune 

-15 

-15 

0 

journalier 

- 

Venus 

460 

15 

10 MPa 

- 

CO 2 , H 2 0 

Mars 

-55 

-60 

600 Pa 

annuel, journalier 

n 2 


moyenne a la surface de la lune qui n’a pas d’atmosphere. L’atmosphere 
legere de Mars ne modifie que tres legerement la temperature de surface. 
Par contre la Terre est nettenrent plus cliaudeparce qu’une fraction du 
rayonnenrent est piegee dans son atmosphere. Get “effet de serre” est 
absolument dramatique sur Venus qui a une atmosphere tres lourde de 
vapeur <PH 2 0 et de C0 2 ■ L’atmosphere et l’hydrosphere attenuent les 
variations journalieres et saisonnieres de la temperature a la surface des 
planetes. 
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Evolution tectonique des planetes internes 


Figure 13.1 Age apparent de la surface des planetes calcule a partir 
du nombre de crateres 

Ages of Planetary Surfaces 

Estimated from Cratering Densities 



Ga 


13.3 Age apparent de la surface 

L’age apparent de la surface des planetes est calcule a partir du nombre 
de crateres de meteorites en tenant compte des variations de flux de 
meteorites au cours du temps. Celles-ci ont ete etablies grace aux etudes 
sur la Lune qui ont montre que le flux a ete tres eleve entre 4.2 et 3.6 Ga. 
Ce sont les processus internes (essentiellemnt le volcanisme) qui forment 
les roches de surface. Les processus de surface (erosion et sedimentation) 
modifient la surface et l’age apparent. La surface de la Lune est la plus 
simple avec formation des roches de la surface entre 4.5 et 3.2 Ga sans 
recyclage ni modification de cette surface. 

Sur Terre, les 3/4 de la surface ont moins de 200 Ma. Les crateres 
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Table 13.3 Age de la surface des planetes internes 



Maximum 

Minimum 

Terre 

4 Ga 

0 

Lune 

4.5 Ga 

3.8 Ga 

Venus 

600 Ma 

400 Ma 

Mars 

4.5 Ga 

1000 Ma 


de meteorites sur les boucliers sont difficilement identifiables parce que 
la surface a ete modifiee par erosion, sedimentation et parfois par la 
tectonique. 

Sur Mars, on voit d’importantes variations de Page de la surface et 
la geologie de Mars a permis de reconnaitre des provinces volcaniques 
dont Page varie entre >4Ga et wlGa et des provinces sedimentaires. II 
y a aussi evidence qu’il y a eu de P erosion sur Mars, bien qu’aujourd’hui 
Patmosphere soit legere et que l’eau n’existe que sous forme de glace 
pres des poles. 

II faut noter que la surface n’est pas visible a cause de Patmosphere 
tres dense. Elle a ete observee depuis la Terre par radar et a partir 
de satellites dans plusieurs gammes spectrales. Sur Venus, Page de la 
surface semble uniformement jeune <600Ma et suggere que la surface a 
ete dramatiquement remplacee il ya 600Ma. Le volcanisme recent s’est 
fait sur une surface qui n’est pas recyclee [Scharber et al., 1992], 


13.4 Gravite et Topographie 

Les anomalies de gravite sur la lune et les planetes ont ete mesurees a 
partir des orbites de satellites 1 . Les mesures a partir de satellites perme- 
ttent aussi une bien meilleure resolution de la topographie, meme pour 
la face non cachee de la Lune (mission CLEMENTINE). 

La distribution de topographie sur Terre est bimodale: un mode cor- 
respond aux bassins oceaniques et l’autre aux continents. Les variations 
de l’anomalie de Bouguer s’expliquent par la compensation isostatique, 
principalement au niveau croute-manteau. 

1 NASA avait envoye un gravimetre sur la Lune qui malheureusement ne 

permettait de mesurer des valeurs aussi elevees que celles qui se trouvaient dans 
la region d’allunissage de l’instrument. 
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La topographie de la Lune est tres marquee, avec des bassins cir- 
culaires relativement plats (Marae) entoures d’anneaux escarpes. Les 
echantillons ramenes par les missions Appolo ont ete dates et des ages al- 
lant jusqu’a 4.5 Ga ont ete obtenus. De tres fortes anomalies positives de 
la gravite(>500mGal), appelees parfois " mascons” , sont associees avec 
certains bassins [Zuber et al . , 1994], La croute de la lune semble s’etre sta- 
bilised entre 3.8 et 3.2 Ga apres une periode d’intenses bombardements 
meteoritiques. La lune est le seul objet pour lequel nous avons quelques 
donnees sismiques, car des sismometres ont pu y etre deposes. Les en- 
registrements sismiques suggerent une croute basaltique d’epaisseur tres 
variable mais en moyenne assez epaisse (60-80km). 

Sur Mars, il y a un contraste entre la partie nord (plate et basse) et 
la partie sud de la planete. On peut voir la meme distribution bimodale 
que sur Terre. Le mont Olympus est un volcan de 28km de hauteur qui 
fait partie d’un tres large edifice volcanique, le bombement de Thar- 
sis (Tharsis bulge). La mission MOLA (Mars Orbiter Laser Altimetry) 
a fourni des images extraordinairement precises et detaillees de la to- 
pographie de Mars (Figure 13.2). Elle aussi permis de determiner avec 
plus de precision le champ de gravite. A partir de ces donnees combinees, 
il a ete possible de calculer l’epaisseur de la croute sur Mars (Figure ??) 
[Zuber, 2001]. 

Sur Venus, il y a quelques chaines avec topographie elevee, mais la 
topographie n’est pas bimodale. Sur toutes les planetes, la correlation 
entre la topographie et les anomalies de Bouguer est negative indiquant 
qu’il y a un certain degre de compensation isostatique. La profondeur de 
compensation isostatique semble assez grande sur Venus, ce qui suggere 
que la topographie est supportee dynamiquement. 
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Figure 13.2 Topographie des hemispheres sud et nord de Mars. La 
topographie de Mars montre une dichotomie entre la partie nord dont 
la surface est plus jeune et lisse et la partie sud beaucoup plus im- 
pactee parce que plus ancienne. 
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Figure 13.3 Carte de l’epaisseur de la croute est estimee en sup- 
posant que la gravite est produite par la topographie et les vari- 
ations d’epaisseur crustale. Les profils sont orientes S-N. On peut 
noter l’amincissement crustal sous Hellas et l’epaississement sous le 
bombement de Tharsis. 
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La Lune et Mars sont inactives aujourd’hui. II y a une activite vol- 
canique sur Venus. Les bassins circulaires sur la lune sont des bassins 
d’epanchement volcanique suivant les impacts meteoritiques. II est clair 
qu’il n’y jamais eu de Tectonique des plaques sur la Lune. 

La partie nord de Mars qui est moins elevee, pourrait representer un 
ancien bassin oceanique. De plus, Mars a eu un champ magnetique dans 
le passe [ Stevenson , 2001]. On y a reconnu des anomalies magnetiques 
lineaires. II est done possible que Mars ait connu un episode de tectonique 
des plaques qui se serait arrete lorsque la planete s’est suffisament re- 
froidie apres 3.5Ga [voir p. ex. le scenario propose par Sleep, 1994]. 

Venus qui est la planete la plus semblable a la Terre semble avoir 
subi une tectonique tout a fait differente. La planete est encore tres 
active avec de nombreux volcans, mais on n’a pas reconnu l’equivalent 
des bassins oceaniques. II semble done qu’il n’y a pas de tectonique 
des plaques sur Venus. Par contre, la surface assez jeune indique un 
renouvellement quasi total il y a 500 Ma. II y aurait done eu un episode 
de subduction generalisee avec destruction et remplacement de toute la 
croute venusienne. Par ailleurs, Venus n’a pas de champ magnetique. 
Comme le moment d’inertie de Venus est inconnu, il a ete longtemps 
impossible d’assurer que Venus a un noyau. Toutefois, les mesures des 
rnarees de la planete solide ont ete effectuees durant la mission Magellan. 
Ces mesures suggerent un noyau semblable au noyau terrestre. L ’absence 
de champ magnetique n’implique pas necessairement que ce noyau est 
solide, mais peut-etre que le noyau interne n’a pas encore commence a 
se former. 


13.6 Le systeme Terre-Lune. Origine de la Lune 

Il y a un ralentissement seculaire de la rotation de la Terre. Ce ralen- 
tissement s’explique par les interactions a travers les rnarees entre la 
Terre et la Lune. La Lune induit des rnarees qui ne sont retardees par 
la friction dans les mers peu profondes. Cette difference de phase en- 
tre l’acceleration de la gravite et la maree cause un transfert de moment 
cinetique de la Terre a la Lune. La vitesse de rotation de la Terre diminue 
tandis que la Lune accelre. Il y a ainsi transfert de moment cinetique de 
la Terre a la Lune. Une partie de l’energie est de rotation est convertie 
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Figure 13.4 Images de la topographie et du champ de gravite de 
Venus. II s’agit des variations de grandes longueurs d’onde. Des im- 
ages locales plus detaillees de la topographie qui montrent de nom- 
breux volcans sont maintenant disponibles. 
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en chaleur par la friction. A mesure que la rotation de la Lune accelere, 
elle s’eloigne de la Terre. 

Ceci repose la question de l’origine de la Lune. La Lune devait etre 
beaucoup plus proche de la Terre dans le passe qu’a present, 

Trois hypothese ont ete proposees pour l’origine de la Lune.: 


• La Lune et la Tere se sont formees dans des regions distinctes du 
syteme solaire et la Lune a ete “capturee” par la Terre. 

• La Lune et la Terre se sont formees en meme temps et dans la meme 
region du systeme solaire. Le systeme binaire existe et a evolue depuis 
la formation des planetes. 

• La Terre s’est formee et la Lune a ete “arrachee” a la Terre peu apres 
sa formation. 


Le consensus semble se porter vers une combinaison des deux dernieres 
hypotheses. La Lune aurait ete arrachee a la Terre par l’impact d’un 
meteore geant (un corps de la taille de Mars) a la fin P accretion de la 
Terre alors que le manteau et le noyau etaient deja separes. Cela irnpli- 
querait que la composition de la Lune est proche de celle du manteau 
terrestre plutot que de celle de la Terre entiere (ce qui explique que 
le noyau est tres petit). La perte d’elements volatiles suivant l’impact 
explique aussi pourquoi la Lune est appauvrie en ces elements (e.g. le 
rapport K/U est moins eleve que sur Terre). 


13.7 Conclusion 

C’est le mode de refroidissement qui determine le style tectonique. Les 
petites planetes ont une production de chaleur faible et se refroidissent 
rapidement. L’activite tectonique s’arrete rapidement. C’est le cas de la 
Lune, inactive depuis 3 Ga, et de Mars, apparement inactive aujourd’lrui 
et ou l’activite est faible depuis 3Ga et la seule activite importante plus 
recente semble etre le volcanisme dans la region de Tharsis (1 Ga). 

Par contre, on comprend toujours mal pourquoi la Terre et Venus, 
qui ont une taille et une composition tres semblables ont evolue si 
differemment. 
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13.8 Lois de conservation 

Le moment d’inertie d’un corps solide par rapport a un axe est defini 
par 

r = J J J p(R)r 2 dV (13.1) 

oil r est la distance a l’axe d’inertie. (C’est done la somme des produits 
des masses par le carre de leurs distances a l’axe). 

Pour un solide en rotation, il faut que l’energie de rotation et le mo- 
ment de la quantite de mouvement (moment cinetique) soient conserves. 
Si oj est la vitesse angulaire de rotation, l’energie de rotation est egale a 
Iuj 2 / 2 ou / est le moment d’inertie par rapport a l’axe de rotation. Le 
moment cinetique est un vecteur parallele a 1’axe de rotation de longueur 
Ilo. 

Durant les recents episodes glaciaires, la redistribution de masse a la 
surface de la Terre a change le moment d’inertie de la Terre. Durant 
les glaciations, le moment d’inertie etait moindre, et la rotation plus 
rapide. A la deglaciation, l’augmentation du momemt d’inertie a cause 
une deceleration de la rotation. 
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Remarques finales 


Lors de la formulation de la theorie de la tectonique des plaques, de 
nombreux geologues, a l’instar des geologues sovietiques, ont souleve 
l’objection que la theorie s’appliquait surtout aux oceans et expliquait 
mal revolution des continents. L’objection que la tectonique des plaques 
ne fonctionnait pas au precambrien ou surtout durant l’Archeen parce 
que la Terre aurait ete plus chaude persiste toujours. 

La concept du cycle de Wilson qui remplagait et expliquait le mysterieux 
cycle geosynclinal repondait a la premiere objection. C’est aux marges 
(actives ou passives) des continents et au frontieres de plaque que se pro- 
duit la quasi totalite de l’activite geologique. Les vestiges de l’activite 
passee se trouvent aujourd’hui a l’inerieur des continents. Les terrains 
allochtones n’etaient pas prevus par la tectonique des plaques, mais leur 
existence s’explique et rentre dans le cadre de la tectonique des plaques. 
Leur decouverte le long des marges actives a perrnis d’expliquer de nom- 
breuses observations, y compris dans les provinces precambriennes. 

L’observation de zones de subduction fossiles dans le manteau a perrnis 
de demontrer que la tectonique des plaques operait a la fin de l’Archeen 
et a vraissemblablement opere pendant la plus grande partie de l’histoire 
de la Terre. La decouverte de mineraux (zircons detritiques) ages de «4.4 
Ma suggere que la croute continentale a commence ‘a se former des 
le debut de l’histoire de la Terre. Ceci est confirme par des anomalies 
d’isotopes extincts. La question reste de savoir comment les cratons se 
sont stabilises et comment le volume de croute continentale a evolue au 
cours du temps. 

Le concept de plaque a perrnis d’evoluer et decouvrir une tectonique 
est infinimement plus complexe que le modele simplifie de plaque rigide 
avec des frontieres bien marquees. Les techniques de geodesie spatiale 
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permettent de mesurer comment la deformation est distribute autour de 
frontieres diffuses. 

II reste aussi a mieux comprendre les interactions entre la surface de 
la Terre, les plaques et le manteau. Les mouvements verticaux et sans 
deformation a l’interieur des plaques ne sont pas prevus mais ne sont 
pas contradictoires avec une tectonique entrainee par la convection du 
manteau. 

Une autre question reste celle des ineractions entre la dynamique in- 
terne et revolution de la circulation oceanique et atmospherique, la com- 
position de l’atmosphere, et le climat. Quelles sont les interactions entre 
climat, tectonique et erosion? Quel est role de la disposition des con- 
tinents dans les cycles climatiques. Comment la tectonique modifie la 
composition de 1’atmosph‘ere et des oceans. 

La Terre est la seule planete ou la tectonique des plaques semble 
operer. D’autres planetes semblent avoir un regime completement different 
Dans quelle la presence de fluide et le recyclage de volatiles dans les zones 
de subduction influencent-ils la dynamique du manteau? 
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Anomalies de gravite. 


L’acceleration de la gravite est mesuree en m s~ 2 . Elle est de l’ordre de 
9.78 m s -2 a l’equateur et de 9.83 m s -2 au pole. L’unite pratique en 
geophysique est le mGal (milligal) 10 _5 m s~ 2 . 

Le champ de gravite terrestre varie en fonction de la latitude, de 
l’altitude, et de la densite des roches. Pour eliminer les effets de la lati- 
tude et de l’altitude, on reduit le champ rnesure au niveau de la mer et 
on le compare au champ de gravite theorique au niveau de la mer. 

La valeur theorique du champ de gravite au niveau de la mer est donne 
par la formule: 

go = 9.7803185(1. + 0.005278895 sin 2 A + 0.000023462 sin 4 A) (A.l) 

ou A est la latitude. La difference entre le champ reduit au niveau de la 
mer et le champ theorique est l’anomalie de gravite. II existe differentes 
fagons de reduire le champ au niveau de la mer qui permettent de definir 
V anomalie d'air libre et V anomalie de Bouguer. 

L’ anomalie d’air libre est calculee en ajoutant au champ rnesure une 
correction pour l’altitude H: 

gAL = gmes 4“ 0.308677 gtheor (A. 2) 

ou H est l’altitude en m et g est rnesure en mGal. 

L’ anomalie de Bouguer tient aussi cornpte de la masse des roches entre 
la surface de rnesure et le niveau de la mer: 

gB = gmes + 0.3083 H - 0.04191/377 - gtheor (A. 3) 

ou p est la densite moyenne des roches en Mg m~ 3 . La valeur standard 
de la densite pour la correction de Bouguer est 2.67Mg nr -3 

En general l’anomalie de Bouguer est negativement correlee avec la 
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Anomalies de gravite. 


topographie tandis que l’anomalie d’air libre est independante de la to- 
pographic. Cette observation a donne lieu au concept cl’isostasie. Le 
geo'ide est une surface equipotentielle qui correspond au niveau moyen 
de la mer. Notez que le champ de gravite varie sur cette surface. 
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Un solide elastique ne peut se deformer de faqon permanente. La deformation 
est proportionnelle a la contrainte et le solide retrouve son etat initial 
en l’absence de contrainte. Par contre un fluide se deforme sous la con- 
trainte et il ne retrouvera pas son etat initial quand la contrainte n’est 
plus appliquee. 


A. 1.1 Loi de Byerlee 

La seule deformation permanente pour un solide est la fracture (regime 
cassant, brittle). Un materiel se fracture lorsque la contrainte de ci- 
saillement sur un plan de fracture excede la friction le long du plan. 

Les experiences ont montre que cette lirnite ne depend que de la pres- 
sion (contrainte normale) et est independante de la composition des 
roches. La loi de Byerlee donne la valeur de la contrainte de cisaillement 
necessaire a s pour qu’il y ait mouvement sur une fracture: 

a s = a(a n - P/) (A. 4) 

ou a n est la contrainte normale, P/ est la pression de fluide et 

a « 0.8 a n < 200 MPa 
a « 0.6 <r n > 200 MPa 

Le gradient de pression lithostatique est a peu pres 30MPakm _1 . La 
pression de fluide est souvent egale a la pression hydrostatique (lOMPa km” 1 ). 


A. 1.2 Deformation continue des roches 

Dans un fluide, la deformation e augmente tant que la contrainte a est 
appliquee. Le cas le plus simple est le cas du fluide newtonien pour lequel: 

(j = ve (A. 5) 

Le taux de deformation est done proportionnel a la contrainte. La con- 
stante de proportionnalite v est la viscosite du fluide. 

En fait, le taux de deformation e depend de la composition, de la 
temperature, et de la contrainte a ; en premiere approximation, l’effet 
de la pression est souvent neglige. Les lois de deformation sont done du 
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type: 


e = Aa n exp (-H/RT) (A.6) 

/A 1 /" e(l/n-l) 

cr=(^-J exp (H/nRT)= ^ 1/n eexp {H/nRT) (A.7) 

oil A, n, et l’enthalpie d’activation H sont des constantes dependant de 
la composition. R = 8.3144J mol -1 K -1 est la constante des gaz, T est 
la temperature thermodynamique. Dans la croute n « 3 tandis que dans 
le manteau n = 1 (fluide newtonien). 
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A. 2 Unites 

Les grandeurs physiques sont mesurees dans le systeme international 
d ’unites (SI). Les trois unites fondanrentales sont l’unite de longueur 
(m), de masse (kg) et de temps (s). Toutes les autres unites mecaniques 
en sont derivees. Pour completer le systeme d’unite, il faut aussi in- 
clure l’unite thermodynamique de temperature K et Punite de courant 
electrique A. Les unites les plus importantes en geologie sont: 

• Vitesse: m s -1 

• Acceleration: nr s -2 . L’unite pratique pour mesurer les variations de 
P acceleration due a la gravite est le nr gal = 10 -5 to s -2 . On utilise 
parfois l’unite de gravite g .u. = 10 -6 to s -2 . 

• Force N (N = kg x m s ~ 2 ). 

• Pression = (Pa = N m -2 ). Les contraintes sont mesurees err Pa. 
Le gradient de pression dans la croute est ~ 30 MPa km -1 . II est 
reconrmande de renrplacer les unites archaiques telles que le k&ar = 
lOOMPa = 10 8 Pa. 

• Travail, energie, clraleur: J = N x m. II est recommande de remplacer 
les unites archaiques telles 1 cal = 4.18J 

• Puissance: W = J s” 1 

• Flux d’energie: W m . On voit parfois les anciennes unites pour 
mesurer le flux de chaleur (h./.it. = fi cal cm~ 2 = 41.8mFF to -2 ). 

Certaines relations entre unites pernrettent de mieux saisir ce qu’elles 
representent. Ainsi, on peut constater que l’energie potentielle k g x 
to s -2 x m et Penergie cinetique kg x (to.s - 1 ) 2 ont bien les memes 
unites. De nreme, la pression (IPa = N to -2 = J to - 3 ) represente une 
densite d’energie (d’orr la relation thermodynamique d W = — PdV ). 

Les forces tectoniques aux frontieres de plaques qui sont souvent ex- 
primees en N to -1 = J to - 2 representent des differences d’energie po- 
tentielle par unite de surface. 
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